MICROCOPY  RESOLUTION  TEST  CHART 

NAIIONAl  BUREAU  Of  SIANUARDb  1963  A 


DOC  FILE  COPY  J&DAO  66038 


NOVA-2S.  A STIFFENED  PANEL  EXTENSION  OF 
THE  NOVA-2  COMPUTER  PROGRAM 

Lawrence  J.  Mente 
William  N.  Lee 

Kaman  AviDyne 
Burlington,  MA  01803 


December  1978 


Final  Report 


Approved  for  public  release;  distribution  unlimited. 


I 

I 


AIR  FORCE  WEAPONS  LABORATORY 
Air  Force  Systems  Command 
Kirtland  Air  Force  Base,  NM  87117 


79  02  23  087 


AFWL-TR-/8-1821* 


This  final  report  was  prepared  by  Kaman  AviDyne,  a division  of  Kaman  Sciences 
Corporation,  Burlington,  Massachusetts,  under  Contract  F29601 -78-C-001 9,  Job 
Order  88090346  with  the  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base, 

New  Mexico.  Mr.  Gerald  M.  Campbell  ( DYV ) was  the  Laboratory  Project  Officer-in- 
Charge. 

When  US  Government  drawings,  specifications,  or  other  data  are  used  for  any 
purpose  other  than  a definitely  related  Government  procurement  operation,  the 
Government  thereby  incurs  no  responsibility  nor  any  obligation  whatsoever,  and 
the  fact  that  the  Government  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be  regarded 
by  implication  or  otherwise,  as  in  any  manner  licensing  the  holder  or  any  other 
person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use, 
or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  has  been  authored  by  a contractor  of  the  United  States  Government. 
Accordingly,  the  United  States  Government  retains  a nonexclusive,  royalty-free 
license  to  publish  or  reproduce  the  material  contained  herein,  or  allow  others 
to  do  so,  for  the  United  States  Government  purposes. 

This  report  has  been  reviewed  by  the  Information  Office  (01)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS, 
it  will  be  available  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


I 


O..Q  Vv\ 


GERALD  M.  CAMPBELL 
Project  Officer 


Major,  USAF 

Chief,  Environment  and  Effects  Branch 


FOR  THE  COMMANDER 


THOMAS  U.  CIAMBRONE 

Lt  Colonel,  USAF 

Chief,  Applied  Physics  Division 


DO  NOT  RETURN  THIS  COPY.  RETAIN  OR  DESTROY. 


SECURITY  CLASSIFICATION  OF  This  PACE  (Whan  Dim  Enfrud) 


REPORT  DOCUMENTATION  PAGE 


1 REPORT  NUM8C 


2.  GOVT  ACCESSION  NO. 


READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

s.  RECIPIENT'S  CATALOG  NUMBER 


AFWL-TR-78-1 82 


1 4.  TITLE  (m\d  Submit) 


RT  » PERIOD  COVERED 


I NC1VA-2S.  A STIFFENED  PANEL  EXTENSION  OF  / M'J  Final  depart  / 

tBT  NOVA-T COMPUTER  PffOGRAff^  / JfunffrAc: 

^ flljAia.-TR.l53 

7 Au  THORfS-"*"'  — ^ C9MPRARP  OR  OptwT  NUMSERfJ 

i^.awrence  J.  Mente  William  nJ  Lee  7 ® '29601 -78-C-001  9 ^ 

9.  PERFORMING  ORGANIZATION  NAMK  ANO  ADDRESS  10.  PROGRAM  ELEMENT,  PROJ ECT.  T 

— j AREA  4 WORK  UNIT  NUMBERS 

SM-TXl.  NSW  Q1 fiJMLL  ®260l  F/88090346  (gj 

II.  CONTROLLING  OFFICE  NAME  ANO  AOORESS  Til  REPORT  BRTE  — y 

Air  Force  Weapons  Laboratory  (DYV)  (' ' ) — Dec ftobay— 1-978  ^ 

Klrtland  Air  Force  Base,  NM  87117  ,,'1  57"'cn  Jr 

14.  MONITORING  AGENCY  NAME  A AOORESSf//  dllUrmu  It on  Coatrolllnd  0111  cl)  IS.  SECURITY  CLASS,  (ot  l hit  r»po rt) 

(fjTXFWLji  1 BTF  / UNCLASSIFIED 


^HnalM^epwrt, 
lKA-TR-1 53  ^ 


NUMOERf  4J 


F29601 -78-C-001 9 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A WORK  UNIT  NUMBERS 

62601  F/88090346  (T*fj 


UNCLASSIFIED 


15«.  0 EC  LAS  SI  PI  CATION/  DOWNGRADING 
SCHEDULE 


I is.  DISTRIBUTION  STATEMENT  (ol  thlo  Report) 


Approved  for  public  release;  distribution  unlimited.  (jTjl-i  t f» 


17.  DISTRIBUTION  STATEMENT  (ot  the  mbotroct  ontorod  In  Btock  20,  It  different  from  Report) 

lil 


IB.  SUPRLEMENTARY  notes 


1 19  KEY  WOROS  (Continue  on  rover  oo  aide  II  no ceeemry  end  identity  by  block  number) 


\ 


Aircraft  Vulnerability 
Digital  Computer  Program 
Structural  Response  Analysis 
Stiffened  Panels 


Elastic-Plastic  Material 
Nuclear  Blast 
Overpressure  Effects 


1 20.  ABSTRACT  ( Continue  on  rev 


••do  II  nocoooory  and  Idontlly  by  block  number) 


p*A  stiffened  panel  extension  is  developed  into  the  N0VA-2S  and  N0VA-2LTS 

computer  programs  which  replace  the  ^rwfws^NOVA-2  and  N0VA-2LT  programs  for 
nuclear  overpressure  vulnerability  and  analysis  of  aircraft.  These  new 
versions  are  capable  of  analyzing  a stiffened  flat  or  cylindrical  panel  in 
both  the  elastic  and  inelastic  response  regions  as  well  as  retaining  the 
capability  of  analyzing  stiffeners  and  pure  panels,  individually.  In  the 
stiffened  panel  analysis  the  stiffeners  are  treated  discretely  in  either  or  

DD  .ET*  w*fr  «0'T'0N  OF  I NOV  .*  II  OMOLBTE  UNCLASSIFIED 

* E CuRitv  Classification  of  This  »*ge  Dm,.  Em ...d; 


79  02  2 3 087  Jj- 


no 


! 


SECURITY  CLASSIFICATION  OF  THIS  RAOEfUTmi  Dmtm  Bntmtmd) 

BLOCK  20.  ABSTRACT  (cont'd) 

■z?  both  coordinate  directions  and  are  allowed  various  eccentric  positions 
relative  to  the  single-layered,  multilayered  and  honeycomb  panel  skin 
configurations.  The  stiffened  panel  analysis  is  evaluated  successfully  by 
comparing  the  analytical  solutions  with  experimental  results  from  -tiiree  -■? 
— separate  test  programs  that  measured  strains  and  displacement  responses  on 
stiffened  panels.  An  evaluation  of  the  stiffened  panel  analysis  versus  the 
approach  of  analyzing  individual  elements  of  the  stiffened  panel  is  made  by 
comparing  the  two  approaches  on  the  basis  of  response  quantities  at  constant 
range  and  slant  ranges  at  constant  damage  level.  Three  stiffened  panels  that 
are  very  similar  to  those  found  on  the  B-52  a1rcraft'T>>e-usod^for--*h4-s-evaluaj 
Tt  is  concluded  from  thfs-  avaluation^that,  in  general,  the  individual 
element  approach  is  not  adequate  and  the  stiffened  panel  analysis  should  be 
used  for  these  types  of  stiffened  panels.  * 


/»  HC ) I i V 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  of  this  FAGEriFh#n  Dmtm  Enfrmd) 


PREFACE 


I 


This  report  represents  continuation  of  work  performed  by  Kaman 
AviDyne,  Burlington,  Massachusetts  and  previously  documented  in  AFWL 
TR-75-262 . The  current  report  contains  a complete  description  of  new 
extensions  and  modifications  of  NOVA-2.  Mr.  Gerald  Campbell  was  project 
officer  for  AFWL  and  Mr.  Lawrence  J.  Mente  was  project  leader  for  Kaman 
AviDyne.  This  work  was  performed  under  Contract  No.  F29601-78-C-0019 
in  the  Structural  Mechanics  Section  of  Kaman  AviDyne  headed  by 
Mr.  Emanuel  S.  Criscione. 


CONTENTS 


Section  Page 

I INTRODUCTION  9 

II  STIFFENED  PANEL  THEORETICAL  FORMULATION  11 

Elastic  Relations  foT  Stiffeners  15 

Elastic-Plastic  Relations  for  Stiffeners  17 

Inertial  Coupling  Matrices  20 

Equations  of  Motion  for  Stiffened  Panels  21 

III  COMPARISON  OF  STIFFENED  PANEL  SOLUTIONS  WITH  EXISTING 

EXPERIMENTAL  RESULTS  25 

Analytical  and  Experimental  Comparisons  for  the 
MIT  Cylindrical  Stiffened  Panel  26 

Analytical  and  Experimental  Comparisons  for 
STRESNO  Test  Specimen  Number  10  35 

Analytical  and  Experimental  Comparisons  for  the 
Stiffened  Fin  Panel  of  the  A-4C  Aircraft  49 

Conclusions  54 


IV 

V 


EVALUATION  OF  THE  STIFFENED  PANEL  ANALYSIS  VERSUS 


INDIVIDUAL  COMPONENT  ANALYSIS  61 

COMPUTER  PROGRAM  DESCRIPTION  86 

Subroutines  and  Common  Blocks  88 

Maximum  Program  Dimensions  93 

Program  Input  93 

Program  Operation  122 

N0VA-2LTS  125 

Example  Problem  129 

REFERENCES  135 

APPENDIX  1 17 


3 


ILLUSTRATIONS 


Figure 

Page 

1 

Stiffened  Panel  Geometry 

13 

2 

Stiffener  Configurations 

14 

3 

Example  Segmentation  of  Stiffener 

18 

4 

Sandwich  (Honeycomb)  Cross  Section 

23 

5 

Schematic  of  the  MIT  Integrally  Stiffened 

Clamped  Cylindrical  Panel 

27 

6 

MIT  Cylindrical  Stiffened  Panel  Model 

29 

7 

Strain  Time  Histories  for  the  MIT  Cylindrical 
Stiffened  Panel  at  Positions  1 and  2 

32 

8 

Strain  Time  Histories  for  the  MIT  Cylindrical 
Stiffened  Panel  at  Positions  3 and  4 

33 

9 

Displacement  Time  Histories  for  the  MIT  Cylindrical 
Stiffened  Panel 

34 

10 

Flat  Stiffened  Panel  Model  for  STRESNO  Test 

Specimen  No.  10 

36 

11 

Reflected  Pressure  Time  Histories  for  Shots  4 and  5 
on  Specimen  No.  10 

37 

12 

Center  Strain  Time  Histories  on  Lower  Flange  of 
Stiffener  for  Shot  No.  4 

41 

13 

Center  Strain  Time  Histories  on  Upper  Panel  Surface 
for  Shot  No.  4 

42 

14 

Center  Strain  Time  Histories  on  Edge  of  Panel  Skin 
for  Shot  No.  4 

43 

15 

Center  Displacement  Time  Histories  of  the  Stiffened 
Panel  for  Shot  No.  4 

44 

16 

Center  Strain  Time  Histories  on  Lower  Flange  of 
Stiffener  for  Shot  No.  5 

45 

17 

Center  Strain  Time  Histories  on  Upper  Panel  Surface 
for  Shot  No.  5 

46 

4 


ILLUSTRATIONS  (Continued) 


Figure  Page 

18  Center  Strain  Time  Histories  on  Edge  of 

Panel  Skin  for  Shot  No.  5 47 

19  Center  Displacement  Time  Histories  of  the 

Stiffened  Panel  for  Shot  No.  5 48 

20  A-4C  Aircraft  Fin  Stiffened  Panel  Model  51 

21  Reflected  Pressure  Time  History  for  Pressure 

Gage  No.  2 53 

22  Edge  Strain  Time  Histories  of  the  A-4C  Aircraft 

Fin  Stiffened  Panel  55 

23  Strain  Time  Histories  0.75  in  from  Edge  of  the 

A-4C  Aircraft  Fin  Stiffened  Panel  56 

24  Center  Strain  Time  Histories  in  Panel  No.  2 of 

the  A-4C  Aircraft  Fin  Stiffened  Panel  57 

25  Center  Inner  Surface  Strain  Time  Histories  in 

Panel  No.  3 of  the  A-4C  Aircraft  Fin  Stiffened  Panel  58 

26  Center  Outer  Surface  Strain  Time  Histories  in 

Panel  No.  3 of  the  A-4C  Aircraft  Fin  Stiffened  Panel  59 

27  Vertical  Fin  Stiffened  Panel  Model  64 

28  Aft  Fuselage  Stiffened  Panel  Model  65 

29  Upper  Wing  Stiffened  Panel  Model  66 

30  Comparison  of  Center  Displacement  Response  on  the 
Stiffener  of  the  Vertical  Fin  Stiffened  Panel 

(Elastic  Solution)  73 

31  Comparison  of  End  Stress  Response  on  the  Stiffener 

of  the  Vertical  Fin  Stiffened  Panel  (Elastic  Solution)  74 

32  Comparison  of  Outer  Surface  Edge  Stress  Response  on 

the  Panel  Skin  of  the  Vertical  Fin  Stiffened  Panel 
(Elastic  Solution)  75 

33  Comparison  of  Center  Displacement  Response  on  the 

Stiffener  of  the  Vertical  Fin  Stiffened  Panel 
(Elastic-Plastic  Solution)  76 


5 


ILLUSTRATIONS  (Concluded) 


Figure 

34  Comparison  of  Center  Strain  Response  on  the 
Stiffener  of  the  Vertical  Fin  Stiffened  Panel 
(Elastic-Plastic  Solution) 

35  Comparison  of  the  Edge  Strain  Response  on  the  Panel 
Skin  of  the  Vertical  Fin  Stiffened  Panel  (Elastic- 
Plastic  Solution) 

36  Comparison  of  Center  Displacement  Response  on  the 
Stiffener  of  the  Aft  Fuselage  Stiffened  Panel 

37  Comparison  of  Maximum  Stress  Response  on  the 
Stiffener  of  the  Aft  Fuselage  Stiffened  Panel 

38  Comparison  of  Outer  Surface  Center  Stress  on  the 
Panel  Skin  of  the  Aft  Fuselage  Stiffened  Panel 

39  Comparison  of  Center  Displacement  Response  on  the 
Stiffener  of  the  Upper  Wing  Stiffened  Panel 

40  Comparison  of  Inner  Flange  Strain  at  the  End  of  the 
Stiffener  of  the  Upper  Wing  Stiffened  Panel 

41  Example  of  Modal  Selection  for  a Panel  Exibiting 
Symmetry  in  both  Coordinate  Directions 

42a  Outer  Stiffener  (SIDEG(B)=+1.0) 

42b  Inner  Stiffener  (SIDEG(B)— 1.0) 

42c  Internal  Stiffener  (SIDEG(B)-+2.0) 

42d  Inner  Stiffener  with  Crimped  Sandwich  Panel 
(SIDEG(B)«+3.0) 

43  Segmentation  Tree  Structure  of  N0VA-2S 

44  Segmentation  Tree  Structure  of  Modified  N0VA-2LTS 
(N0VA-2L) 

45  Example  Problem  Input  Card  Listing 

46  Example  Problem  Output  at  1.75  Milliseconds 

47  Summary  Output  for  Example  Problem 


77 

78 

79 

80 

81 

82 

83 

109 

111 

111 

112 

112 

124 

127 

130 

132 

134 


6 


TABLES 


Page 


Geometrical  and  Material  Properties  30 

Geometrical  and  Material  Properties  38 

Pressure  Models  for  Specimen  10  40 

Geometric  and  Material  Properties  52 

Pressure  Model  Based  on  Gage  No.  2 52 

Geometric  and  Material  Properties  of  the  Vertical 

Fin  Panel  67 

Geometric  and  Material  Properties  of  the  Aft 

Fuselage  Panel  68 

Geometric  and  Material  Properties  of  the  Upper 

Wing  Panel  69 

Response  Comparison  at  Constant  Slant  Range  71 

Difference  in  Slant  Range  Between  Individual  Element 
Approach  and  Stiffened  Panel  Analysis  72 

List  of  Subprograms  of  N0VA-2S  89 

Common  Blocks  and  Subprograms  Using  Them  in 

NOVA- 2 S 90 

Dimensions  of  Variables  for  NOVA  Routine  94 

Dimension  of  DEPROB  Variables  95 

Dimensions  of  DEPROP  Variables  96 

Program  Changes  Required  by  Dimension  Changes  97 

New  DEPROP  Variables  98 

Revision  of  NOVA  Data  Group  10  103 

Revision  of  NOVA  Data  Groups  27  and  28  104 

Groups  27  and  28  Data  Options  for  Loading  Fuselage 
Elements  105 


ir 


7 


TABLES  (Concluded) 


Table 

21  DEPROP  Input 

22  Segmentation  Directives 

23  List  of  Subprograms  for 

24  Segmentation  Directives 


Page 


114 

for  N0VA-2S 

125 

N0VA-2LTS  (N0VA-2L) 

126 

for  N0VA-2LTS  (N0VA-2L) 

128 

8 


SECTION  I 
INTRODUCTION 


The  computer  code  NOVA-2  (Nuclear  Overpressure  Vulnerability 
Analysis,  Version  2)  given  in  reference  1 was  developed  to  increase  the 
level  of  sophistication  in  analyzing  nuclear  overpressure  effects  on 
aircraft.  This  program  provided  a technique  for  predicting  the  dynamic 
response  of  individual  aircraft  structural  elements,  such  as  stringers, 
frames,  and  panels,  to  the  transient  pressure  loads  associated  with  the 
blast  wave  from  a nuclear  burst.  The  NOVA-2  dynamic  response  analysis 
included  both  geometric  and  physical  nonlinearities  inherent  in  the 
behavior  of  these  structural  elements  in  the  response  range  bounded  by 
threshold  of  permanent  damage  and  catastrophic  damage.  Although  the 
NOVA-2  code  represented  a significant  improvement  over  prior  static 
solutions  coupled  with  dynamic  load  factors  to  assess  the  dynamic  response 
of  the  individual  structural  elements,  the  structural  coupling  between 
mutually  flexible  elements  was  ignored  in  analyzing  the  stiffened 
panels  of  an  aircraft.  Furthermore,  in  the  individual  element  concept 
used  in  NOVA-2,  the  pressure  loading  on  skin  panels  is  assumed  to  be 
transmitted  directly  to  adjacent  stringers  and  frames,  ignoring  the 
panel  response  effects  on  the  transmitted  loads. 

Stiffened  panels  are  those  skin  panels  of  the  aircraft  that  are 
stiffened  by  several  stringers  and/or  frames  between  the  more  rigid 
boundaries  represented  by  bulkheads,  large  longerons,  spars  and  ribs. 

To  satisfy  the  need  for  an  overall  stiffened  panel  analysis,  NOVA-2 
has  been  extended  to  include  discrete  stiffeners  within  the  cylindrical 
or  flat  panels  in  both  coordinate  directions  for  both  elastic  and 
inelastic  deformation  regions.  The  extended  code  developed  herein  is 
designated  as  N0VA-2S.  Similar  revisions  have  been  made  to  the  companion 
N0VA-2LT  code  (ref.  2),  which  provides  various  general  pressure  loading 


1.  Lee,  W.  N.,  Mente,  L.  J.,  NOVA-2  - A Digital  Computer  Program  for 
Analyzing  Overpressure  Effects  on  Aircraft,  Air  Force  Weapons 
Laboratory,  Kirtland  AFB , AFWL-TR-75-262 , Parts  1 and  2,  August,  1976. 

2.  Lee,  W.  N.,  A User's  Manual  for  N0VA-2LT,  Kaman  AviDyne,  Burlington, 
MA,  KA-TM-114 , January,  1978. 
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^ options  for  the  structural  elements;  the  extended  version  is  designated 

as  N0VA-2LTS.  The  stiffened  panel  capability  has  been  developed  within 
the  analysis  framework  that  existed  in  NOVA-2  and  N0VA-2LT  for  unstiffened 
panels,  so  that  the  past  capability  for  individual  structural  elements 
have  been  retained  in  N0VA-2S  and  N0VA-2LTS  with  the  addition  of  the 
stiffened  panel  option.  Thus,  the  stiffened  panel  response  analysis  is 
compatible  with  the  existing  blast,  aerodynamic  and  criteria  subroutines 
i of  NOVA-2 . 

Although  the  stiffened  panel  option  offers  a significant  increase 
in  the  level  of  sophistication  in  analyzing  stiffened  panel  structures, 
there  are  still  limitations  when  applied  to  some  aircraft-type  structures. 

The  boundaries  of  the  stiffened  panel  must  be  some  combination  of  clamped 
and  simply  supported.  The  geometry  of  the  stiffened  panel  is  limited  to 
flat  and  cylindrical,  although  introducing  initial  imperfections  from 
1 these  shapes  allows  approximations  for  other  shapes. 

Section  2 of  this  report  presents  the  theoretical  formulation  for 
the  stiffened  panel  analysis.  Since  the  analysis  uses  most  of  the 
unstiffened  panel  theory  given  in  NOVA-2,  only  the  additional  formulation 
for  including  stiffeners  is  presented.  To  gain  confidence  in  the 
stiffened  panel  analysis,  comparisons  between  existing  experimental 
dynamic  response  results  from  several  stif fened-panel  tests  and  corres- 
ponding analytically  determined  responses  from  N0VA-2LTS  are  presented 
| in  Section  3.  The  quality  of  the  stiffened  panel  tests  used  for  this 

comparison  ranges  from  a very  well  defined  laboratory  experiment  to  a 
field  test  with  several  uncertainties  involved.  Section  4 presents  an 
evaluation  of  the  stiffened  panel  analysis  versus  the  individual  element 
analysis  using  N0VA-2S  to  compare  the  two  approaches  based  on  response 
characteristics  and  slant  range.  This  evaluation  used  selected  stiffened 
( panels  similar  to  those  found  in  the  B-52  aircraft  and  subjected  them  to 

a simulated,  nominal  nuclear  encounter.  Response  levels  corresponding 
to  threshold  of  permanent  damage  and  catastrophic  damage  were  considered 
using  the  criteria  as  given  in  NOVA-2  and  N0VA-2S . Section  5 contains 

j 

the  computer  program  description  changes  made  to  NOVA-2  to  incorporate 
the  stiffened  panel  option.  A sample  problem  is  also  given  in  Section  5.  • 

li 
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SECTION  II 

STIFFENED  PANEL  THEORETICAL  FORMULATION 

The  addition  of  stiffeners  to  the  panel  skin  is  accomplished  within 
the  theoretical  framework  of  DEPROP  which  is  a response  routine  contained 
in  NOVA-2  (ref.  1).  The  stiffeners  in  NOVA-2S  are  treated  discretely  in 
the  analysis  and  are  not  smeared  out  over  the  stiffener  spacing.  Thus, 
the  stiffened  panel  analysis  handles  as  few  as  one  intermediate  stiffener 
or  as  many  stiffeners  spaced  over  the  panel  as  the  computer  program 
dimensions  allow.  The  stiffeners  must  be  oriented  parallel  to  either  or 
both  spatial  coordinate  directions  of  the  flat  or  cylindrical  panel,  and 
stiffener  locations  are  restricted  to  coincide  with  spatial  integration 
grid  lines.  The  stiffeners  in  the  circumferential  coordinate  direction 
can  have  variable  cross  sections. 

The  eccentricity  of  the  stiffeners  either  above  or  below  the  panel 
skin  is  taken  into  account  in  the  analysis.  Both  bending  and  membrane 
deformations  causing  normal  strains  and  stresses  in  the  stiffener's 
coordinate  direction  are  included,  but  lateral  bending  of  the  stiffener 
is  ignored.  Thus,  in  the  analysis  the  stiffener  is  assumed  symmetrical 
about  the  plane  of  bending.  The  torsional  stiffnesses  of  the  stiffeners 
are  included  in  a limited  manner  by  assuming  that  the  twisting  is  always 
elastic.  Therefore,  the  shear  stress  associated  with  torsion  of  the 
stiffener  is  assumed  small  compared  to  the  normal  stresses  and  is  neg- 
lected in  the  elastic-plastic  formulation. 

The  theoretical  development  for  the  stiffened  panel  analysis  is  an 
extension  of  the  virtual  work  theory  used  to  establish  the  equations  of 
motion  for  the  unstiffened  or  pure  panel  in  Section  4.2  of  reference  1. 
Thus,  the  internal  work  of  the  stiffeners  undergoing  infinitesimal 
virtual  displacements  is  added  to  that  of  the  skin  portion  of  the  panel 
in  the  formulation.  While  the  spatial  integration  for  the  pure  panel  is 
a surface  integral,  the  integration  for  the  stiffeners  are  line  integrals 
taken  over  the  length  of  the  stiffeners  in  the  appropriate  coordinate 
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direction.  The  position  of  the  coordinate  surface  for  the  flat  or 
cylindrical  stiffened  panel  is  defined  in  the  same  manner  as  for  the 
pure  panel  in  Section  4.2.3  of  reference  1.  This  coordinate  surface  is 
shown  in  figure  37  of  reference  1.  Thus,  the  membrane  elongation  and 
shear  strains  and  the  change  of  curvature  quantities  defined  by  the 
strain-displacement  relations  in  Section  4.2.2  of  reference  1 are  based 
on  this  coordinate  surface  for  the  stiffened  panel. 

The  geometry  of  the  stiffened  panel  is  illustrated  in  figure  1 and 
depicts  discrete  integral  stiffeners  located  along  various  integration 
grid  lines  in  both  of  the  nondimens ional  y and  3 coordinate  directions. 

Figure  2 illustrates  the  various  types  of  skin-stiffener  configurations 
that  are  treated  by  the  analysis.  The  analysis  accommodates  stiffeners 
with  any  shape  whose  cross  section  can  be  represented  as  a series  of 
connected  rectangular  segments.  Since  lateral  bending  of  the  stiffeners 
is  ignored,  stiffeners  such  as  channels  and  z-sections  are  treated 
symmetrically  as  I-sections.  Configurations  A and  B of  figure  2 show 
the  stiffener  attached  to  the  outer  and  inner  surfaces,  respectively,  of 
any  panel  skin  construction  designated  in  NOVA-2,  i.e.,  single-layered, 
multilayered  and  sandwich  (honeycomb).  Configuration  C shows  the 
stiffener  located  in  the  interior  of  a sandwich  panel  and  configuration  D 
shows  the  stiffener  attached  to  the  inner  surface  of  a sandwich  panel 
that  has  been  crimped  for  connection  purposes. 

In  the  stiffened  panel  analysis,  the  segmented  stiffener  is  treated 
as  a multilayered  configuration  with  variable  widths  for  each  segment. 

These  segments  are  referenced  to  the  selected  coordinate  surface  located 
within  the  panel  skin  in  the  same  manner  as  was  done  for  a multilayered 
skin.  The  depths  of  these  stiffener  segments,  defined  by  h^  in  figure  2, 

* are  referenced  to  the  inner  skin  surface  which  is  also  the  reference 

surface  used  for  the  multilayered  skin.  The  manrer  by  which  these 
segment  depths  are  specified  for  the  various  stf  fener  configurations 
are  discussed  and  illustrated  in  more  detail  in  Section  5.  The  provision 
for  a gap  between  the  skin  and  stiffener  is  included;  a stiffener  which 

t 

* H 
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is  not  directly  attached  to  the  skin  is  still  assumed  to  be  acting 
integrally  with  the  panel  skin.  This  gap  parameter  provides  a means  of 
defining  stiffeners  which  are  attached  to  the  stiffened  panel  through 
orthogonal  stiffeners  that  are  in  direct  contact  with  the  skin  or 
stiffeners  which  are  attached  only  to  the  interior  surface  of  the  upper 
face  sheet  of  a sandwich  panel.  The  variable  cross  section  option  in 
the  6-direction  is  accomplished  by  allowing  a varying  cross  sectional 
definition  at  each  integration  point  along  the  stiffener. 

2.1  ELASTIC  RELATIONS  FOR  STIFFENERS 


For  the  elastic  solution,  additional  membrane,  bending,  cross 
coupling  and  torsional  stiffness  coefficients  are  determined  for  the 
stiffeners  and  added  to  the  corresponding  panel  skin  stiffness  coeffi- 
cients (equation  114,  ref.  1)  at  all  spatial  integration  points  at  which 
stiffeners  are  located.  Since  the  stiffness  coefficients  for  the  panel 
skin  are  used  in  a surface  integration,  the  stiffness  coefficients  for 
the  stiffeners  which  use  a line  integration  require  a correction  factor 
when  combined  together.  When  the  center  of  gravity  of  the  stiffener  is 
above  that  of  the  skin,  the  stiffener  configuration  is  defined  as  "outer" 
and  when  below  it  is  defined  as  "inner".  Whether  the  stiffener  is 
"outer"  or  "inner"  requires  some  sign  changes  in  the  stiffness  coefficients. 
The  stiffness  coefficients  for  stiffeners  in  the  y and  g directions  are 
given  by 
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For  3-stiffeners: 
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For  both  y and  3 stiffeners: 
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where 


(2) 


(3) 


E^,  Eg  are  the  moduli  of  elasticity  for  the  y and  3 stiffeners, 
respectively 


a&o  is  the  arc  length  of  the  cylindrical  panel  and  for  a flat 
panel  is  replaced  by  width  b 

A is  the  length  of  the  panel 

Mi  N are  the  number  of  spatial  integration  points  in  the  y and 
3 directions,  respectively 


V 


are  the  weighting  values  for  Simpson's  quadrature  formula 
in  the  y and  3 directions,  respectively 


b^  is  the  width  of  the  ic^  stiffener  segment 

h^  is  the  distance  from  the  inner  panel  skin  surface  to  the 
furthest  edge  of  the  ic^  stiffener  segment 


H is  the  distance  from  the  inner  panel  skin  surface  to  the 

coordinate  surface  (defined  by  equations  115  and  116  of  ref.  1) 
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NSEG  is  the  total  number  of  stiffener  segments 

, Gg  are  the  shear  moduli  for  the  y and  B stiffeners,  respectively 

J , Jg  are  the  torsional  constants  for  the  y and  B stiffeners, 
respectively 

Where  double  signs  (+  or  +)  are  indicated  in  equations  1 and  2,  the 
upper  sign  is  used  for  "outer"  stiffeners  and  the  lower  sign  is  used  for 


"inner"  stiffeners.  The 

strains 

and 

stresses  in 

the  stiffeners 

determined  from 
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(6-stiffener) 

where  z is  the  distance  from  the  coordinate  surface  to  a designated 
position  on  the  stiffener. 

2.2  ELASTIC-PLASTIC  RELATIONS  FOR  STIFFENERS 

The  elastic-plastic  solution  applies  to  stiffened  single-layered 
and  sandwich  skin  panels.  In  NOVA-2  the  sandwich  panel  was  approximated 
by  an  equivalent  single-layered  panel,  but  in  this  current  version  the 
inelastic  response  is  determined  directly  using  the  two  thin  face  sheets 
of  the  sandwich  panel.  For  the  elastic-plastic  solution,  the  stiffeners 
are  divided  into  a sufficient  number  of  segments  so  that  the  stress 
distribution  across  the  cross  section  is  accurately  represented.  In 
N0VA-2S  and  N0VA-2LTS  the  selection  of  this  segmentation  of  the  stiffeners 
is  left  to  the  discretion  of  the  user.  Figure  3 gives  an  illustration 
of  the  segmentation  used  for  elastic  and  elastic-plastic  solution  of  a 
typical  stiffener.  The  elastic  solution  required  only  three  segments 
for  this  stiffener  while  for  the  elastic-plastic  solution  six  segments 
are  selected  to  give  a reasonable  representation  of  the  stress  distribution. 
The  constitutive  relations  for  the  stiffener's  material  are  based  on 
those  used  for  the  skin  panel  in  reference  1,  except  they  are  reduced  to 
the  uni-axial  case. 
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Elastic  Model 
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indicates  segment  number 


Figure  3.  Example  Segmentation  of  Stiffener 
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For  the  stiffeners  the  secant  modulus  is  defined  by 
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where 

aQ,  are  the  yield  stress  and  strain  of  the  stiffener,  respectively, 

and  a = Ee 
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Et  is  the  strain  hardening  slope  of  the  stiffener 
For  stiffeners  in  the  y-direction, 
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and  for  stiffeners  in  the  6-direction, 
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Sij’  ®ij  are  de^ined  by  equation  110  in  reference  1 
The  general  stress-strain  relations  for  stiffeners  are  given  by 
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The  initial  elastic,  initial  plastic,  elastic  unloading  and  reyielding 
regions  of  response  are  defined  the  same  as  given  in  equation  113  of 
reference  1.  The  integrand  quantities  for  the  stiffener  are  given  by 
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where  z is  the  distance  from  the  coordinate  surface  to  the  center  of 
the  itr‘  stiffener  segment  and  is  expressed  as  z^  - + — (h^+h^  ^)-H 

where  the  plus  sign  is  used  for  "outer"  stiffeners  and  the  minus 
sign  for  "inner"  stiffeners.  The  trapezoidal  rule  is  used  for  the 
numerical  integration  through  the  depth  of  a stiffener  in  the  equation 
of  motion. 

I 

2.3  INERTIAL  COUPLING  MATRICES 


In  the  spatial  surface  integration  of  the  kinetic  energy  the  addition 
of  the  line  integrals  in  the  y and  8 directions  to  include  the  mass  of 
the  stiffeners  leads  to  inertial  coupling  of  the  modes.  The  M 

PQ 

coefficients  associated  with  the  w-equations  of  motion  of  the  inertial 
coupling  matrix  [MJ  are  determined  from 
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where 

pq  extends  over  all  the  modal  combinations  selected  for 
the  solution 

p is  the  mass  density  of  the  panel  skin 
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p1  is  the  mass  density  of  the  it^1  stiffener 
s 

is  the  area  of  the  ic^  stiffener 
s 

h is  the  total  thickness  of  the  panel  skin 


6 , 6 

mr  ns 


are  Kronecker  deltas 


k^  = /2k^,  where  k^  is  defined  in  equation  120  of  reference  1 

k^  » /2kg,  where  kg  is  defined  in  equation  120  of  reference  1 

r,s  are  particular  values  of  m and  n,  respectively 

^m’  ^n  are  8iven  in  equations  118  and  119  of  reference  1 

Y,,  8,  are  defined  in  equation  124  of  reference  1 
J k 

NSG,  NSB  are  the  number  of  y and  8 stiffeners,  respectively 
In  general  matrix  form  the  w-equations  of  motion  are  given  by 


t2wk„!  - - 


For  the  solution  of  these  equations  in  N0VA-2S,  equation  11  is  placed 
in  the  form 


{"r.}  ' ' 72  M'1  KJ 


It  should  be  noted  from  equation  12  that  the  inertial  coupling  matrix 
has  been  inverted.  In  order  to  accomplish  this  operation,  a matrix 
inversion  subroutine  has  been  placed  in  the  new  N0VA-2S  and  N0VA-2LTS 
programs.  Although  the  above  derivation  is  only  demonstrated  for  the 
normal  motion  of  the  stiffened  panel,  similar  inertial  coupling  matrices 
have  been  established  in  the  program  for  the  inplane  motions  of  the 
panel  (u  and  v-equations  of  motion) . 

2.4  EQUATIONS  OF  MOTION  FOR  STIFFENED  PANELS 

The  Inclusion  of  the  stiffeners  necessitated  modifications  of  the 
equations  of  motion  for  the  pure  panel  in  reference  1.  These  modifications 
required  the  line  integrals  from  the  stiffeners  be  integrated  with  the 
surface  integral  of  the  panel  skin.  For  the  elastic  case  the  form  of 
the  equations  of  motion  were  not  altered  since  the  stiffness  coefficients 
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of  the  stiffeners  were  integrated  directly  with  the  stiffness  coefficients 
of  the  panel  skin.  However,  for  the  elastic-plastic  solution  of  a 
stiffened  single- layered  panel  the  equations  of  motion  given  in  equation  124 
of  reference  1 are  modified  into  the  form: 
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where  the  nomenclature  has  been  given  in  equation  124  of  reference  1 
and  in  Section  2 of  this  report. 

In  NOVA-2  the  elastic-plastic  response  of  sandwich  (honeycomb) 
panels  was  handled  by  an  equivalent  single-layered  panel.  It  was 
apparent  that  the  same  technique  used  to  include  stiffeners  into  the 
elastic-plastic  solution  could  be  used  to  solve  the  inelastic  response 
of  the  honeycomb  panels  without  reducing  the  three-layered  panel  section 
to  an  equivalent  single-layered  panel.  It  is  assumed  the  core  of  the 
sandwich  or  honeycomb  always  remains  undamaged  and  the  normal  stresses 
are  carried  just  by  the  face  sheets.  It  is  further  assumed  that  the 
stress  across  each  face  sheet  is  constant.  Figure  4 shows  the  nom- 
enclature for  the  sandwich  section.  In  the  equations  of  motion  given  by 
equation  13,  the  single  layered  expression  in  the  first  brackets 
(associated  with  the  first  summation  over  i)  is  replaced  for  the  sandwich 
panel  by 
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where 


are  defined  in  equation  123  of  reference  1 


|(h3  + h2)  - H 


In  stiffened  honeycomb  panels  the  honeycomb  core  is  often  crimped 

where  the  panel  skin  intersects  the  various  stiffeners  for  attachment 

purposes  (see  configuration  D of  figure  2) . It  is  assumed  that  the  core 

is  fully  removed  over  the  stiffener,  so  that  the  bending  resistance  of 

the  panel  along  the  stiffener  line  is  negliable  compared  with  that  of 

the  uncrimped  honeycomb  panel.  Therefore,  to  account  for  this  loss  of 

bending  resistance  from  the  crimped  honeycomb  panel,  F and  D stiff- 

ij  ij  b 

ness  coefficients  are  set  equal  to  zero  for  elastic  solutions  and  f ^ is 
set  equal  to  zero  for  inelastic  solutions  at  all  integration  points 
along  the  stiffeners.  Thus,  only  membrane  stresses  in  both  coordinate 
directions  are  transmitted  through  the  honeycomb  face  sheets  at  positions 
along  the  stiffener  lines  for  this  type  of  construction. 

Stresses  and  strains  in  the  stiffeners  are  computed,  for  printout, 
at  the  extreme  outer  and  inner  fibers  for  the  elastic  solutions  and  in 
the  center  of  the  first  and  last  segments  for  the  inelastic  solutions. 
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SECTION  III 


COMPARISON  OF  STIFFENED  PANEL  SOLUTIONS 
WITH  EXISTING  EXPERIMENTAL  RESULTS 

To  evaluate  the  stiffened  panel  option  contained  in  NOVA-2S  and 
NOVA-2LTS,  comparison  of  calculated  displacement  and  strain  time  histories 
are  made  with  existing  experimental  results  from  tests  performed  on 
stiffened  panels.  In  this  initial  evaluation  of  the  stiffened  panel 
program,  three  stiffened  panels  are  analyzed,  each  from  a different  test 
program.  The  sources  for  the  three  stiffened  panels  are  from  tests 
performed  by  the  MIT  Aeroelastic  and  Structures  Research  Laboratory 
(ref.  3),  Boeing-Wichita  in  their  Structural  Response  to  Simulated 
Nuclear  Overpressure  (STRESNO)  test  program  (ref.  4)  and  the  Naval 
Weapons  Evaluation  Facility  tests  on  A-4C  aircraft  in  the  DICE  THROW 
event  (ref.  5). 

These  three  test  sources  represent  a wide  range  in  the  overall 
quality  of  the  test  results  from  the  standpoint  of  definitions  of  both 
structure  and  loading.  The  test  in  reference  3 represents  a well  controlled 
laboratory  test  in  which  the  geometry  and  boundary  conditions  of  the 
stiffened  panel  were  well  defined  and  the  implusive  loading  with  a known 
spatial  distribution  was  carefully  calibrated  for  magnitude.  The  test 
selected  from  reference  4 was  conducted  in  a large  shock  tube  in  which 
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the  pressure  was  measured  at  one  position  on  the  panel  and  the  boundary 
conditions  of  the  panel  had  some  uncertainties.  The  test  on  the  A-4C 
aircraft  was  a field  test  in  which  pressures  were  measured  outside  the 
test  panel  area  and  some  uncertainties  also  existed  in  the  geometry  of 
the  stiffened  panel. 

In  the  following  subsections,  comparisons  are  made  between  measured 
displacements  and  strain  time  histories  from  the  three  selected  stiffened 
panel  tests  and  the  corresponding  analytical  response  obtained  using  the 
N0VA-2LTS  stiffened  panel  code.  In  order  to  obtain  better  accuracy, 
the  normal  dimensions  of  N0VA-2LTS  were  expanded  to  accommodate  more 
integration  points  to  cover  the  large  areas  of  the  stiffened  test  panels. 

3.1  ANALYTICAL  AND  EXPERIMENTAL  COMPARISONS  FOR  THE  MIT  CYLINDRICAL 

STIFFENED  PANEL 

In  reference  3 well-defined  response  tests  were  performed  on 
stiffened  cylindrical  panels  with  clamped  edges.  To  assure  reliable 
structural  geometry  and  clamped-edge  boundary  conditions,  the  test 
specimens  were  machined  from  solid  blocks  of  6061-T6  aluminum.  Figure  5 
shows  a sketch  of  the  test  specimen  which  is  nominally  a 60-degree 
cylindrical  panel,  0.1  inch  thick,  6.0  inches  long  and  6.0  inches  in 
radius.  The  boundaries  of  the  specimen  are  thick  and  massive  to  simulate 
an  ideally-clamped  edge  and,  furthermore,  they  are  attached  by  bolts  to 
a thick  steel  plate.  The  integral  inner  stiffener  in  the  circumferential 
direction  is  located  in  the  center  of  the  panel  and  is  nominally  0.1  inch 
thick  and  0.4  inch  deep.  All  boundary  edges  of  the  cylindrical  panel 
and  the  curved  edges  where  the  stiffener  intersects  the  panel  were 
machined  with  1/8-inch  filets  to  reduce  the  threat  of  premature  cracking 
due  to  stress  concentrations.  All  dimensions  of  this  panel  were  care- 
fully measured  to  determine  the  actual  geometric  properties  after  fabri- 
cation. These  actual  average  dimensions  were  used  for  the  N0VA-2LTS 
analytical  model. 

The  stress-strain  curves  in  tension  and  compression  were  determined 
experimentally  for  this  particular  6061-T6  aluminum  and  the  data  are 
given  in  reference  3.  The  impulsive  loading  was  obtained  by  placing  a 
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high  explosive  (HE)  sheet  with  a foam  buffer  over  a prescribed  area  of 
the  panel.  The  impulse  imparted  to  the  panel  by  the  HE  sheet  was  care- 
fully calibrated  from  experimental  test  data.  A sufficient  impulse  was 
imparted  to  produce  significant  permanent  deformations  of  the  stiffened 
cylindrical  panel.  Strain  time  histories  were  obtained  at  four  positions 
on  the  panel  from  high  elongation  annealed  constantan  foil-type  strain 
gages. 

This  clamped  cylindrical  stiffened  panel  subjected  to  an  impulsive 
loading  was  modeled  using  the  N0VA-2LTS  program  to  predict  the  analytical 
response  for  comparison  with  the  experimental  results.  Figure  6 illustrates 
the  geometry  of  the  analytical  model  for  the  stiffened  panel.  The 
circumferential  stiffener  located  midway  along  the  length  of  the  panel 
is  rectangular  in  cross  section  and  has  1/8  inch  filets  where  it  inter- 
sects the  panel  skin.  The  shaded  area  indicates  the  portion  of  the 
panel  that  was  impulsively  loaded  by  the  HE  sheet.  The  geometrical  and 
material  properties  of  the  stiffened  panel  are  given  in  table  1.  The 
plasticity  parameters  represent  average  values  obtained  from  the  tension 
and  compression  material  test  data.  As  shown  in  figure  6 the  stiffener 
is  divided  evenly  into  five  segments  for  the  elastic-plastic  solution. 

The  width  of  the  first  segment  is  increased  to  account  for  the  area  of 
the  filets.  In  the  skin  of  the  stiffened  panel,  five  integration  points 
are  used  in  the  z coordinate  direction.  Since  the  geometry  and  spatial 
loading  distribution  are  symmetrical  in  both  spatial  coordinate  directions, 
only  one-quarter  of  the  panel  is  modeled  in  the  N0VA-2LTS  analysis  by 
using  a 21  by  21  integration  net.  For  the  elastic-plastic  response 
solution  of  the  stiffened  panel,  28  modes  are  used.  For  the  temporal 
numerical  integration  a 0.5  microsecond  time  step  is  used. 

The  magnitude,  I,  of  the  impulsive  loading  applied  to  the  panel  is 

0.162067  psi-sec.  Through  load  option  3 in  N0VA-2LTS  this  loading  is 

applied  as  a triangular  pressure  load  over  the  first  time  step.  The 

21 

peak  pressure  (p  ) is  given  by  — and  is  equal  to  648268.0  psi.  Spatially 
in  At 

on  the  quarter  panel,  the  loading  is  zero  for  x * 0 to  1.35225  inches 
and  8 • 0 to  13.4235  degrees  and  defined  by  p for  x » 1.65275  inches 
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TABLE  1 

GEOMETRICAL  AND  MATERIAL  PROPERTIES 


Length  (in) 

6.01 

Subtended  Angle  (deg) 

59.66 

Radius  to  center  of  skin  (in) 

6.002 

Thickness  of  panel  (in) 

0.099 

Depth  of  Stiffener  (in) 

0.398 

Width  of  Stiffener  (in) 

0.093 

Material 

6061-T6 

2 4 

Mass  Density  (lb-s  /in  ) 

0.25383  x 10”3 

Modulus  of  Elasticity  (psi) 

10  7 

Poisson's  Ratio 

0.33 

Yield  Stress  (psi) 

46000 

Strain  Hardening  Slope  (psi) 

6.84  x 104 

to  3.005  inches  and  0 = 16.4065  to  29.83  degrees.  Thus,  the  edges  of 
the  HE  sheet  are  smeared  out  over  two  grid  spacings  as  an  approximation 
in  the  loading  distribution. 

Comparisons  between  analytical  and  experimental  results  for  the  MIT 
stiffened  cylindrical  panel  are  made  for  the  four  measured  inner  surface 
strain  time  histories  and  two  permanent-set  displacement  measurements. 

The  approximate  spatial  locations  of  these  strain  measurements  are 
indicated  in  figure  6 by  small  circles.  For  the  strain  measurement  the 
small  straight  line  segments  indicate  whether  the  orientation  was  axially 
or  circumferentially.  Figures  7 and  8 show  the  comparisons  for  the  four 
strain  positions  where  the  solid  lines  are  analytically  determined  from 
N0VA-2LTS  and  the  dashed  lines  are  experimentally  measured.  At  position  1 
the  experimental  strain  trace  terminated  just  after  reaching  the  peak 
and  at  position  3 the  strain  trace  briefly  went  out  of  recording  range 
during  the  peak  portion  of  the  response.  The  comparisons  were  very  good 
for  the  two  larger  strain  responses  at  positions  1 and  2 given  in  figure  7. 
The  two  lower  level  strain  responses  at  positions  3 and  4 in  figure  8 
show  reasonable  comparison,  particularly  in  phasing,  but  the  analytical 
responses  are  higher  than  the  experimental.  This  stiffened  panel  under- 
went large  plastic  deformations  throughout  the  panel  skin  and  the  stiffener. 
Figure  9 illustrates  the  analytically  determined  displacement  time 
histories  at  two  positions  on  the  panel.  The  measured  permanent  set 
values  at  these  positions  are  compared  to  the  level  of  oscillation  near 
the  end  of  the  analytical  time  histories.  The  projected  analytical 
permanent  sets  are  slightly  lower  than  the  measured  values.  This  might 
be  expected  since  the  stiffener  exhibited  plastic  lateral  buckling  over 
a small  region  near  the  ends  of  the  stiffener  which  can  not  be  represented 
in  the  analytical  model.  This  plastic  lateral  buckling  occurred  approxi- 
mately between  2.5  to  9.2  degrees  from  each  end  and  would  have  the 
tendency  to  retard  the  displacement  recovery  of  the  panel.  The  analytical 
strain  results  confirmed  the  severity  of  deformations  of  the  stiffener 
in  this  region.  In  fact,  the  maximum  analytical  compressive  strain  in 
the  stiffener  occurs  at  9 degrees  at  a magnitude  of  0.19  in/in. 
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Figure  7.  Strain  Time  Histories  for  the  MIT  Cylindrical  Stiffened  Panel  at  Positions  1 and 
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placement  Time  Histories  for  tiie  MIT  Cylindrical  Stiffened  Panel 


3.2  ANALYTICAL  AND  EXPERIMENTAL  COMPARISONS  FOR  STRESNO  TEST  SPECIMEN 

NUMBER  10 

The  STRESNO  test  program  of  reference  4 tested  a variety  of  unstiffened 
and  stiffened  panels  in  the  Sandia  Thunderpipe  shock  tube.  These  test 
specimens  were  specially  fabricated  and  attached  to  a support  frame  in  a 
manner  to  simulate  either  clamped  or  simply  supported  boundary  conditions. 
However,  it  is  believed  that  the  support  frames  had  some  compliance  in 
the  inplane  direction  and  did  not  provide  boundaries  fixed  from  inplane 
movement.  Thus,  it  was  expected  that  the  measured  normal  displacements 
would  be  larger  than  those  analytically  predicted.  From  the  stiffened 
panels  tested  in  reference  4,  specimen  number  10  was  selected  for  the 
comparison  with  the  stiffened  panel  analysis  because  this  specimen 
seems  to  have  the  best  defined  boundary  conditions  that  fit  into  the 
ideally  clamped  or  simply  supported  category.  Specimen  10  is  a flat 
36-  by  36-inch  skin  panel  stiffened  by  three  z-shaped  inner  stiffeners 
spaced  9 inches  apart.  The  ends  of  the  stiffeners  are  pinned  to  the 
support  frame.  The  skin  panel  boundaries  parallel  to  the  stiffeners  are 
hinged  while  the  skin  panel  boundaries  perpendicular  to  the  stiffeners 
are  unattached  except  for  being  riveted  to  the  stiffeners  at  their  three 
locations.  Figure  10  illustrates  the  geometry  of  the  stiffened  panel 
and  the  locations  of  the  strain  and  pressure  measurements.  The  displace- 
ment time  history  was  also  measured  at  the  center  of  the  stiffened 
panel.  The  dimensions  of  the  cross  section  of  the  stiffeners  are  also 
given  on  figure  10.  The  material  of  the  0.0625-inch  panel  skin  is 
2024-T3  aluminum  while  the  material  of  the  stiffeners  is  2024-T3511 
aluminum.  The  general  geometric  and  material  properties  are  given  in 
table  2. 

Analytical  and  experimental  comparison  were  made  for  shots  4 and  5 
on  specimen  number  10.  Shot  4 was  a purely  elastic  response  while 
shot  5 was  at  a level  of  response  in  the  threshold  of  yielding  region. 

The  outer  surface  reflected  pressure  time  histories  for  shots  4 and  5 
are  given  in  figure  11  in  which  the  approximated  pressure  time  histories 
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Figure  11.  Reflected  Pressure  Time  Histories  for  Shots  4 and 
on  Specimen  No.  10 
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Length  (in) 


Width  (in) 

Thickness  of  skin  panel  (in) 
2 4 

Mass  Density  (lb-s  /in  ) 
Modulus  of  Elasticity  (psi) 
skin  panel 
stiffeners 

Poisson's  Ratio 
Yield  Stress  (psi) 

Strain  Hardening  Slope  (psi) 


36 

0 

0 


9 

10 

0 

50000 

2 


0625 

259  x 10 

8 x 106 
8 x 106 
33 

2 x 105 


are  shown  by  solid  straight  line  segments.  Table  3 gives  the  values  of 
the  pressure  model  for  shots  4 and  5,  corrected  for  the  slight  internal 
pressure  generated  by  the  panel  response  within  the  enclosed  support 
frame  box.  The  spatial  distribution  of  the  pressure  is  assumed  to  be 
uniform  and  is  inputted  into  the  N0VA-2LTS  program  by  pressure  option  2. 

Since  the  geometry  and  spatial  loading  distribution  are  symmetrical 
in  both  spatial  coordinate  directions,  only  one-quarter  of  the  stiffened 
panel  is  modeled.  The  stiffeners  are  oriented  in  the  y-direction  and 
all  edges  of  the  stiffened  panel  are  assumed  to  be  simply  supported.  A 
15-by-23  spatial  integration  grid  and  a 4-microsecond  time  step  in  the 
temporal  integration  are  used  in  the  analysis.  For  the  elastic-plastic 
response  in  shot  5,  five  integration  points  through  the  thickness  of  the 
panel  skin  are  used  and  the  webs  of  the  stiffeners  are  divided  evenly 
into  4 segments.  For  the  analytical  solutions  of  the  two  shots,  36  modes 
were  selected  out  of  a 5 by  9 matrix  of  symmetric  modes. 

Comparisons  between  analytical  and  experimental  strain  and  displace- 
ment time  histories  are  made  for  the  various  response  locations  published 
in  reference  4.  Strain  responses  at  the  center  of  the  stiffener  are 
compared  at  the  lower  surface  of  the  lower  flange  (S10-6)  and  lower 
surface  of  the  upper  flange  (S10-4) . Strains  on  the  panel  skin  at  the 
center  of  the  stiffened  panel  are  compared  at  the  upper  skin  surface  in 
the  y-direction  (S10-2c)  and  in  the  3-direction  (S10-2a).  Normal  displace- 
ments are  compared  at  the  center  of  the  stiffened  panel.  Figures  12-15 
show  these  comparisons  for  shot  4 and  figures  16-19  for  shot  5 where  the 
solid  trace  is  the  analytical  results  and  the  dashed  trace  is  the 
measured  results. 

The  largest  strains  occurred  on  the  lower  flange  of  the  stiffener 
as  shown  in  figures  12  and  16,  respectively  for  shots  4 and  5.  For  this 
strain  the  comparison  in  magnitude  and  phasing  are  good.  The  strains  in 
the  upper  flange  of  the  stiffener  and  the  upper  surface  of  the  panel 
skin  are  shown  in  figures  13  and  17,  respectively,  for  shots  4 and  5. 

For  these  strains,  which  are  much  lower  than  those  on  the  lower  flange. 


39 


Center  Strain  (iiin/in) 


Figure  12.  Cencer  Strain  Time  Histories  on  Lower  Flange 
of  Stiffener  for  Shot  No.  4 
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the  phasing  is  good,  but  the  magnitude  of  the  measured  strain  is  approxi- 
mately double  that  of  the  analytical  strain.  This  difference  is  mis- 
leading in  that  the  distribution  of  the  strains  across  the  cross  section 
is  pretty  good  if  the  strain  distribution  is  represented  by  the  combination 
of  pure  bending  and  membrane  strains.  For  example,  if  the  extreme  inner 
and  outer  strains  are  3000  and  400  yin/in  analytically  and  3000  and 
800  yin/in  experimentally,  the  corresponding  pure  bending  and  membrane 
strains  on  the  cross  section  are  ±1700  and  1300  yin/in,  analytically  and 
±1900  and  1100  yin/in,  experimentally.  Thus,  the  overall  comparison  for 
the  strain  distribution  on  the  cross  section  is  pretty  good,  even 
though  the  smaller  analytical  and  experimented  strains  at  the  outer 
position  differ  by  a factor  of  two. 

Figures  14  and  18  show  the  comparison  of  the  edge  strains  at  the 
upper  skin  surface  of  the  local  panel  between  stiffeners,  respectively, 
for  shots  4 and  5.  The  strain  comparisons  are  good  in  peak  magnitude 
but  the  phasing  does  not  compare  as  well.  Figures  15  and  19  illustrate 
the  comparisons  of  the  center  displacement  time  histories,  respectively, 
for  shots  4 and  5.  Although  the  phasing  of  the  analytical  and  experi- 
mental displacement  responses  are  good  the  experimental  displacements 
are  larger  than  the  analytical  ones.  This  is  expected  since  it  is 
believed  the  support  system  for  the  stiffened  panel  did  not  provide 
enough  rigidity  to  prevent  inplane  movement  of  the  boundaries.  This 
inplane  movement  is  magnified  into  significant  additional  normal  dis- 
placements of  the  panel. 

3.3  ANALYTICAL  AND  EXPERIMENTAL  COMPARISONS  FOR  THE  STIFFENED  FIN  PANEL 

OF  THE  A-4C  AIRCRAFT 

An  instrumented  A-4C  aircraft  was  blast  tested  in  the  DICE  THROW 
project  at  an  approximate  free-field  overpressure  level  of  6 psi.  The 
results  of  this  test  are  reported  in  reference  5.  Two  areas  of  the 
vehicle  were  instrumented  where  the  construction  was  of  the  stiffened 
panel  type,  namely,  a nearly  flat  stiffened  panel  on  the  vertical  fin 
and  a curved  stiffened  panel  on  the  aft  fuselage.  The  curved  fuselage 
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panel,  which  contains  the  upper  longerons  as  stiffeners,  was  eliminated 
from  consideration  due  to  uncertainties  in  its  initial  geometry  and 
probable  presence  of  coupling  between  the  direct  overpressure  loading 
and  the  overall  bending  of  the  fuselage.  The  assumed  flat  panel  on  the 
vertical  fin,  which  contains  two  intermediate  stiffeners  separating 
three  panel  bays,  was  the  better  defined  structure  from  which  to  estab- 
lish an  analytical  model. 

The  boundaries  of  this  stiffened  panel  are  supported  by  ribs  and 
bulkheads  which  are  continuous  through  the  fin.  These  boundaries  are 
assumed  to  be  clamped  in  the  analytical  model.  The  geometry  of  the 
stiffened  fin  panel  is  shown  in  figure  20  along  with  the  approximate 
cross  sections  of  the  stiffeners.  The  actual  stiffened  panels  are  only 
approximately  rectangular  but  the  analytical  model  is  assumed  rectangular. 
The  material  of  the  skin  panel  and  stiffeners  are  7075-T6.  There  were 
three  experimental  pressure  time  histories  taken  on  the  fin  outside  this 
panel  area,  but  in  the  vicinity  of  the  stiffened  panel.  Pressure  gage 
number  2 was  selected  to  represent  the  assumed  uniform  pressure  load. 
Strains  were  measured  only  on  the  skin  panel  at  the  approximate  locations 
indicated  in  figure  20  for  panels  designated  as  2 and  3.  The  general 
geometric  and  material  properties  of  the  analytical  model  are  given  in 
table  4 and  the  pressure  model  is  given  in  table  5 as  determined  from 
the  pressure  time  history  in  figure  21.  The  structural  response  to  this 
pressure  load  remains  in  the  elastic  range.  Loading  option  2 is  used  in 
N0VA-2LTS  to  input  the  segmented  pressure  time  history. 

The  two  inner  stiffeners  are  oriented  in  the  y-direction  and  the 
actual  dimensions  of  the  panel  are  changed  slightly  to  achieve  the 
desired  stiffener  spacing  within  the  selected  integration  grid  in  the 
S-direction.  Since  there  is  only  symmetry  in  the  y-direction,  half  the 
stiffened  panel  is  modeled  for  the  analysis.  All  edges  of  the  stiffened 
panel  are  assumed  clamped.  A 15-by-38  spatial  integration  grid  and  a 
2.25-microsecond  time  step  in  the  temporal  integration  are  used  in  the 
analysis.  From  a 5 by  7 matrix  of  symmetric  modes,  34  modes  were  used 
in  the  solution. 
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TABLE  4 


GEOMETRIC  AND  MATERIAL  PROPERTIES 


Length  (in) 

15 

Width  (in) 

14.16406 

Thickness  of  Skin  Panel  (in) 

0.04 

2 4 

Mass  Density  (lb-s  /in  ) 

0.2617  x 10-3 

Modulus  of  Elasticity  (psi) 

10.5  x 106 

Poisson's  Ratio 

0.33 

Shear  Modulus  (psi) 

3.9  x 106 

Yield  Stress  (psi) 

68500  psi 

4 

Torsion  Constant  (in  ) 

Stiffener  No.  1 

0.002475 

Stiffener  No.  2 

0.00103 

TABLE  5 


PRESSURE  MODEL  BASED  ON  GAGE  NO.  2 


Time  (ms) 

Pressure  (psi) 

0.0 

16.58 

0.3623 

21.75 

0.725 

16.58 

5.0 

8.56 

7.25 

7.92 

9.05 

6.99 
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Comparisons  between  analytical  and  experimental  strain  time  histories 
are  made  for  measurements  at  gages  11,  14,  15,  18  and  19  as  shown  in 
figure  20.  All  these  strain  measurements  are  in  the  8-direction  and  are 
at  or  very  close  to  the  center  of  the  panel  in  the  y-direction.  Because 
of  the  slight  distortion  of  the  panel  dimensions  in  the  8-direction  and 
the  evenly  spaced  integration  points,  the  analytical  strain  positions  do 
not  exactly  coincide  with  the  measured  positions,  but  are  as  close  as 
possible,  within  the  analytical  geometric  restrictions  imposed  by  placing 
the  stiffeners  on  grid  lines.  Figures  22-24  illustrate  the  strain 
comparison  for  panel  2 at  the  edge  over  stiffener  2 (gage  15),  0.75  inches 
from  this  edge  (gage  14)  and  the  center  (gage  11),  respectively.  All 
these  strain  positions  are  on  the  outer  skin  surface.  Figures  25  and  26 
illustrate  the  strain  comparison  for  panel  3 at  the  center  on  the  inner 
and  outer  skin  surface  (gages  18  and  19),  respectively.  The  solid 
traces  on  these  figures  are  analytical  determined  strains  from  N0VA-2LTS 
and  the  dashed  traces  are  experimentally  measured  strains.  In  general, 
these  comparisons  are  not  as  good  as  the  previous  comparisons  in  Sections  3.1 
and  3.2.  The  analytical  results  exhibited  a higher  frequency  response 
than  the  experimental  results.  The  general  behavior  of  the  time  history 
were  similar  in  most  cases,  that  is,  whether  the  response  was  primarily 
compression,  tension  or  oscillatory  between  tension  and  compression. 

The  peak  magnitude  comparisons  for  several  of  the  plots  (figs.  22, 

23  and  26)  were  fair.  Strain  responses  are  very  sensitive  throughout 
the  panels  and  can  change  very  rapidly  over  short  distances.  Considering 
the  uncertainties  between  the  analytical  stiffened  panel  model  and  a 
field  tested  actual  aircraft  structure,  the  comparisons  are  considered 
reasonable. 

’.4  CONCLUSIONS 

The  stiffened  panel  analysis  using  N0VA-2LTS  has  been  compared  to 
experiments  on  stiffened  panels  that  varied  in  the  quality  of  the 
testing  techniques  employed.  The  N0VA-2LTS  stiffened  panel  analysis 
comparisons  with  the  well  defined  laboratory  tests  are,  generally,  very 
good.  Comparisons  are  good  for  the  lesser  controlled  large  shock  tube 
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tests  on  special  structural  models  and  are  at  best  fair  for  the  least 
controlled  field  test  of  an  actual  aircraft.  More  comparisons  between 
the  stiffened  panel  analysis  and  tests  on  actual  aircraft  stiffened 
panel  structures  are  needed  in  which  the  loading  and  structure  are  well 
defined.  Tests  recently  performed  on  the  fuselage  section  of  the 
KC-135  aircraft  and  currently  being  performed  on  the  B-52  fuselage 
section  in  the  Sandia  THUNDERPIPE  Shock  Tube  are  prime  candidates  for 
further  correlation  of  the  N0VA-2LTS  stiffened  panel  analysis. 
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SECTION  IV 

EVALUATION  OF  THE  STIFFENED  PANEL  ANALYSIS 
VERSUS  INDIVIDUAL  COMPONENT  ANALYSIS 


The  stiffened  panel  analysis  option  of  NOVA-2S  is  evaluated  relative 
to  the  analyses  of  individual  components  of  the  stiffened  panel  (stiffener 
and  panel  between  stiffeners)  using  the  beam  and  pure  panel  options  of 
NOVA-2S  which  have  been  retained  from  the  original  NOVA-2.  Three 
stiffened  panel  configurations  have  been  selected  from  the  B-52  aircraft 
structure  in  this  evaluation  and  subjected  to  the  overpressure  loading 
from  a nominal  nuclear  encounter  with  the  B-52  aircraft.  The  evaluation 
is  based  on  the  two  levels  of  damage  considered  in  N0VA-2S,  namely, 
threshold  of  permanent  damage  and  catastrophic  damage.  Analyses  are 
performed  using  the  response-only  option  in  N0VA-2S  for  the  entire 
stiffened  panel,  the  individual  stiffener,  and  the  pure  panel  between 
stiffeners  at  the  same  load  level  (equal  ranges) . The  ranges  at  which 
the  response  comparisons  are  to  be  made  are  determined  by  the  weakest 
structure  reaching  the  two  damage  levels.  For  this  response  evaluation, 
comparisons  of  the  critical  response  levels  are  made  for  the  three 
structures.  Secondly,  analyses  are  performed  using  the  iteration 
option  in  N0VA-2S  to  determine  the  critical  slant  range  at  which  thres- 
hold of  permanent  damage  and  catastrophic  damage  occurs  for  the  stiffened 
panel  and  the  governing  individual  structural  component  or  element. 
Comparisons  of  the  slant  ranges  are  made  for  this  evaluation.  The 
objective  of  the  response  and  range  evaluations  is  to  show  the  degree  of 
error  introduced  by  analyzing  the  components  of  a stiffened  panel 
individually  as  is  done  in  the  beam  and  panel  options  of  NOVA-2  rather 
than  analyzing  the  entire  stiffened  panel  by  N0VA-2S.  If  individual 
structural  components  analyses  are  acceptable  compared  with  the  complete 
stiffened  structure  analysis,  there  is,  generally,  an  advantage  of  less 
computer  cost  using  the  individual  component  analyses. 
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Three  stiffened  panels  were  selected  for  analysis  from  the  B-52 
aircraft  structure  which  exhibit  skin-stringer-frame  type  of  construction. 
These  three  stiffened  panels  were  selected  from  reference  6 and  are 
generally  described  as  follows: 

1)  a three  bay  skin-stringer  panel  in  the  vertical  fin  bounded 
between  rudder  stations  2 and  44  and  between  the  aft  auxilary 
spar  and  closure  beam  (see  figs.  19  and  20  of  ref.  6); 

2)  a twelve  bay  skin-frame  panel  in  the  aft  fuselage  bounded 
between  the  upper  and  lower  longerons  and  between  stations  1357 
and  1477  (see  figs.  8 and  11  of  ref.  6);  and 

3)  . an  eighteen  bay  skin-stringer  panel  in  the  upper  wing  surface 

bounded  between  the  front  and  rear  spars  and  between  WS402  and 
372  (see  fig.  31  of  ref.  6). 

The  purpose  of  the  models  selected  is  only  to  serve  as  example 
structures  to  evaluate  the  stiffened  panel  analysis  versus  individual 
component  analysis  and  not  to  analyze  the  vulnerability  of  the  B-52 
aircraft.  Therefore,  the  model  geometry  and  loading  distributions  will 
be  idealized  to  produce  symmetry  in  both  coordinate  directions  in  order 
to  obtain  more  accurate  solutions  for  the  stiffened  panels  with  the 
available  integration  points  and  modes  for  comparison  with  the  individual 

•A 

component  solutions. 

The  general  information  required  by  N0VA-2S  involving  the  principal 
dimensions  of  the  B-52  aircraft  are  given  in  reference  6.  Location 
dimensions  of  the  selected  stiffened  panels  are  also  given  in  reference  6. 
The  nuclear  burst  orientation  relative  to  the  aircraft  for  the  vertical 
fin  and  fuselage  panels  is  from  the  side  at  orientation  number  15  while 
for  the  wing  panel  the  orientation  is  from  above  at  orientation  number  9. 
The  pressure  loadings  on  these  panels  were  assumed  to  be  uniform. 


6.  Leang,  L.  T.  and  Swaney,  T.  G.,  Analytical  Models  for  the  B-52H, 
EC-135A  and  747-200B  Aircraft.  Air  Force  Weapons  Laboratory, 
Kirtland  AFB,  AFWL-TR-72-197 , Vol.  XI  (B-52H  Aircraft) , July  1974. 
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The  idealized  structural  models  for  the  fin,  fuselage  and  wing 
stiffened  panels  are  illustrated  in  figures  27  through  29.  The  geometry 
of  the  stiffened  panel  models  and  the  cross  section  dimensions  of  the 
stiffeners  are  given  in  these  figures.  The  general  geometric  and 
material  properties  for  the  three  stiffened  panels  are  listed  in  tables  6-8 
for  the  individual  stiffener  model,  the  pure  panel  model  and  the  entire 
stiffened  panel  model.  In  addition,  the  N0VA-2S  analysis  parameters, 
such  as,  integration  grid,  modes,  time  increment  and  number  of  masses 
are  given  in  these  tables.  For  all  elastic-plastic  panel  solutions, 
five  integration  points  are  used  through  the  thickness  of  the  skin 
panel,  and  the  webs  of  the  various  stiffeners  are  evenly  divided  into 
usually  four  segments. 

For  elastic  solutions  of  the  individual  stiffener,  pure  panel  and 
stiffened  panel,  the  boundaries  of  these  structures  are  clamped  for  the 
evaluations  and  the  range  is  keyed  on  threshold  of  permanent  damage. 

For  the  evaluation  keyed  on  catastrophic  damage  both  boundary  conditions 
of  clamped  and  simply  supported  are  used  for  all  the  structural  models. 

The  reason  for  using  both  boundary  conditions  in  the  elastic-plastic 
* solutions  is  the  uncertainty  associated  with  the  strain  criterion  at 
the  clamped  ends  for  the  beam  analysis.  In  the  beam  analysis  for  large 
inelastic  deformation  solutions,  the  special  technique  used  to  predict 
the  strain  right  at  the  boundary  discontinuity  results  in  extremely  large 
strains  being  determined  at  the  ideal  clamped  end  boundary.  The  strain 
gradient  in  this  small  local  region  near  the  end  is  extremely  steep  for 
large  inelastic  deformations.  The  practically  of  this  idealized  strain 
calculation  at  the  clamped  end  of  the  beam  analysis  for  establishing 
catastrophic  damage  for  real  aircraft  structure  is  uncertain  at  the 
present  time.  In  the  comparison  between  the  individual  stiffener 
solutions  for  catastrophic  damage  with  pure  panel  and  stiffened  panel 
solutions,  the  differences  are  distorted  by  this  strain  criterion  at  the 
clamped  end  of  the  beam.  However,  the  comparison  is  still  partly 
meaningful  because  that  is  what  is  used  in  the  NOVA-2  code.  To  obtain  a 
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Figure  29.  Upper  Wing  Stiffened  Panel  Model 
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TABLE  6 

GEOMETRIC  AND  MATERIAL  PROPERTIES  OF  THE  VERTICAL  FIN  PANEL 




Individual 

Pure 

Stiffened 

Stiffener 

Panel 

Panel 

Length  (in) 

42.0 

42.0 

42.0 

Width  (in) 

- 

10.9 

32.7 

Stiffener  Spacing  (in) 

10.9 

- 

10.9 

Thickness  of  Skin  Panel  (in) 

0.032 

0.032 

0.032 

Effective  Skin  Width  (in) 

3.2 

- 

- 

2 A 

Mass  Density  (lb-s  /in  ) 

0.259  x 10“3 

0.259  x 10~3 

0.259  x 10“3 

Torsional  Constant  (in  ) 

- 

- 

0.393  x 10-3 

Number  of  Masses 

10 

- 

- 

Integration  Grid 

- 

19  x 

19 

19  x 19 

Number  of  Modes 

- 

30 

30 

Time  Increment  (s) 

9 x 10“ 6 

2 x 

10“6 

4 x 10“6 

Skin  Panel  Material 

Stiffener  Material 

(2024-T3  AL) 

(7075 

-T6  AL) 

Modulus  of  Elasticity  (psi) 

10.5  x 106 

10.4 

x 106 

Shear  Modulus  (psi) 

4.0  x 106 

4.0 

x 106 

Poisson's  Ratio 

0.33 

0.33 

Yield  Stress  (psi) 

50000 

70500 

Strain  Hardening  Slope  (psi) 

1.24  x 105 

5.9 

x 104 

Ultimate  Strain  (in/in) 

0.15 

0.1 

TABLE  7 

GEOMETRIC  AND  MATERIAL  PROPERTIES  OF  THE  AFT  FUSELAGE  PANEL 


Length  (in) 

Subtended  Angle  (deg) 

Radius  (in) 

Stiffener  Spacing  (in) 

Thickness  of  Skin  Panel  (in) 

Effective  Skin  Width  (in) 

2 4 

Mass  Density  (lb-s  /in  ) 

4 

Torsional  Constant  (in  ) 
Number  of  Masses 
Integration  Grid 
Number  of  Modes 
Time  Increment  (s) 


Modulus  of  Elasticity  (psi) 
Shear  Modulus  (psi) 

Poisson's  Ratio 
Yield  Stress  (psi) 

Strain  Hardening  Slope  (psi) 
Ultimate  Strain  (in/in) 


Individual 

Stiffener 

Pure 

Panel 

Stiffened 

Panel 

- 

10.0 

120.0 

57.5 

57.5 

57.5 

106.5 

108.0 

108.0 

10.0 

- 

10.0 

0.064 

0.064 

0.064 

1.92 

- 

— 

0.259  x 10-3 

0.259  x 10-3 

0.259  x 10-3 

- 

- 

0.959  x 10-3 

15 

- 

- 

- 

15  x 23 

19  x 19 

- 

26 

30 

10  x 10"6 

2 x 10-6 

8 x 10-6 

Skin  Panel  and  Stiffener  Material 
(7075-T6  AL) 

10.4  x 106 

4 x 106 

0.33 

70500 

5.9  x 104 

0.1 

TABLE  8 
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GEOMETRIC  AND  MATERIAL  PROPERTIES  OF  THE  UPPER  WING  PANEL 


I 

Individual 

Stiffener 

Pure 

Panel 

r 

Stiffened 

Panel 

Length  (in) 

30.0 

30.0 

30.0 

Subtended  Angle  (deg) 

- 

1.005201 

18.093624 

Radius  (in) 

- 

404.695 

404.695 

Stiffener  Spacing  (in) 

7.1 

- 

7.1 

Thickness  of  Skin  Panel  (in) 

0.271 

0.271 

0.271 

Effective  Skin  Width  (in) 

5.6 

- 

- 

2 A 

Mass  Density  (lb-s  in  ) 

0.259  x 10-3 

0.259  x 10~3 

0.259  x 10“3 

Torsional  Constant  (in  ) 

- 

- 

0.0131 

Number  of  Masses 

15 

- 

- 

Integration  Grid 

- 

19  x 15 

15  x 28 

Number  of  Modes 

- 

28 

28 

Time  Increment  (s) 

2 x 10"6 

1.5  x 10'6 

4 x 10-6 

Modulus  of  Elasticity  (psi) 

Shear  Modulus  (psi) 

Poisson's  Ratio 

Yield  Stress  (psi) 

Strain  Hardening  Slope  (psi) 

Ultimate  Strain  (in/in) 

Skin  Panel  and  Stiffener  Material 
(7075-T6  AL) 

10.4  x 106 

4 x 106 

0.33 

70500 

5.9  x 104 

0.1 

,,  ■ . ...... 
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more  valid  comparison  which  is  not  shadowed  by  uncertainty  in  the 
critical  strain  criteria,  the  beam,  pure  panel  and  stiffened  panel 
models  also  are  used  with  simply  supported  boundary  conditions  which 
shifts  the  critical  strain  location  away  from  the  boundaries.  With  the 
simply  supported  boundaries  all  three  structural  models  used  in  the 
evaluation  are  compatible  on  a criteria  basis. 

The  results  of  this  evaluation  based  on  response  and  slant  range 
comparisons  are  given  in  tables  9 and  10,  respectively.  Table  9 gives 
elastic  and  inelastic  response  comparisons  at  constant  slant  range 
(equal  loading)  between  the  individual  structural  elements  and  the 
stiffened  panel  analysis  for  the  three  structural  configurations.  The 
evaluation  is  based  on  the  comparison  of  the  CRIT  values  determined  at 
the  range  at  which  CRIT  is  approximately  unity  for  the  weakest  structure. 
CRIT  is  the  ratio  of  the  critical  stress  or  strain  response  parameter  in 
the  structure  to  the  yield  stress  value  for  threshold  of  permanent 
damage  or  ultimate  strain  value  for  catastrophic  damage. 

In  all  cases  considered  in  this  evaluation,  the  individual  stiffener 
was  the  weakest  structure.  Therefore,  the  range  (or  loading)  which 
produced  yielding  for  the  elastic  response  or  fracturing  for  the  elastic- 
plastic  response  in  the  individual  stiffener  is  used  as  the  basis  for 
the  comparison  between  the  individual  element  analysis  approach  and  the 
stiffened  panel  analysis.  The  percentage  difference  tabulated  in 
table  9 indicates  the  error  introduced  by  using  the  individual  element 
approach  instead  of  the  more  correct  stiffened  panel  analysis. 

To  further  illustrate  the  differences  in  the  structural  response 
for  the  two  analysis  approaches  under  the  same  loading,  figures  30-40 
show  selected  comparisons  for  displacement,  stress,  and  strain  time 
histories.  For  elastic  solutions,  comparisons  were  made  for  the  center 
displacement  response  on  the  central  stiffener,  the  end  stress  or  strain 
response  of  the  central  stiffener,  and  the  stress  or  strain  response  at 
some  position  on  the  skin  panel  adjacent  to  the  central  stiffener.  It 
should  be  noted  that  skin  panel  comparisons  are  also  influenced  by 
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TABLE  9 

RESPONSE  COMPARISON  AT  CONSTANT  SLANT  RANGE 


Vertical  Fin  Panel 

Type  of  Response 

Elastic 

Elastic-Plastic 

Elastic-Plastic  j 

Boundary  Condition 

Clamped 

Clamped 

1 

Simply  Supported 

Analysis  Approach 

CRIT 

% Diff. 

CRIT 

Z Diff. 

CRIT 

% Diff. 

Individual  Elements 

18.7 

822 

227 

a)  Stiffener 

1.0 

1.005 

0.99 

b)  Panel 

0.781 

0.0486 

0.0986 

Stiffened  Panel 

0.843 

0.109 

0.31 

Aft  Fuselage  Panel 

Analysis  Approach 

CRIT 

% Diff. 

CRIT 

% Diff. 

CRIT 

% Diff. 

Individual  Elements 

51.4 

1187 

2094 

a)  Stiffener 

0.952 

1.03 

1.108 

b)  Panel 

0.794 

0.117 

0.0623 

Stiffened  Panel 

0.629 

0.08 

0.0505 

Upper  Wing  Panel 

Analysis  Approach 

CRIT 

% Diff. 

CRIT 

% Diff. 

CRIT 

% Diff. 

Individual  Elements 

-10.5 

321 

891 

a)  Stiffener 

0.98 

0.94 

1.03 

b)  Panel 

0.81 

0.275 

0.074 

Stiffened  Panel 

1.095 

0.223 

0.104 
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TABLE  10 


DIFFERENCE  IN  SLANT  RANGE  BETWEEN  INDIVIDUAL 
ELEMENT  APPROACH  AND  STIFFENED  PANEL  ANALYSIS 


Damage  Criteria 

Boundary  Condition 

TPD 

Clamped 

CD 

Clamped 

CD 

Simply  Supported 

Structural  Configuration 

Percentage  Difference 

Vertical  Fin  Panel 

15.9 

138 

32.7 

Aft  Fuselage  Panel 

24.3 

44.1 

34.4 

Upper  Wing  Panel 

-4.04 

11.1 

56.6 

TPD  denotes  threshold  of  permanent  damage 
CD  denotes  catastrophic  damage 
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Time  (ms) 


N0VA-2S  Elastic  Solutions 

Stiffened  Panel  Analysis 

Individual  Stiffener  Analysis 


Figure  31.  Comparison  of  End  Stress  Response  on  the  Stiffener 

of  the  Vertical  Fin  Stiffened  Panel  (Elastic  Solution) 


Figure  33.  Comparison  of  Center  Displacement  Response  on  the  Stiffener  of  the 
Vertical  Fin  Stiffened  Panel  (Elastic-Plastic  Solution) 


Inner  Flange  Strain  (ln/in) 


Outer  Surface  Strain  (ln/in) 
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Center  Displacement  (in) 


N0VA-2S  Elastic  Solution 


Time  (ms) 


Figure  36.  Comparison  of  Center  Displacement  Response 
on  the  Stiffener  of  the  Aft  Fuselage 
Stiffened  Panel 
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Figure  37.  Comparison  of  Maximum  Stress  Response  on  the 
Stiffener  of  the  Aft  Fuselage  Stiffened  Panel 
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9 

Figure  38.  Comparison  of  Outer  Surface  Center  Stress  K 

on  the  Panel  Skin  of  the  Aft  Fuselage 
Stiffened  Panel 
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Comparison  of  Center  Displacement  Response 
on  the  Stiffener  of  the  Upper  Wing  Stiffened 


Stress  (1000  psl) 


solution  accuracies,  since,  for  the  multi-bay  models,  only  a few  inte- 
gration points  are  available  within  each  panel  bay  compared  to  solutions 
for  individual  panels.  These  figures  show  response  comparisons  mainly 
for  the  elastic  response  of  the  fin,  aft  fuselage  and  wing  panels. 

In  the  elastic-plastic  response  comparisons,  the  differences, 
generally,  were  too  large  to  provide  meaningful  time  history  comparisons. 
Comparisons  for  elastic-plastic  response  are  shown  in  figures  33-35 
only  for  the  vertical  fin  panel  with  simply  supported  boundaries.  The 
very  large  differences  indicated  for  the  elastic-plastic  response  com- 
parisons are  somewhat  misleading,  since  the  plastic  collapse  and  in  some 
cases  the  buckling  of  the  curved  panels  accelerate  the  response  rapidly 
near  the  failure  loading. 

A better  indication  for  these  comparisons,  especially  for  catas- 
trophic damage,  are  given  in  table  10  where  differences  in  slant  range 
are  given  for  the  same  damage  level.  Table  10  shows  the  percentage 
differences  in  slant  range  between  the  individual  element  approach  and 
the  stiffened  panel  analysis  based  on  threshold  of  permanent  damage  and 
catastrophic  damage.  From  the  response  and  slant  range  comparisons 
given  in  tables  9 and  10  and  figures  30-40,  the  following  general 
observations  are  made: 

1.  There  are  significant  differences  between  the  solutions  from 
the  stiffened  panel  analysis  and  the  individual  structural 
element  analyses,  whether  the  comparisons  are  based  on  response 
parameters  or  slant  range.  These  differences  are  less  for 
elastic  response  than  inelastic  response.  The  percentage 
difference  is  reduced  when  compared  on  the  basis  of  slant 
range. 

In  all  three  stiffened  panel  configurations  the  stiffeners 
were  the  critical  structural  members  in  both  regions  of 
response. 


2. 


The  individual  stiffener  structural  model  was  always  weaker 
than  the  stiffened  panel  model,  except  for  the  elastic  solutions 
of  upper  wing  panel  configuration.  This  exception  occurred 
because  the  skin  of  the  panel  is  very  thick,  so  that  the 
computed  effective  width  of  the  skin  produced  a significant 
skin  segment  in  the  individual  stiffener  model. 

From  figures  30-38,  which  show  selected  comparisons  for  the 
fin  and  aft  fuselage  panels,  the  response  time  histories 
indicate  that  the  time  of  peak  response  is  less  for  the 
stiffened  panel  models.  Thus,  as  might  be  expected,  the 
stiffened  structured  models  are  higher  frequency  than  the 
corresponding  individual  stiffener  models. 

In  the  use  of  NOVA-2S,  the  slant  range  is  the  more  important 
parameter  on  which  to  draw  a conclusion  from  this  evaluation 
of  stiffened  panel  analysis  versus  individual  element  analysis. 
From  table  10  the  percentage  difference  in  slant  range  are 
between  4 percent  and  24  percent  for  threshold  of  permanent 
damage  and  between  11  percent  and  138  percent  for  catastrophic 
damage.  These  differences  are  significant  and  become  even 
more  significant  in  terms  of  a volume  envelope.  It  is  there- 
fore concluded  that  the  stiffened  panel  analysis  should  be 
used  instead  of  the  individual  element  technique  for  stiffened 
panels  as  described  in  this  report.  The  individual  element 
technique  is  still  useful  for  many  aircraft  structures,  such 
as  pure  panels  bounded  by  ribs,  spars  or  bulkheads,  ribs 
analyzed  for  buckling,  and  configurations  with  free  or  spring 
supported  boundary  conditions. 


85 


X 


SECTION  V 


COMPUTER  PROGRAM  DESCRIPTION 


This  section  outlines  the  changes  made  in  the  computer  program, 
including  dimension  changes  and  preparation  of  the  input  data.  The  user 
is  referred  to  reference  1 for  full  documentation  of  the  NOVA-2  computer 
code. 

The  postscript  "S"  in  N0VA-2S  refers  to  the  addition  of  discrete 
stiffeners  in  the  DEPROP  model,  as  described  earlier.  Although  the  most 
significant  change,  this  was  not  the  only  modification.  A brief  summary 
of  the  changes  follows: 

In  the  NOVA  routines,  a fuselage  loading  option  was  added  to 
permit  either  a circumferentially  uniform  or  nonuniform  blast  load  for 
beam  or  panel  elements.  Previously  only  frames  and  radome  elements 
received  the  nonuniform  load.  The  change  necessitated  modifying  the 
NOVA  input  slightly. 

In  DEPROB,  the  maximum  allowable  flanges  in  the  cross-sectional 
model  (NLK)  was  increased  from  20  to  21  in  order  to  provide  better 
representation  of  certain  elements. 

The  summary  output  for  DEPROB  and  DEPROP  iterative  runs  was  modified 
to  indicate  the  type  of  damage  corresponding  to  the  CRIT  used.  For 
example,  threshold-of-damage  criteria  for  a frame  can  be  either  tensile 
or  compressive  yielding  of  the  material,  or  a compressive  buckling  of 
the  outstanding  leg. 

Several  major  changes  were  made  in  DEPROP.  The  most  significant 
was  the  addition  of  stiffeners  in  either  coordinate  direction  in  the 
panel  model.  These  stiffeners  must  be  located  at  the  grid  lines  in  the 
spatial  integration  model,  but  can  be  located  either  on  the  inside  or 
outside  of  a multilayered  panel  and  in  the  interior  of  a sandwich  panel. 
The  cross  section  is  modelled  in  a manner  similar  to  the  DEPROB  models. 
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The  cross  section  can  vary  along  the  length  of  the  beta  (6)  stiffeners 
(stiffeners  running  parallel  to  the  beta  axis),  but  not  for  gamma  (y) 
stiffeners  (see  figure  1) . Stiffeners  can  be  of  different  construction 
and  material,  but  are  assumed  to  be  attached  at  each  grid  point  in  the 
system. 

The  inner  and  outer  flanges  are  monitored  for  maximum  CRIT  at  each 
grid  point  and  the  results  printed  out  at  the  conclusion  of  the  run. 

The  criteria  used  are  the  same  as  for  stringer  elements  in  the  DEPROB 
models. 

In  general,  the  stiffeners  introduce  coupling  in  the  mass  matrix 
and  this  capability  was  added.  It  was  made  optional,  however,  because 
the  matrix  algebra  requires  considerable  storage  and  computer  time,  and 
may  not  be  significant  for  all  problems.  By  rejecting  the  option,  only 
the  diagonal  terms  of  the  mass  matrix  are  included. 

In  addition  to  including  stiffeners,  the  modal  representation  of 
the  panel  response  was  expanded  to  include  non-symmetric  mode  shapes. 

Thus,  either  a non-symmetric  panel  or  a panel  subjected  to  a non-symmetric 
load  can  be  analyzed.  This  change  necessitated  changing  the  numbering 
system  of  the  modes,  since  the  even-numbered  modes  had  been  automatically 
excluded  in  the  old  system. 

Printout  of  the  stresses,  strains  and  displacements  at  user-selected 
spatial  locations  has  been  added,  whereas  before,  the  program  auto- 
matically printed  out  every  third  spatial  point  in  both  coordinate 
directions  and  the  points  along  the  boundaries  and  lines  of  symmetry. 

The  user  now  has  complete  control  over  the  printout,  making  for  more 
efficient  use  of  output.  For  maximum  CRIT,  however,  the  program  con- 
tinues to  automatically  check  every  spatial  location. 

A formulation  similar  to  that  employed  to  treat  the  discrete 
stiffeners  was  used  to  replace  the  "equivalent  layer"  treatment  of 
honeycomb  metal  panels  undergoing  catastrophic  damage  (KTYPE=3,  KDAM=1 

or  101) . This  method  involves  two  integration  points  through  the  ■* 

thickness  (LBAR=2),  one  in  each  face  sheet. 
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The  deck  structure  of  DEPROP  was  modified  somewhat  with  the  addition 
of  three  routines  associated  with  the  stiffeners,  MATXIN,  SIGMAB,  and 
STIFF.  Separately,  the  routine  DERV2  was  broken  down  into  two  routines, 
DERV1  and  DERV2.  Several  common  blocks  were  also  changed  in  DEPROP. 

Subsection  5.5  of  this  report  deals  with  a special  version  of 
N0VA-2S  called  N0VA-2LTS.  The  blast  and  aerodynamic  subroutines  have 
been  replaced  by  analytical  and  tape-supplied  pressure  data  to  permit 
correlation  with  experimentation. 

The  final  subsection  documents  an  example  problem  intended  to 
provide  the  user  with  both  an  example  of  program  input  and  modelling, 
and  a check  on  the  computer  program. 

5.1  SUBROUTINES  AND  COMMON  BLOCKS 

Three  new  routines  were  added  to  DEPROP:  MATXIN,  SIGMAB  and  STIFF. 
Subroutine  STIFF  sets  up  all  the  constants  associated  with  stiffened 
panels.  If  the  inertia  coupling  is  included,  it  calls  MATXIN  which 
inverts  the  mass  matrix.  SIGMAB  calculates  the  inelastic  stresses  for 
option  NDERV=2  associated  with  the  stiffeners. 

Subroutine  DERV2  was  separated  into  two  routines,  DERV1  and  DERV2, 
because  of  the  length  of  the  original  routine  and  the  logical  differences 
which  exist.  DERV1  calculates  displacements,  strains,  and  stresses; 

DERV2  calculates  accelerations. 

Table  11  lists  the  107  routines  of  NOVA- 2 S and  table  12  lists  all 
the  associated  common  blocks.  Common  block  IFIRST  was  added  so  that  the 
first  storage  location  (101)  contains  an  integer  variable  monotonically 
increasing  in  value  as  long  as  the  program  is  running  normally.  This 
can  be  checked  by  the  operator. 

Two  versions  of  N0VA-2S,  representing  different  dimensions  for  the 
program  DEPROP,  are  documented.  The  smaller  version  can  be  run  on  the 
Control  Data  Corporation  (CDC)  6600;  the  larger  version  can  only  be  run 
using  LCM  on  the  CDC  176,  or  an  equivalent  system. 
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TABLE  11 

LIST  OF  SUBPROGRAMS  OF  N0VA-2S 


NOVA 

DEPROP 

DEPROB 

NOVA 

WFDZR 

DEPROP 

DEPROB 

BLOCK 

WELL 

BOLT 

COMP1 

IODUM 

WFPRMT 

DERV1 

COMP  2 

SEC 

WFVZR 

DERV2 

COMSET 

NIN 

WFVRMT 

DSET1 

CYCLE 

NEWSL 

AIR 

DSET2 

DAB 

NOVSUM 

WFPKOP 

DSET3 

DEFORM 

RITC 

REFRA 

DTSTEP 

DP  UR 

RITER 

OPT1 

HIM 

EQUILP 

C SETUP 

OPT2 

LEGEND 

EQUILX 

INTP 

OPT3 

LIST1 

FB 

PINIT 

ADVANC 

LIST2 

FBCTL 

SOLVE 

BISH 

MATXIN 

FBSET 

BLAST 

READ 

RELAX? 

FINAL 

XBLAST 

POSTAP 

SIGMA 

FSOL 

HYDRA 

SKIP 

SIGMAB 

PRINT1 

IOPT1 

FPRES 

STIFF 

READ1 

IOPT2 

INTSLO 

RESD 

IOPT3 

PFUSE 

RESET 

ATMOS 

PJUMP 

RLAXB 

MATM62 

POSTW1 

RLAXF 

SHOCK 

POSTW2 

SLAY 

TP1NT 

POSTW3 

STRESS 

INTI 

P0STW4 

STRESX 

INT2 

POSTW5 

STRN1 

WFZR 

POSTW6 

STRN2 

WFPKOD 

POSTW7 

STSET 

WFPR 

PRESS 

TSTEP 

WFPKV 

WFDRMT 

PREW 

SETW 

WPRES 

VCS 

TABLE  12 

COMMON  BLOCKS  AND  SUBPROGRAMS  USING  THEM  IN  NOVA- 2 S 


Common 

Block 

Length 

(Decimal) 

Subprograms 

CDC  6600 

CDC  176 

FIRST 

1 

1 

NOVA*,  DEPROP,  DEPROB,  DEFORM 

CNOVA 

546 

605 

NOVA,  NIN,  NEWSL,  NOVSUM,  RITC, 

CSETUP , BLAST,  XBLAST , PINIT,  FPRES , 
PRESS,  PREW,  WPRES , DEPROP,  DSET1, 

DSET2,  DSET3,  DERV1,  DERV2 , DTSTEP , 

LIST1 , LIST2,  SIGMA,  STIFF,  DEPROB, 
COMP1,  COMP 2,  COMSET,  CYCLE,  DEFORM, 
EQUILP,  EQUILX,  FB,  FINAL,  PRINT1, 

READ1 , STRESS,  TSTEP 

DNOVA 

2858 

2858 

NOVA,  BLOCK,  NIN,  NEWSL,  NOVSUM, 

RITC,  BLAST,  XBLAST,  PINIT,  FPRES, 
POSTW1,  POSTW2,  POSTW3,  POSTW4, 

POSTW5,  POSTW6,  POSTW7,  PRESS,  PREW, 
SETW,  WPRES 

CTLX 

2 

2 

NOVA,  BLAST,  REFRA,  FPRES,  WPRES 

CONSTC 

15 

15 

HYDRA,  IOPT1 , IOPT2,  IOPT3 

SCALEC 

5 

5 

HYDRA,  IOPT1 , IOPT2,  IOPT3,  SHOCK 

WFRT 

13 

13 

SHOCK,  WFPKOD,  WFPR,  WFPKV,  WFDRMT, 

WELL,  WFPRMT,  WFVRMT 

REFRAC 

7495 

7495 

REFRA,  OPT1 , OPT 2,  ADVANC,  READ, 

SKIP 

PW1 

23 

23 

POSTW1 , POSTW2,  POSTW3,  POSTW4, 

POSTW5,  POSTW6,  POSTW7 , SETW, 

WPRES 

CBLK1 

894 

1159 

DEPROP,  BOLT,  DSET1,  DSET2,  DSET3, 

DERV1 , DERV2,  DTSTEP,  LEGEND,  LIST1, 
LIST2,  SIGMA,  SIGMAB,  STIFF 

CBLK2 

4547 

5487 

DEPROP,  BOLT,  DSET1 , DSET2,  DSET3, 

DERV1 , DERV2,  DTSTEP,  STIFF 

*Underlined  routine  in  each  group  owns  that  common  block  in 
segmentation  setup. 
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TABLE  12  (Continued) 


1 

Common 

Block 

Length 

(Decimal) 

Subprograms 

6600 

176 

CBLK3 

12 

12 

DEPROP,  DSET1,  DSET2,  DSET3,  DERV2, 
LEGEND,  SIGMA 

CBLK4 

589 

589 

DEPROP,  DSET1 , DSET2,  DSET3,  DERV1 , 
DERV2,  SIGMA 

CBLK5 

1185 

1185 

DEPROP,  DSET1 , DSET2,  DSET3,  STIFF 

CBLK6* 

25270 

29400 

SIGMA 

CBLK7 

23 

23 

DEPROP,  DSET1 , DSET2,  DSET3,  LIST2, 
SIGMA,  SIGMAB,  STIFF 

CBLK8 

148 

148 

DEPROP,  HIM 

CBLK9 

163 

184 

DEPROP,  DSET1 , DSET2,  DSET3,  LIST1, 

LIST2 

CBLK10 

5415 

6300 

DEPROP,  DSET1 , DSET2,  DSET3,  DERV1, 
DERV2,  LIST1 , LIST2 

CBLK11 

12 

12 

DEPROP,  DSET1 , DSET2,  DSET3 , DERVl 

CBLK12 

22638 

22638 

RELAXP 

CBLK13 

9 

9 

DEPROP,  DSET1 , DSET2,  DSET3,  DTSTEP , 
STIFF 

CBLK15 

5755 

10547 

DEPROP,  DERVl,  DERV2,  DSET1 , LIST1, 
LIST2,  SIGMAB,  STIFF 

CBLK16* 

2944 

5376 

SIGMAB 

CBLK17* 

7203 

7203 

DEPROP,  DERV2,  STIFF 

CBLANK* 

1 



14259 

16501 

DEPROP,  DERVl,  DERV2,  DSET1 , DSET2, 
DSET3,  LIST1,  LIST2,  SIGMA,  SIGMAB, 

STIFF 

*Assigned  to  Level  2 storage  on  CDC  176. 


91 


. 

I ! 

f 

'f  TABLE  12  (Concluded) 


Common 

m 

Block 

Subprograms 

6600 

176 

BLK2 

12717 

12717 

DEPROB,  C0MP1 , COMP 2.  COMSET.  CYCLE. 

DAB,  DEFORM,  DPUR,  EQUILP,  EQUILX, 

FB,  FBCTL,  FBSET,  FINAL,  FSOL,  PRINT1, 
READ1 , RESD,  RESET,  SLAY,  STRESS, 

STRESX,  STRN1 , STRN2,  STSET,  TSTEP , 

VCS 

BLK3 

466 

466 

DAB,  DEFORM,  DPUR,  FSOL,  RESD.  RESET. 
STSET 

BLK4* 

7216 

7216 

RLAXB 

BLK5 

21 

21 

RLAXF 

BLK6 

2369 

2369 

COMP1,  COMP 2,  COMSET,  DEFORM,  FB, 

FBSET,  PRINT1 , STRESS,  STRN1 , STRN2 

*Assigned  to  Level  2 storage  on  CDC  176. 
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5.2  MAXIMUM  PROGRAM  DIMENSIONS 


Nearly  all  of  the  dimensioned  variables  appear  in  labelled  common, 
and  the  current  maximum  dimensions  are  indicated  in  tables  13  through  15. 
The  variable  associated  with  each  dimension  is  listed  in  case  the  user 
should  want  to  change  program  dimensions.  These  tables  should  also  be 
consulted  when  making  up  input  for  the  program  to  be  sure  the  dimensions 
are  not  exceeded. 

There  are  a few  other  changes  to  be  made  when  the  dimensions  are 
changed,  and  these  are  listed  in  table  16.  The  new  integer  variables 
which  represent  maximum  dimensions,  along  with  the  list  of  dimensions  of 
the  new  program  variables,  make  up  table  17. 

5 . 3 PROGRAM  INPUT 

Input  instructions  remain  the  same  (reference  1) , except  for  two 
minor  changes  in  the  NOVA  input,  and  a complete  overhaul  of  DEPROP. 

Groups  10,  27  and  28  of  the  NOVA  data  change  (see  the  new  instructions 
in  tables  18  and  19).  Group  20  permits  the  user  to  make  a response  run, 
yet  still  receive  output  indicating  maximum  response,  or  CRIT.  For 
KDAM  = 100  (or  101)  the  program  executes  as  if  KDAM  = 0 (or  1) , except 
that  only  one  iterative  trial  is  permitted;  otherwise  there  is  no 
difference.  For  KDAM  = 2,  a response  run  without  any  iterative  information 
is  made. 

An  extra  input  parameter,  NU,  is  added  to  Group  27.  This  parameter 
gives  the  user  the  choice  of  either  a uniform  load  or  a circumferentially 
varying  load  on  certain  fuselage  elements.  Table  20  lists  the  new 
options.  Previously,  radomes  and  frame  elements  received  the  varying 
load;  panels,  stringers  and  longerons  received  a uniform  load. 

The  parameter  NFP  locates  the  element  longitudinally  on  the  air- 
craft, while  THETAR  (0R)  locates  the  element  circumferentially.  Both 
parameters  should  correspond  to  the  center  of  the  structural  element,  or 
that  point  on  the  structure  at  which  the  loading  is  desired.  See 
figure  19  of  reference  1 for  the  definition  of  9R  in  the  Aircraft 
Axis  System  (AAS),  and  note  that  the  DEPROB  coordinates  (V,  W)  in  the 
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TABLE  13 

DIMENSIONS  OF  VARIABLES  FOR  NOVA  ROUTINE 


VARIABLE 

DIMENSION 

MMAS1 

28 

NEL 

20 

NFS 

20 

NLE2 

5 

NLEHT 

5 

NLEVT 

5 

NLEW 

5 

NMASS3 

40 

NORMAX 

30 

NORMAX*NEL 

100 

NTE2 

5 

NTEHT 

5 

NTEVT 

5 

NTEW 

5 

NTP1+1 

1000 

HlBAR  in  DEPROP 

^LE  must  be  largest  of  NLEHT,  NLEVT,  NLEW  and  NTE  must  be  largest  of 
NTEHT,  NTEVT , NTEW 

NMASS  must  be  the  largest  of  N in  DEPROB  and  NBAR  in  DEPROP 
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TABLE  14 

DIMENSION  OF  DEPROB  VARIABLES 


if 


VARIABLE 

DIMENSION 

NX1 

10 

N2 

40 

NL3 

8 

NLK 

21 

NSL 

5 

NSSC4 

5 

NSSCT4 

5 

The  program  automatically  assigns  NX  or  fewer  flanges  to  each  layer, 
so  NX  should  usually  be  four  or  six  since  the  sum  of  all  flanges  must 
not  exceed  NLK 

2NMASS  in  NOVA 

3 

For  a uniform  beam,  the  program  may  add  one  layer,  so  the  actual 
limit  on  input  would  ordinarily  be  seven 

4 

The  number  of  distinct  slopes  defined  by  NSSC  and  NSSCT,  excluding 
zero  slope  segments,  must  not  exceed  NSL 
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TABLE  15 

DIMENSIONS  OF  DEPROP  VARIABLES 


VARIABLE 

DIMENSION 

6600  - 

176 

LBAR 

6 

6 

MB 

13 

13 

MBAR1 

23 

28 

MBAR*NBAR*LBAR2 

1805 

2100 

MG 

13 

13 

MGMB3 

49 

49 

NBAR4 

23 

28 

NGNBT»MBAR*NBAR 

361 

420 

NKP 

46 

46 

NL 

8 

8 

NSG 

4 

9 

NSB 

4 

6 

NSMAX 

8 

8 

NSG*MBAR4-NSB*NBAR 

92 

168 

HlMASS  in  NOVA 

2 

This  constraint  is  only  significant  for  an  elastic-plastic  run 
(NDERV*2) . Four  possible  combinations  using  maximum  dimensions 
on  the  6600  are:  (17x17x6),  (19x19x5),  (23x15x5),  (15x23x5). 
Possible  combinations  on  the  176  include:  (28x15x5)  and  (18x18x6). 

3 

The  total  number  of  modes  selected  (MGMB)  from  the  total  possible 
(MB*MG)  cannot  exceed  49 

4NMASS  in  NOVA 
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TABLE  16 

PROGRAM  CHANGES  REQUIRED  BY  DIMENSION  CHANGES 


When  Changing 
the  Dimensions 
Corresponding  to: 

Also  Change  the  Fixed-Point  Number  in  the 
Indicated  Statement 

Routine 

Subroutine 

Location^ 

NORMAX 

NOVA 

NOVA 

5"1 

NORMAX*NEL 

NOVA 

NOVA 

570+6 

NTP1+1 

NOVA 

WPRES 

43"2 

LBAR 

DEPROP 

LEGEND 

300-11 

MG*MB*3 

DEPROP 

RELAXP 

20+1 

NL 

DEPROB 

C0MP1 

Ln 

O 

1 

^4 

NSL*NLK*2* (N+l) 

DEPROB 

COMP1 

50“6 

NLK 

DEPROB 

COMP1 

50-5 

NSL 

DEPROB 

COMP1 

50-4 

N*2+2 

DEPROB 

RLAXB 

40+3 

L +TI 

The  location  code  is  read  as  follows:  s refers  to  the  nth  line 

after  statement  number  s. 
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TABLE  17 

NEW  DEPROP  VARIABLES 


AASB(NSEGB,NBAR,NSB) 

AASG(NSEGG.NSG) 

ALX(LXHAX) 

ASB (NBAR, NSB) 
ASG(NSG) 

AU (MGMB, MGMB) 

AV (MGMB , MGMB ) 
AW(MGMB.MGMB) 

*BETC (NBAR, NSB) 
BEX(LXMAX) 

*BIGJB (NBAR, NSB) 
*BIGJG(NSG) 

*BSTB (NSEGB , NBAR , NSB) 
*BSTG(NSEGG,NSG) 

CA1 

CA2 (NSG) 

CA3 (NSB) 

CllG(NSG) 

C22B (NBAR, NSB) 

DUG  (NSG) 


2 

Area  of  segment  in  B-stiffener,  in  . 

2 

Area  of  segment  in  y-stiffener,  in  . 
Storage  for  stiffener  stress-strain,  a. 

2 

Area  of  cross-section  in  B-stiffener,  in  . 

2 

Area  of  cross-section  in  y-stiffener,  in  . 

Inverse  of  inertia  matrix  in  u-direction, 
[-M]"1,  in4/lb-s2. 


Inverse  of  V-inertia  matrix,  in4/lb-s^. 

4 2 

Inverse  of  W-inertia  matrix,  in  /lb-s  . 


B-position  for  B-stiffener  input,  in  or  deg. 

Storage  for  stiffener  stress-strain,  B. 

4 

Torsion  constant  for  B-stiffener,  in  . 

4 

Torsion  constant  for  y-stiffener,  in  . 


Width  of  segment  in  B-stiffener,  in. 


Width  of  segment  in  y-stiffener,  in. 

2 

Constant  equal  to  2L  R. 

2 — 

Constant  equal  to  6L  (N-l)/a9  hl^  at 
the  fcth  B-position  of  a y-sti?fener. 


2 — 

Constant  equal  to  6L  (M-l)/£hH  at  the 
jth  y-position  of  a B-stiffener. 

* 

Stiffness  constant.  Cn/a,  for  a 
Y-stiffener,  lb/in2. 


Stiffness  constant,  C.?/a 
B-stiffener,  lb/in2. 


for  a 


* 3 

Stiffness  constant.  D,,/®  » for  a 
Y-stiffener,  lb/in2. 
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TABLE  17  (Continued) 


*r  ■ 


D22B(NBAR,NSB) 

D33B(NBAR,NSB) 

D33G(NSG) 

EPOB(NSB) 

EPOG(NSG) 

*EPSB (NSB) 

*EPSG(NSG) 

*ESTRB (NSB) 

*ESTRG(NSG) 

*ETSTRB(NSB) 

*ETSTRG(NSG) 

EXl(LXMAX) 

FllG(NSG) 

F22B(NBAR,NSB) 

*GBARB(NSB) 

*GBARG(NSG) 

*HOB(NBAR,NSB) 

*HOG(NSG) 

*HSTB (NSEGB , NBAR , NSB) 


* 3 

Stiffness  constant,  D?_/a  * 
B-stiffener,  lb/in2. 


for  a 


* , 3 

Stiffness  constant,  D.Ja  , for  a 
B-stiffener,  lb/in2. 


* 3 

Stiffness  constant,  D_,/a  for  a 
y-stiffener,  lb/in2. 


Yield  strain  for  B-stiffener,  in/ in. 


Yield  strain  for  y-stiffener,  in/in. 

Ultimate  tensile  strain  for  B-stiffener, 
in/ in. 

Ultimate  tensile  strain  for  y-stiffener, 
in/ in. 

Elastic  modulus,  E,  for  B-stiffener, 
lb/ in^. 


Elastic  modulus,  E,  for  y-stiffener, 
lb/in2. 


Strain- hardening  slope,  E , for  B-stiffener, 
lb/in2.  C 


Strain- hardening  slope,  E. , for  y-stiffener, 
lb/ in2.  C 

Storage  for  stiffener  stress-strain,  e^. 

* , 2 

Stiffeners  constant,  F../a  , for  a 
Y-stiffener,  lb/in2. 

* 2 

Stiffness  constant,  F_2/a  » ^or  a 
B-stiffener,  lb/in*. 

— 2 

Shear  modulus,  G,  for  a 6-stiffener,  lb/in  . 

_ 2 

Shear  modulus,  G,  for  a y-stiffener,  lb/ in  . 

Gap  between  8-stiffener  and  panel,  hQ,  in. 

Gap  between  y-stiffener  and  panel,  hQ,  in. 

Distance  from  inner  panel  surface  to  the 
fcth  segment  of  B-stiffener,  hQ,  in. 
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TA3LE  17  (Continued) 


*HSTG (NSEGG , NSG) 

*KCOUP 

*KPB(NKP) 

*KPG(NKP) 

*KSB(NSB) 

KSBX (MBAR) 
*KSG(NSG) 
KSGX(NBAR) 

KSTIF 

KSUMB(NSGMB) 

*KSUPB(NSB) 


Distance  from  inner  panel  surface  to 
the  fcth  segment  of  y-stiffener,  hQ,  in. 

Code  indicating  whether  the  full  inertia 
coupling  is  to  be  included:  0,  only 
diagonal  terms;  1,  yes. 

Mesh-point  number  ((3),  when  paired  with 
KPG,  specifies  printout  locations. 

Mesh-point  number  (y) , when  paired  with 
KPB,  specifies  printout  locations. 

Gamma-point  location  of  B-stiffener 
(y  grid-point  number) . 

1 

Beta-stiffener  number  corresponding  to 
each  y grid  point  (zero  for  no  stiffener) . 

Beta-point  location  of  y-stiffener  ; 

(B  grid-point  number).  1 

Gamma-stiffener  number  corresponding  to 
each  B grid  point  (zero  for  no  stiffeners). 

Total  number  of  stiffeners  in  model. 

Number  of  z points  in  stiffeners  which 
have  not  yielded  during  response. 

Support  code  for  outstanding  leg  of 
B-stiffener  (0,  1,  or  2). 


*KSUPG(NSG) 

KYX(LXMAX) 


LXMAX 


Support  code  for  outstanding  leg  of 
y-stiffener  (0,  1,  or  2). 

Code  in  elastic-plastic  response  indicating 
number  of  times  a stiffener  integration 
point  has  yielded,  unloaded,  etc. 


in 

NSEGB (J) 

|» 

100 

; 

! 


Total  number  of  integration  points 
stiffeners,  equal  to 

NSG  NSB 

MBAR  • Z NSEGG(I)  + NBAR  • Z 

1-1  J-l 


TABLE  17  (Continued) 


MFIRST 

NFIRST 


Code  indicating  whether  any  stiffener  has 
yielded  (0,  no;  1,  yes). 

Code  indicating  first  pass  through 
routine  SIGMAB. 


*NKP 


*NSB 


Number  of  spatial  grid  points  for  which 
printout  of  strains,  stresses,  displacements, 
and  pressures  is  required. 

Number  of  B-stiffeners  (stiffener  parallel 
to  the  B-axis). 


*NSEGB(NSB) 


Number  of  segments  (layers)  in  the 
B-stiffener . 


*NSEGG(NSG) 


Number  of  segments  (layers)  in  the 
Y-stiffener . 


*NSG 


Number  of  y-stiffeners  (stiffeners  parallel 
to  the  y-axis) . 


NSGMB 


Total  number  of  grid  points  involving 
stiffeners,  equal  to  NSG*MBAR+NSB*NBAR. 


NSMAX 

*NSTB (NSB) 

*NSYMB 

*NSYMG 

NUSE (NBAR.MBAR) 

PRLU(MGMB) 

PRLV(MGMB) 


Maximum  number  of  segments  in  any  stiffener  - 
gamma  or  beta. 

Number  of  B-stiffeners  which  define  cross- 
section  for  B-stiffener  (<NBAR) . 

Symmetry  code  for  panel  model  in  B-direction: 
0-symmetric;  1-not  symmetric. 

Symmetry  code  for  panel  model  in  y -direction: 
0-symmetric;  1-not  symmetric. 

Use-code  for  the  spatial  integration 
stations:  0-not  used;  1-printout  only; 
2-integration  only;  3-both. 

Diagonal  terms  of  u-stiffness  matrix, 
inVlb-s2. 

Diagonal  terms  of  v-stiffness  matrix, 
in4/lb-s2. 
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TABLE  17  (Concluded) 


PRLW(MGMB) 

*RHOSTB(NSB) 

*RHOSTG(NSG) 

*SIDEB(NSB) 


*SIDEG(NSG) 

*SIGOBC (NSB) 

*SIGOBT(NSB) 

*SIGOGC (NSG) 

*SIGOGT(NSG) 

SIXl(LXMAX) 

SX(LXMAX) 

ZFB (NSMAX , NSG+NSB*NBAR) 
ZSTB(2 ,NBAR,NSB) 

ZSTG(2,NSG) 


Diagonal  terms  of  w-stiffness  matrix, 
inVlb-s2. 

2 

Density  of  6-stiffeners,  pg,  lb-s  /in 

2 

Density  of  y-stif feners,  p , lb-s  /in 


4 

4 


Input  code  designating  location  of  6-stiffener: 
-1.0,  inner  (exterior  to  panel);  +1.0, 
outer;  +2.0,  internal  (honeycomb  only); 

+3.0,  inner  with  panel  crimped  at  stiffener 
locations.  After  input,  variable  takes  on 
a value  of  1.0  unless  crimped,  when  it  is  0.0. 


Same  as  SIDEB,  only  for  y-stiff eners . 

Compressive  yield  stress  for  6-stiffener, 
lb/ in2. 

Tensile  yield  stress  for  6-stiffener, 
lb/in2. 

Compressive  yield  stress  for  y-stif f eners , 
lb/ in2. 

Tensile  yield  stress  for  y-stiffeners, 
lb/ in2. 

Storage  for  stiffener  stress-strain,  o^. 

2 

Stress  in  stiffeners,  lb/in  . 

z-position  in  stiffener  for  integration,  in. 

z-position  on  inner  and  outer  surfaces  of 
6-stiffener,  in. 

z-position  on  inner  and  outer  surfaces  of 
y-stif fener,  in. 


Asterik  (*)  indicates  an  input  variable. 
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TABLE  18 

REVISION  OF  NOVA  DATA  GROUP  10 


Group  10:  (2112)  KDAM,  KALT 

Range  iteration/ damage  code  (KDAM) 

0,  iterate  to  determine  range  at  which 
permanent  damage  first  occurs. 

1,  iterate  to  determine  range  at  which 
catastrophic  damage  occurs. 

2,  determine  structural  response  only  at 
specified  range. 

100,  same  as  KDAM  * 0,  except  only  1 trial 
in  iteration. 

101,  same  as  KDAM  * 1,  except  only  1 trial 
in  iteration. 

Constant  altitude  (KALT) 

0,  no  restriction  on  iteration. 

1,  iteration  restricted  to  constant  altitude. 

Note:  KALT  is  not  necessary  for  KDAM  » 2.  Otherwise 
KALT  must  be  1 if  both  KB  and  KGRD  are  1. 


If  KDAM  - 3.  skip  to  GROUP  12. 


Group  27 : 


(2112)  NFP , NU 


The  n_  ber  of  the  fuselage  section  (from  the  table 
of  values  supplied  in  Group  24)  at  which  pressures 
are  desired,  i.e.,  the  section  at  which  the  structural 
element  if  located.  (NFP) 

Code  for  circumferential  variation  of  load: 

0 - circumferentially  varying  load. 

1 - uniform  load. 


Note:  If  the  structural  element  is  a radome  (KTYPE>7), 


to  GROUP  29. 


Group  28:  (F12.1)  THETAR 

Angular  location  of  the  center  of  the  structural 
element  on  the  circumference  of  the  equivalent,  circular 
section  for  the  fuselage  (figure  19  of  reference  1).  For 
a uniform  load  (NU  = 1)  this  locates  the  point  at  which 
the  pressures  are  applied.  (-ir  <_  0R  <_  it)  (THETAR),  rad. 


TABLE  20 

GROUPS  27  AND  28  DATA  OPTIONS  FOR  LOADING  FUSELAGE  ELEMENTS 


Fuselage 

Element 

KTYPS 

Group  28  data 
for  Uniform  Load 

(NU  « 1; 

Group  28  data  for 
non-uniform  Load 

(NU  = 0) 

Panel 

< 5 

en 

e„ 

R 

R 

Stringer 

6 

9r 

Not  possible 

Frame 

7 

6r 

0 can  be  anything* 

Radome 

8,9 

Not  possible 

0^  not  inputted* 

*For  these  cases,  the  V,  W 

coordinates  (or  0^,  0 

2)  (LAS)  in  DEPROB 

must  locate  the  beam  elements  circumferentially. 


! 
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Local  Aircraft  System  (LAS)  align  with  (y,  z)  in  the  AAS.  The  DEPROB 

angle  6,  though,  is  not  defined  in  the  same  manner  as  6 (page  274  of 

R 

reference  1) . 

The  only  input  instruction  which  is  modified  in  the  DEPROB  section 
deals  with  the  KSUP  parameters  of  Group  1.  KSUP  is  not  needed  for  cases 
where  there  is  not  a threshold-of-permanent  change  requirement;  i.e., 
for  KDAM  - 1,  2 or  101. 

The  DEPROP  input  has  been  changed  significantly,  although  much  of 
it  remains  the  same.  Even  so,  the  entire  set  of  input  instructions  is 
documented  in  table  21  to  facilitate  the  preparation  of  an  input  deck. 
Specific  input  instructions  follow  several  paragraphs  of  general  remarks. 
The  user  is  reminded  to  compare  all  input  variables  with  the  maximum 
dimension  provided  in  the  program,  as  delineated  in  table  15.  This  is 
very  important  since  the  program  does  not  attempt  to  check  the  input  for 
such  violations. 

Group  1 contains  the  number  of  modes  to  be  used  in  the  solution  and 
the  number  of  integration  points  to  be  used.  The  accuracy  of  the 
solution  is  based  on  the  degree  of  convergence  of  stress  and  strain 
quantities.  These  quantities  converge  less  rapidly  than  the  radial 
displacement.  Also,  cases  involving  a clamped  edge  condition  will 
converge  less  rapidly  than  simply-supported  cases.  Since  both  computer 
time  and  accuracy  increase  with  more  modes  and  points,  a trade-off 
usually  becomes  necessary.  Although  the  program  allows  up  to  13  gamma 
modes  and  13  beta  modes  to  be  used,  only  a small  number  of  modal  com- 
binations are  normally  required,  as  will  be  discussed  shortly. 

The  actual  mode  numbers  are  specified  in  Groups  2 and  3.  The 
maximum  value  that  the  mode  numbers  can  assume  in  the  program  is  19. 

When  symmetry  is  taken  in  either  direction  (Group  4,  or  if  the  pressure 
loading  is  symmetrically  oriented,  only  the  odd  numbered  modes  (1,  3, 

5,  ...)  are  required  in  that  direction. 
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In  general,  a minimum  of  sixteen  modal  combinations  should  be  used 
for  a symmetric  panel,  and  it  is  recommended  that  at  least  25  be  used 
for  clamped  panels  where  edge  stresses  and  strains  are  important.  The 
maximum  number  of  combinations  permitted  is  49  (see  the  discussion  of 
data  groups  6 and  7) . 

Spatially,  the  optimum  number  of  integration  points  (MBAR  and  NBAR) 
for  a full  panel  should  be  approximately  two  times  the  maximum  mode 
number  used  in  that  direction,  plus  three.  However,  when  NBN  or  MGM  is 
large,  this  condition  may  not  be  satisfied  for  nonsymmetrical  panels, 
since  MBAR  and  NBAR  are  dimensioned  at  23  (28  on  the  CDC  176)  in  the 
program  (see  table  15).  For  symmetric  solutions,  MBAR  (or  NBAR)  need 
only  be  approximately  one-half  the  value  for  a full  panel  since  only 
one-half  (or  one  quarter)  of  the  panel  is  actually  analyzed  in  the 
solution.  For  a nonsymmetric  condition,  MBAR  (or  NBAR)  must  be  an  odd 
number.  For  an  elastic-plastic  solution,  a minimum  of  four  integration 
points  through  the  thickness  is  recommended,  and  a maximum  of  six  is 
provided  in  the  program.  The  exception  is  a metal  honeycomb  panel  where 
only  two  points  are  used. 

In  Group  5,  the  user  is  given  the  option  of  a purely  elastic 
solution,  or  an  elastic-plastic  solution.  The  elastic-plastic  option 
will  tend  to  be  slower  and  require  more  computer  memory.  The  second 
option  (elastic-plastic)  must  be  used  for  metal  panel  solutions  which 
iterate  to  a catastrophic  damage  level  (KTYPE  = 1,  3;  KDAM  ■ 1,  101). 

Group  5 also  specifies  the  number  of  stiffeners  in  the  model  and 
whether  the  full  coupled  mass  matrix  is  used,  or  only  the  diagonal 
terms.  The  advantage  of  only  using  the  diagonal  terms  is  to  reduce  the 
computer  time;  however,  the  savings  is  not  that  much  and  it  is  recom- 
mended that  the  full  matrix  be  used  for  a more  accurate  solution 
(KCOUP  - 1). 

A gamma  stiffener  is  defined  as  a stiffener  running  parallel  to  the 
(x,  y)  axis  and  the  beta  stiffener  is  defined  similarly.  These  stiffeners 
must  be  located  on  either  a y or  (3  grid  line  and  are  assumed  to  be 
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attached  at  each  grid  line  intersection.  If  all  the  gamma  stiffeners 
are  of  identical  construction,  NSG  should  be  inputted  as  a negative 
number  (i.e.,  -3  for  three  identical  stiffeners)  which  will  eliminate 
unnecessary  input.  And  the  same  instruction  applies  to  NSB  for  beta 
stiffeners. 


Groups  6 and  7 provide  a mechanism  for  selecting  a maximum  of  49 
modal  combinations  from  a 13  by  13  combination  array  (MG»MB=13) . Thus, 
the  more  significant  modal  combinations  for  an  optimal  solution  with 
respect  to  accuracy  and  computer  time  can  be  selected  and  the  other 
combinations  eliminated.  A general  rule  of  thumb  is  to  eliminate  the 
higher  frequency  modes  which  are  usually  associated  with  modal  com- 
binations having  the  larger  MG+MB  values.  An  example  of  this  would  be 
the  selection  of  MG=*MB»7  for  a symmetric  problem,  but  eliminating 
24  combinations  as  indicated  in  figure  41.  The  relative  importance  of 
each  modal  combination  can  be  evaluated  by  examining  the  response  output 
and  comparing  the  magnitudes  of  the  displacement  coefficients. 

Groups  8 and  9 are  responsible  for  selecting  the  points  in  the 
integration  grid  for  which  printout  of  strains,  stresses,  displacements, 
and  pressures  is  required.  Strains  and  stresses  are  computed  at  the 
inner  and  outer  surfaces  of  the  panel  layers  and  stiffeners.  Each  point 
in  the  grid  is  designated  by  a pair  of  integers,  the  first  integer 
referring  to  the  gamma- posit ion,  the  second  to  the  beta-position. 

Actual  positions  are  found  from 


x _ L (i-D 
2 — 
(M-l) 


I - 1 M 

(symmetric  in  x-direction) 


x - l 

M-l 


I » 1, . . . ,M 

(full  in  x-direction) 


and  similar  expressions  for  y (or  9).  Vor  example,  the  corner  point 
in  a symmetric  panel  would  be  numbered  v.1,1);  the  center  (MBAR.NBAR)  . 
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It  should  be  noted  that  the  maximum  response  values  associated  with 
determining  CRIT  for  an  iterative  run  are  calculated  at  every  grid 
point,  independent  of  the  printout  selected  in  groups  8 and  9. 

The  length  and  width  (XLP  and  THETAO)  selected  in  Group  11  represent 
the  total  dimensions  of  the  panel,  even  if  only  1/2  or  1/4  is  analyzed 
in  a symmetric  case. 

1 Group  15  provides  the  data  required  for  computing  allowable 

stresses  for  honeycomb  panels.  The  core  cell  size  (DC)  is  defined  as 
the  distance  between  opposite  flat  sides  of  the  honeycomb  cell. 

Groups  18-23  describe  the  gamma  stiffeners  and  groups  24-29 
describe  the  beta  stiffeners.  The  code  SIDEG  (or  SIDEB)  indicates 
whether  a stiffener  is  attached  on  the  inside,  outside,  or  within  the 
panel,  as  illustrated  in  figure  42  (a-d) . 

The  stiffener  cross  section  is  modelled  in  a manner  similar  to  that 
used  in  DEPROB  for  beam  elements.  The  section  is  broken  down  into 
rectangular  layers  called  segments.  These  segments  also  serve  the 
purpose  of  flanges  for  integration  through  the  thickness,  so  the  analyst 
must  assign  enough  segments  to  be  able  to  adequately  model  the  cross 
section  for  an  inelastic  problem.  A maximum  of  eight  segments  is 
permitted  (NSEGG,  NSEGB) . 

Although  stiffeners  can  be  of  different  material,  any  one  stiffener 
is  assumed  to  be  of  homogeneous  construction;  i.e.,  each  segment  is 
composed  of  the  same  material.  It  is  further  assumed  that  if  the  panel 
is  made  of  plastic  material,  the  stiffeners  are  plastic,  and  similarly 
if  the  panel  is  metal  (or  metal  face  sheets  on  honeycomb) , the  stiffeners 
are  metal. 

Gamma  stiffeners  also  must  be  uniform  in  the  spanwise  direction, 

i 

but  beta  stiffeners  can  have  variable  cross  section.  This  is  accom- 
plished by  specifying  the  cross-sectional  shape  at  one  or  more  arbitrary 
beta  locations  (BETC) . The  program  linearly  interpolates  between  points, 
if  necessary,  to  provide  data  at  every  spatial  grid  point.  Obviously, 


if  only  one  cross  section  is  specified,  the  stiffener  is  assumed  to  be 
uniform  and  it  makes  no  difference  what  beta  location  (BETC)  is  used. 
Constant  cross  section  is  assumed  for  beta  points  outside  the  domain  of 
BETC. 

Figure  42  Illustrates  the  four  types  of  stiffeners.  Case  (a) 
represents  an  "outer"  stiffener  with  a spacer  of  thickness  hQ  (HOG, HOB) 
located  between  the  panel  and  the  stiffener.  This  space  might  represent 
insulating  material  or  a gap  created  by  another  stiffener  running  orthogonal 
to  the  one  being  modelled.  For  simplicity  the  stiffener  is  modelled  as 
an  "I"  section  with  three  segments;  in  an  elastic-plastic  model  the  web 
(segment  2)  should  probably  be  broken  down  into  three  or  four  segments. 

The  input  parameters  HSTG(B)  correspond  to  the  h^,  hj  and  h^  shown  and 
are  measured  relative  to  the  inner  surface  of  the  panel. 

Case  (b)  is  nearly  identical,  except  that  the  stiffener  is  called 
an  "inner"  stiffener  as  it  is  located  on  the  inside  of  the  panel.  The 
h^  are  again  referenced  relative  to  the  inner  surface  of  the  panel,  and 
are  also  inputted  as  positive  numbers,  as  is  hQ. 

Cases  (c)  and  (d)  represent  a panel  of  sandwich  (or  honeycomb) 
construction.  In  the  first  case  the  stiffener  lies  within  the  panel  and 
the  input  parameters  are  defined  as  before.  The  second  case,  however, 
considers  the  case  when  the  panel  is  crimped  in  order  to  attach  the 
stiffener,  and  thus  the  stiffener  may  or  may  not  lie  totally  within  the 
panel  section.  If  not,  the  h^'s  must  be  defined  in  a different  manner, 
as  shown  in  figure  42  (d) . To  begin  with,  the  segments  must  be  defined 
so  that  there  is  a division  between  two  segments  at  the  imaginary  inner 
panel  surface.  The  segments  are  ordered  from  that  point  outward  as  far 
as  possible  until  a switch  over  is  required,  as  shown.  In  order  to  flag 
this  switch  over,  which  occurs  at  the  Jl  - segment,  the  parameter  h^ 

is  made  negative.  In  this  case,  h^  would  be  inputted  as  a negative 
number  - all  the  others  are  positive. 

Group  30  contains  the  modal  components,  5 , for  the  Initial 

mn 

radial  imperfections.  The  analyst  must  compute  the  6 's  from  measured 

mn 
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TABLE  21 


Group  1 


Group  2 


Group  3 


Group  4 


Group  5 


DEPROP  INPUT 


(5112)  MG,  MB,  MBAR,  NBAR,  LBAR 

Number  of  gamma  modes  to  be  used.  (MG) 

Number  of  beta  modes  to  be  used.  (MB) 

Number  of  gamma  integration  points  actually  used  over 
the  portion  of  the  panel  analyzed.  Must  be  an  odd 
number  for  full  panel  (see  Group  4) . (MBAR) 

Number  of  beta  integration  points  actually  used  over 
the  portion  of  the  panel  analyzed.  Must  be  an  odd 
number  for  full  panel  (see  Group  4) . (NBAR) 

Number  of  z integration  points  used  through  the  thick- 
ness. (LBAR)  Should  be  2 for  KTYPE-3. 

[Not  needed  for  NDERV*1  (see  Group  5)] 


(6112)  (MGM(I) , 1-1,  MG) 
Gamma  mode  numbers,  m. 


(6112)  (NBN(I) , 1-1,  MB) 
Beta  mode  numbers,  n. 


(2112)  NSYMG,  NSYMB 

Symmetry  code  in  gamma  direction  (NSYMG) : 

0,  symmetry  assumed  (0  ± Y ± ir / 2 ) 

1,  no  symmetry  (0  £ y _<  ir) 

Symmetry  code  in  beta  direction  (NSYMB): 

0,  symmetry  assumed  (0  < 8 < ir/2) 

1,  no  symmetry  (0  ± 8 n) 


: (6112)  NPLT,  NBND,  NDERV,  NSG,  NSB,  KCOUP 

Panel  type  (NPLT) : 

0,  flat  panel 

1,  cylindrical  panel 

Boundary  condition  code  (NBND).’ 


' 
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DEPROP  INPUT  (Continued) 


Note: 


Note: 


Note: 


Note: 


y-direction 

1,  clamped-clamped 

2,  simple-simple; 

3,  clamped-clamped; 

4,  simple-simple 

5,  clamped-simple; 

6,  clamped-clamped; 

7,  clamped-simple; 

8,  simple-simple; 

9,  clamped-simple; 


8-direction 

clamped-clamped 

simple-simple 

simple-simple 

clamped-clamped 

clamped-clamped 

clamped-simple 

simple-simple 

clamped-simple 

clamped-simple 


Whenever  a clamped-simple  condition  is  selected,  the  full 
panel  is  analyzed  in  that  direction,  and  NSYMG,  NSYMB,  MBAR 
and  NBAR  should  reflect  this. 

Response  option  (NDERV) : 

1,  elastic  only 

2,  elastic-plastic 

NDERV  must  be  2 for  KTYPE-1  or  3 whenever  KDAM-1  or  101. 


Number  of  y stringers  (NSG) 

Number  of  6 stringers  (NSB) 

A negative  value  for  NSG  (or  NSB)  means  that  all  y (or  8) 
stringers  are  identical.  Only  one  set  of  input  data  will 
be  required  in  that  case. 

Mass-matrix  coupling  code  (KCOUP) 

0,  no  compiling 

1,  compiling 

KCOUP  will  be  zero  if  NSG*NSB-0. 


Group  6:  (112)  NNOUT 

Number  of  modal  combinations  to  be  el iminated  from 
solution  (NNOUT). 

(0  £ NNOUT  < MG*MB) 

If  NN0UT"0,  skip  to  Group  8. 


Group  7:  (2112)  MOUT(I),  NOUT(I) 

Gamma  mode.  (MOUT(I)) 

Beta  mode.  (NOUT(I)  ) 

Repeat  Group  7 for  I»l,  NNOUT.  The  cards  in  Group  7 may  be  arranged 
in  any  order. 
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DEPROP  INPUT  (Continued) 


Group  8:  (112)  NKP 

Number  of  spatial  points  at  which  printout  of  stresses, 
strains,  displacements,  reactive  forces  and  pressures 
are  requested.  If  NKP-0,  all  of  the  above  information 
will  be  suppressed.  (NKP) 


If  NKP-0, 


to  Grout 


Group  9:  (2112)  KPG(I),  KPB(I) 

Integration  point  in  gamma-direction  at  which  printout 
is  requested.  Points  are  ordered  1-MBAR,  beginning  at 
Y-0,  and  evenly  spaced  from  there.  (KPG(I)) 

Integration  point  in  beta-direction  at  which  printout 
is  requested.  Points  are  ordered  1-NBAR,  beginning  at 
8-0,  and  evenly  spaced  from  there.  (KPB(I)) 


Note:  These  two  indices  are  taken  as  pairs  where  each  pair  desig- 

nates a particular  spatial  point.  The  pairs  may  be  specified 
in  any  order. 


Repeat  Group  9 for  1-1, 


Group  10:  (112)  NL 

Number  of  layers.  (NL) 

(NL  must  be  1 for  KTYPE-1 , and  3 for  KTYPE-3) 


Group  11:  (3F12.1)  XLP,  THETAO , A 

Full  length  of  panel,  i,  in.  (XLP) 

Full  width  of  flat  panel,  b (short  direction),  \ 

in.  (NPLT-0)  f 

or  / (THETAO) 

Full  subtended  angle  of  cylindrical  panel,  0-,  1 

deg.  (NPLT-1)  U ' 

Radius  of  cylindrical  panel,  in.  (A) 

(Not  needed  for  NPLT-0) 


If  NDERV-2. 


to  Grout 
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DEPROP  INPUT  (Continued) 

( 

] 

Group  12:  (2F12.1)  HM(I),  RHOM(I) 

Distance  (h)  from  the  inner  panel  surface  to  the  outer 
surface  of  layer  I,  in.  (HM(I)) 

Mass  density  of  layer  I,  lb-s^/in^.  (RHOM(I)) 

Group  13:  (5F12.1)  EX(I),  ET(I),  XXNU(I),  THNU(I),  GXT(I) 

Modulus  of  elasticity  in  the  x-direction,  psi.  (EX(I)) 

Modulus  of  elasticity  in  the  theta-direction,  psi.  (ET(I)) 

Poisson's  ratio  in  the  x-direction.  (XXNU(I)) 

Poisson's  ratio  in  the  theta-direction.  (THNU(I)) 

Shear  modulus,  psi.  (GXT(I)) 

Group  14:  (2F12.1)  SAT(I),  SAC(l) 

Tensile  yield  stress  of  metal  panels;  tensile  ultimate 
stress  for  plastic  panels,  psi.  (SAT(I)) 

Absolute  value  of  compressive  yield  stress  for  metal 
panels;  absolute  value  of  compressive  ultimate  stress 
for  plastic  panels,  psi.  (SAC(I)) 

If  KTYPE»1,2,  or  5,  skip  to  Group  18. 

If  KDAM-1,2,  or  101,  skip  to  Group  18. 

Group  15:  (3F12.1)  EC,GC,DC 

Core  modulus  of  elasticity  parallel  to  core  depth, 
psi.  (EC) 

Shear  modulus  of  core,  psi.  (GC) 

Core  cell  size,  in.  (DC) 

Skip  to  Group  18. 


DEPROP  INPUT  (Continued) 


Group  16:  (3F12.1)  HM(I) , RHOM(I) , EM(I) 

Distance  (h)  from  inner  shell  surface  in  the  outer 
surface  of  layer  I,  in.  (HM(I)) 

9 / 

Mass  density  of  layer  I,  lb-s  /in  . (RHOM(I)) 

Modulus  of  elasticity,  psi.  (EM(I)) 

Note:  EM(2)  need  not  be  specified  for  a metal  honeycomb  material 

(KTYPE-3) 

Repeat  Group  16  for  1-1,  NL. 


Group  17:  (4F12.1)  TNU,  SIGO,  EP,  EPSIF 

Poisson's  ratio.  (TNU) 

Yield  stress  for  a metal  panel,  psi.  (SIGO) 

Strain  hardening  modulus  (Efc) , psi.  (EP) 

Ultimate  strain,  in/in.  (EPSIF) 

(Not  necessary  for  KDAM“2) 

If  NSG“0,  skip  Groups  18-23,  which  pertain  to  gamma  stiffeners. 


Group  18:  (6112)  KSG(I),  1-1,/NSG/ 

Beta-point  number  corresponding  to  gamma  stiffener  location. 


Group  19:  (6F12.1)  SIDEG(I),  ESTRG(I) , GBARG(I) , RHOSTG(I), 

SIGOGT (I) , SIGOGC(I) 

Code  designating  type  of  stiffener  (SIDEG) : 

-1.0,  Inner 
+1.0,  Outer 

+2.0,  Internal  (honeycomb  panel  construction) 
+3.0,  Inner  (crimped  honeycomb  panel) 

2 

Elastic  modulus^  E,  lb/ in  (ESTRG) 

Shear  modulus.  G,  lb/in2  (GBARG) 

Density,  ps,  lb-s2/in4  (RHOSTG) 

Tensile  yield  stress,  lb/ in2  SIGOGT) 

Compressive  yield-stress,  lb/in2  (SIGOGC) 

If  NDERV-1,  skip  Group  20. 
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DEPROP  INPUT  (Continued) 


Group  20:  (2F12.1)  ETSTRG(I),  EPSG(I) 

— 2 

Strain-hardening  slope,  E , lb/ in  (ETSTRG) 

Ultimate  tensile  strain,  Ce  , in/ in  (EPSG) 

Group  21:  (2112)  NSEGG(I) , KSUPG(I) 

Number  of  segments  (NSEGG) 

Support  code  for  outstanding  leg  (KSUPG) : 

0,  no  outstanding  leg. 

1,  outstanding  leg  supported  at  one  end. 

2,  outstanding  leg  supported  at  both  ends  or  has 
two  corners . 


Group  22:  (2F12.1)  BIGJG(I),  HOG(I) 

4 

Torsion  constant  for  stiffener,  J,  in  (BIGJG) 

Gap  between  stiffener  and  panel,  hQ,  in  (HOG) 

Group  23:  (2F12.1)  HSTG(L,I),  BSTG(L,I) 

Distance  from  inner  panel  surface  to  the  furthest  edge  of 
the  segment,  h^,  in.  (HSTG) 

Width  of  Jtch  segment,  b^,  in.  (BSTG) 

Note  - HSTG  is  always  a positive  number  except  for 
SIDEG-3,  for  the  fl,*L  segment  which  causes  HSTG  to  switch 
directions.  See  figure  42. 

th 

Repeat  Group  23  for  all  segments  in  the  I stiffener. 

Unless  NSG  was  read  in  as  a negative  number,  repeat  Groups  19-23  for 
each  gamma  stiffener. 

If  NSB~0,  skip  Groups  24-29,  which  pertain  to  the  beta  stiffeners. 


¥ 


Group  24:  (6112)  KSB(I),  1-1,/NSB/ 

Gamma-point  number  corresponding  to  beta  stiffener  location. 


Group  25:  (6F12.1)  SIDEB(I),  ESTRB(I),  GBARB(I) , RHOSTB(I) , 

SIGOBT(I) , SIGOBC(I) 
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DEPROP  INPUT  (Continued) 


! 

i 


Code  designating  type  of  stiffener  (SIDEB) : 

-1.0,  Inner 

( +1.0,  Outer 

s +2.0,  Internal  (honeycomb  panel  construction) 

+3.0,  Inner  (crimped  honeycomb  panel) 

7 — , 2 

Elastic  modulus^  E,  lb/in  (ESTRB) 

Shear  modulus,  G,  lb/in2  (GBARB) 

Density,  ps,  lb-sZ/in*  (RHOSTB) 

Tensile  yield  stress,  lb/ in2  (SIGOBT) 

Compressive  yield  stress,  lb/in^  (SIGOBC) 

If  NDERV°1,  skip  Group  26. 


Group  26:  (2F12.1)  ETSTRB (I) , EPSB(I) 

_ 2 

Strain-hardening  slope,  E , lb/in  (ETSTRB) 

Ultimate  tensile  train,  e^,  in/in  (EPSB) 

Group  27:  (3F12.1)  NSEGB(I) , KSUPB(I),  NSTB(I) 

Number  of  segments  (NSEGB) 

Support  code  for  outstanding  leg  (KSUPB) : 

0,  no  outstanding  leg 

1,  outstanding  leg  supported  at  one  end. 

2,  outstanding  leg  supported  at  both  ends  or  has  two 
corners. 

Number  of  8-stations  used  to  define  cross  section.  (NSTB) 

Notes:  NSTB  should  be  1 for  a uniform  cross  section  and  BETC  can  be 

anything . 

NSTB  must  be  < NBAR. 


Group  28:  (3F12.1)  BIGJB(K.I),  H0B(K,I),  BETC(K,I) 

Torsion  constant  for  station,  J,  in^.  (BIGJB) 
Gap  between  stiffener  and  panel,  h , in.  (HOB) 

Beta  position  of  the  station,  ?n  or  deg.  (BETC) 
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DEPROP  INPUT  (Concluded) 


Group  29:  (2F12.1)  HSTB(L,K,I),  BSTB(L,K,I) 

Distance  from  inner  panel  surface  to  the  furthest  edge  of 
the  1th  segment,  h^,  in.  (HSTB) 

Width  of  segment,  bn,  in.  (BSTB) 

Note:  HSTB  is  always  a positive  number  except  for  SIDEB=3.0,  for 

the  H-L  segment,  which  causes  HSTB  to  switch  directions. 

See  figure  42. 

Repeat  Group  29  for  all  segments  at  the  kth  station  of  the  Ith  stiffener. 

Unless  NSB  was  read  in  as  a negative  number,  repeat  Groups  25-29 
for  each  fl-stiffener. 


Group  30:  (6F12.1)  ((FG(N,M) , N-l.MB),  M-l.MG) 

Modal  displacement  coefficients  for  initial  radial 
imperfections,  in.  (FG (N,M)) 


Group  31:  (3F12.1)  DELTIM,  TSTOP,  PRINT 

Integration  time  increment,  sec.  If  DELTIM-0.0,  the 
program  determines  the  time  increment  required  for 
stability.  (DELTIM) 

Integration  stop  time,  sec.  (TSTOP) 

Print  frequency  (integration  steps  per  printout) . If 
PRINT-0.0,  printout  of  intermediate  data  will  be  sup- 
pressed. (PRINT) 
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data  using  the  integration  technique  applied  to  Fourier  series  coeffi- 
cients. Generally,  such  data  will  not  be  available,  and  zero  values 
should  be  specified  for  the  6 's.  The  capability  of  considering 

initial  imperfections  also  enables  the  analyst  to  determine  the  sensi- 
tivity of  panel  response  to  initial  imperfections. 

Group  31  provides  the  integration  time  increment,  the  response  stop 
time,  and  printout  interval.  If  the  user  specifies  a zero  time  incre- 
ment, the  program  computes  an  appropriate  At  which  in  most  cases  will 
give  a stable  solution.  It  should  be  noted,  however,  that  stiffeners 
are  neglected  in  the  computation,  so  it  is  possible  a smaller  At  will  be 
required  for  some  stiffened  panels.  Because  the  At  is  approximate,  the 
analyst  may  want  to  make  comparable  runs  using  different  At's.  In 
general,  an  elastic  solution  which  is  numerically  stable  will  be  accurate. 
Hence,  the  optimum  At  is  the  largest  which  remains  stable.  For  an 
elastic-plastic  solution,  however,  the  accuracy  of  the  solution  may 
deteriorate  slightly  as  the  point  at  which  the  solution  diverges  is 
approached.  Once  a time  increment  is  selected,  it  should  be  valid  for 
other  orientations  and  moderate  changes  in  response  level. 

Although  the  stop  time  can  vary  a great  deal,  the  total  number  of 
integration  steps  required  to  capture  peak  response  will  be  roughly 
between  500  and  1500.  One  exception  to  this  may  be  a curved  panel 
experiencing  "snap-through"  buckling,  in  which  case  considerably  larger 
response  times  may  be  required.  A printout  frequency  of  once  every 
20  steps  is  usually  adequate  for  monitoring  the  response  time  history. 

The  program  checks  response  values  every  ten  time  steps. 

5.4  PROGRAM  OPERATION 

Two  versions  of  N0VA-2S  have  been  assembled  due  to  the  difference 
in  core  allocation  between  the  Control  Data  Corporation  (CDC)  6600  and 
the  CDC  176  computer  systems.  The  6600  version  has  smaller  dimensions 
and  only  uses  small  core  memory  (SCM) . The  176  version  has  somewhat 
larger  dimensions  in  DEPROP  and  makes  use  of  large  core  memory  (LCM)  by 
assigning  6 common  blocks  to  LEVEL  2 (see  table  12) . 
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Otherwise,  the  deck  structure  is  the  same  for  each  version,  using 
segmentation  as  in  the  past  (figure  43  and  table  22).  Logical  files 
TAPE5  and  TAPE6  are  used  for  input  and  output,  respectively,  and  file 
TAPE1  is  reserved  for  internal  use.  When  the  REFRA  near-ground  reflection 
(blast)  model  is  selected  and  ground  reflection  is  to  be  included  in  the 
analysis  (KB**1  ,KGRD«1) , the  REFRA  data  must  be  available  on  logical  file 
TAPE10.  This  data  is  unformatted  (binary)  with  record  type  S and  block 
type  C.  This  record-type-blocking  combination  is  compatible  with 
SCOPE  3.3  and  possesses  the  important  feature  that  the  system  copy 
utility  COPYBF  can  be  used  with  SCOPE  3.4  to  transfer  the  data  to  disk, 
tape,  etc.  When  used  with  the  NOVA  program  under  SCOPE  3.4,  however,  a 
FILE  card  is  required  prior  to  loading  to  specify  the  record  type  and 
block  type.  The  REFRA  routine  will  usually  operate  more  efficiently 
when  the  information  is  on  disk,  so  a transfer  to  disk  is  recommended 
whenever  possible. 

Using  the  FTN  compiler  on  the  176,  approximately  60,000g  cells  of 
SCM  and  15  seconds  of  CP  time  is  required,  specifying  the  fast  compile 
mode  (OPT-1).  , 

The  6600  version  of  the  program  requires  approximately  260,000- 

O 

cells  of  SCM  to  load;  the  176  version  needs  approximately  172,000g  cells 
of  SCM  and  255,000g  cells  of  LCM. 

5.5  NOVA-2LTS 

This  special  version  includes  the  modifications  made  in  previous 
versions  N0VA-2L  and  N0VA-2LT , where  the  "L"  refers  to  the  fact  that  all 
of  the  blast  and  aeronautical  loading  routines  have  been  replaced  by 
user-generated  functions  or  tape-supplied  data  (reference  2). 

No  other  changes  were  made  here,  except  to  adapt  the  up-to-date 
response  codes  in  N0VA-2S  to  this  special  deck.  Table  23  lists  all  the 
subroutines  in  the  deck,  and  figure  44  and  table  24  document  the  proper 
segmentation  directives.  Core  requirements  for  this  deck  are  approxi- 
mately the  same  as  for  N0VA-2S , except  for  an  additional  LCM  requirement 
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WFPR 

WFPKV  WFDZR  AIR 

Level  2 WFPKOD  WELL  WFPVZR  

WFDRMT  WFVRMT  WFZR  BISH 

SHOCK  WFPRMT  WFPKOP  SKIP 


Figure  43.  Segmentation  Tree  Structure  of  NOVA- 2 S 


TABLE  22 

SEGMENTATION  DIRECTIVES  FOR  NOVA-2S 


TWEE 

•'4  0 V A 

VOVA 

TNCLUDE 

LEVEL 

PIVJT. IODUM.sEC 

TPEE 

BLOCK 

BLOCK 

INCLUDE 

•VI  v.F'JE^SL.mOVSUN.PITC.bITEP 

TPEE 

OEPPDP 

T»EE 

UEPWOB 

OEPHOB 

INCLUDE 

OEFOWM.EB.FBCTl.PWINTl.WLAXF.STWr  P 

TPEE 

FPRES 

FPRES 

T MCL  une 

PFiJSF. , P .lU'^P 

TWEE 

wPWPS 

RPRE  S 

INCLUDE 

I VTSLO.PRE/V,  SETi*i.  POST,.!  , p Fi.ST  V3,  PDST/.a,  POST  v5.  P 

0ST*6, 

pos  rm 

LEVEL 

TWEE 

03ET1 

OSt  T 1 

INCLUDE 

D3ETP. LEGEND, QTSTEP 

TREE 

05ET3 

ose  rs 

INCLUDE 

BOLT 

TWEE 

STIFF 

STIFF 

INCLUDE 

M A r X I M 

TREE 

0EWV1 

OF  R V 1 

INCLUDE 

SIGna.SIGmAB.hIV.DEWVP.LISH.I.ISTP 

TREE 

WELAXP 

TREE 

CO^PI 

COmpi 

INCLUDE 

CUMPP. slay .STRV1 , VCS 

twee 

OAB 

DAB 

INCLUDE 

DPUW.EOUlLX,FBSET,FSUL,«tSO,RFRF  T.RLAXB.STWESX.STSET 

TREE 

CYCLE 

cycle 

INCLUDE 

E'lvJlLP, STRESS, FINAL 

TWEE 

CORSET 

COMSET 

INCLUDE 

READ  1 . TSTEP 

TREE 

level 

BL  AST  - ( XBL  AST  - ( TNTg,  TPI  '.T  ) , WEFWA  - f DP  It  , OPT?,  (>P  I <J) 

TWEE 

CStTuP 

CSF  TUP 

include 

1 vTl 

TREE 

SO.LVF 

TREE 

PRESS 

P»ESS 

INCLUDE 

I v TP 

TWEE 

H Y 0 W A 

HYDRA 

INCLUDE 

<A  r.Kh?,  I OPT  1 , I OPT?,  I OP  I i.  -'ELL.  Fr  RV'T  . -Fn  7.  W , >.F  P-v  3* , -„FP 

<00,  AIR./vFP*V.ftFPW,  yF2k,<*FPwmT,  *E vW^t , 'f  J H . amJC*  » a f '•»  0 S 

TWEE 

pjstap 

POST  AP 

INCL  IDE 

aDVANC.wEAU.SkIP.BISh 

VOVA 

GLOBAL 

FIRST.C  'LiVA,OV(jvA,rTL<“^Al/E 

HYOR  A 

GL  -1BAL 

CUVSTC.SCAL ft 

*F’WES 

global 

P« 1-SAVE 

WFF  RA 

global 

REFpaC-SA VF 

OFPPu« 

GLOBAL 

Blkp.blk  s,bl<(?-SAvE 

Dfc  PROP 

GL  ORAL 

CBLK1,rBLKP,CBLK3,CPL<-‘.C^LKS,rBLK7,CBL,,P.r“L<P.C-'L<1 

0 » CBLK 1 1 ,CRL<1  3 

•CBL<15,CPL*  1 7,CBLANK-SAVt 

OEWVl 

global 

CBLKW.CBLx 1 w-SAVE 

PEL  a XP 

global 

END 

CBLkip-Savp  »„»TTTYmCIICABiil 

nova  ihis  ?ao®  is  BEST  VJM'1  ^ 

OOPY  TUKlLLSHBC  I®  PftQ  - 
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TABLE  23 


LIST  OF  SUBPROGRAMS  FOR  NOVA-2LTS  (NOVA-2L) 


NOVA 

DEPROP 

DEPROB 

NOVA 

deprop 

DEPROB 

SEC 

BOLT 

COMP1 

RITER 

DERV1 

COMP  2 

CSETUP 

DERV2 

COMSET 

PINIT 

DSET1 

CYCLE 

SOLVE 

DSET2 

DAB 

INTI 

DSET3 

DEFORM 

PRESS 

DTSTEP 

DPUR 

HIM 

EQUILP 

LEGEND 

EQUILX 

LIST1 

FB 

LIST2 

FBCTL 

MATXIN 

FBSET 

RELAXP 

FINAL 

SIGMA 

FSOL 

SIGMAB 

PRINT1 

STIFF 

READ1 

RESD 

RESET 

RLAXB 

RLAXF 

SLAY 

STRESS 

STRESX 

STRN1 

STRN2 

STSET 

TSTEP 

VCS 
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Level  2 


Figure  '44.  Segmentation  Tree  Structure  of  Modified  N0VA-2LTS  (N0VA-2L) 


TABLE  24 

SEGMENTATION  DIRECTIVES  FOR  NOVA-2LTS  (NOVA-2L) 


TK  EE 

111  V A 

nova 

INCLUDE 

LEVEL 

PI  JIT. SEC.  1 00 UM 

THEE 

PI  TEH 

THFE 

OEPHiJP 

DEP«OP 

INCLUDE 

DEEUHvi.FR.FbCTL.  PPINTI  , 'JL  WF,  srP'  ? 

THEE 

LEVEL 

OEPHQP 

TBfeE 

D St  T 1 

dsf  r i 

IMCL JOE 

D5ETP. LEGEND, OTSTEP 

T«E£ 

OSEM 

osem 

INCLUDE 

bolt 

TBEE 

STIFF 

STIFF 

INCLUDE 

M A T X I M 

THEE 

HELA  kp 

thee 

OtHVI 

DEPV1 

INCLUDE 

Hl‘A,oEPv?,L  IST1  .LlSTE.blGJA,  SIRRAH 

r pee 

CO^lPI 

COfJ|Pl 

INCLUDE 

ClJMPP.  SL  ay,  STHnI  . VCS 

TPgE 

0 A * 

0 4* 

INCLUDE 

DPuH.EOUILX.FHSET.FSuL.HESD.PFFE T.BLAXB.STPESX.SrSET 

TWEE 

CYCLE 

CYCLE 

include 

EQJILP.STH6SS. FINAL 

TBEE 

COMSET 

comsf r 

I INCLUDE 

level 

HEAD  t . TSTEP 

TBEE 

C SETUP 

CSETijP 

INCLUDE 

i»rt 

THEE 

P H t S S 

THEE 

SOLVE 

NOVA 

GLOBAL 

CNQVA.CLOAO.CBLM-SAVE 

NOV  * 

GLOBAL 

CD-41  .CUV? 

OEPPOP 

GLOBAL 

C0LKP.CBLK1.  CDLK<I.CBLS,3.C*L*7.r<L'>.A.C*L!'O.C'<L'<l'),Cr<i_<; 

. 1 1,C*LKH,CBL«1S 

. CBL< 1 7 .CHLANK-SAVE 

PEL  4 y P 

Gl OB  A L 

COLK 1 P-SA  VE 

DEPV1 

global 

CBLKB. CBLS 1 S - 5 6 V E 

OEP»OS 

Gl.  OPAL 

E 4.0 

RLKH. ELM. *L*b- SAVE 

IBIS  PACE  IS  BEST  QUALITY  FRACTICABUI 
t&M  OCffY  rURJUSHED  TO  DDG  
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when  the  tape  option  is  used  for  supplying  pressure  data,  as  before. 
Appendix  A contains  a card  listing  of  the  appropriate  UPDATE  changes  to 
transform  N0VA-2S  into  N0VA-2LTS. 


5.6  EXAMPLE  PROBLEM 

An  example  problem  is  presented  in  this  section  to  provide  the  user 
with  a test  case  for  exercising  N0VA-2S,  and  to  also  indicate  the 
modelling  technique  used  in  analyzing  stiffened  panels. 

The  panel  modelled  is  the  fin  panel  on  the  vertical  tail  of 
the  B-52  discussed  in  section  4.  Figure  45  lists  the  input  data 
for  an  inelastic  response  run  (KDAM  * 101) . The  model  contains  one 
inner  gamma  stiffener  located  at  the  13th  beta  position,  or 
y ■ ^ “ 10.9  inches,  of  a clamped  panel  exibiting 

symmetry  in  both  coordinate  directions. 

Figure  46  contains  the  time-history  printout  at  time  1.75  milli- 
seconds, approximately  the  time  of  peak  response.  A summary  of  the  run, 
including  the  maximum  response  compared  with  allowables,  follows  the 
response  output  and  is  shown  in  figure  47.  In  this  case,  a tensile 
strain  at  the  clamped  edge  of  the  stiffener  produced  a CRIT  of  0.539. 

The  maximum  CRIT  of  0.445  in  the  panel  was  also  due  to  a tensile  strain 
at  the  same  edge,  at  y ■ 11.81  inches. 

Computer  time  for  this  relatively  large  structural  model  (30  modes 
and  361  spatial  points)  was  403  cp  seconds  and  2 EC  seconds  on  the 
CYBER  176. 


< 
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8-5?H  FIN 

PANEL 

1 

0 

1 000. 

525. 

1000. 

0.0 

0.0 

0.0 

15 

15 

7000. 

2 

0 

101 

0 

.5 

2 

2 

-426. 

0.0 

-1299. 

1 110. 

-783. 

0.0 

-1407 

1110. 

2 

2 

-1493.5 

0.0 

-1774.5 

309.81 

-1825.8 

0.0 

-1859.3 

309.81 

2 

2 

-1493.5 

0.0 

-1731.5 

278.8 

-1795.5 

0.0 

-1851.7 

278.8 

9 

-45. 

-1935.87 

0.0 

-89.1 

-75. 

30. 

-202.3 

64.4 

-538. 

69.1 

-847. 

69.1 

- 1 2°  7 . 

69.1 

-1307. 

66.2 

-1417. 

59.5 

-1462. 

56.71 

-1846.5 

44.93 

4 

1 

-1758. 

0.0 

35.0 

0.0 

7 

7 

1 

3 

13 

1 

3 

1 3 

0 

0 

0 

1 

19 

5 

11 

5 

13 

7 

9 

7 

11 

7 

1 3 

IBIS  PACE  IS  BEST  QUALITY  PRACTICABLE 
rau*  COPY  JfUKWISHJSD  TO  DOC 


0.0 

0.0 


(BLANK  ) 


1 

19  19  5 

5 7 9 1 1 

5 7 9 1 1 

8 1 0 1 


Figure  45.  Example  Problem  Input  Card  Listing 
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0 

9 

7 

9 

fHf  8 P48*  is  BEST  QUALITY  PRACTICABLE 

WW#  t>Ur*  FURNISHED  TO  HOC 

9 

i 1 

9 

1 3 

1 1 

5 

1 1 

7 

1 1 

9 

1 1 

1 1 

1 1 

1 3 

1 3 

5 

t 3 

7 

1 3 

9 

13 

1 1 

1 3 

1 3 

1 1 

1 

7 

1 

1 i 

1 

19 

5 

1 3 

9 

13 

1 3 

1 3 

1 7 

1 3 

19 

1 

1 9 

7 

1 9 

1 3 

1 9 

t 9 

l 

42.0 

33.7 

,03a 

.259 

£-3  1.05 

E 7 

.33 

5 o o n o . 

1 3 

1 .29 

P 5 .15 

-1.0 

1 . 04 

£7  9.0 

Pfe  .259 

£-3  70500.  70500. 

5.Q 

M .1 

0 

0 

.000393 

0.0 

.0*4 

1 .97? 

.2945 

.120 

.535 

.120 

. 7555 

.12" 

.9*6 

. 1 20 

1.05 

1.3® 

0.0 

0.0 

0.0 

o.o 

0.0  0.0 

0.0 

0.0 

0.0 

O.o 

0.0  O.o 

O.o 

0.0 

o.O 

O.o 

o.o  n.o 

0.0 

O.o 

o.o 

O.o 

■).!<  0.0 

0.0 

o . ;> 

o.O 

0.0 

o.o  0.0 

O.o 

o 0 

o.O 

0.0 

o.o  o.o 

0.0 

O.f 

0 . 0 

0.'' 

0 . o.o 

0.0 

.0 . 0 

0.0 

O.o 

9.0  0.0 

0.0 

O.o 

o.O 

O.o 

0.0  0.0 

0.0 

5. 

E mh  2.0 

F.-S  50. 

(9L ) 

Figure  45.  (Continued) 


131 


17S«lOOnE-oa  SEC 


mis  PACE  IS  BEST  QUALITY  FRACTICABU 
H» ■ OQPY  JTUR1USHJSD  TO  DOC _ 


** 

o 

e 

Oj 

o 

*- 

rv 

rv 

rv 

nr. 

rv 

— 

rv 

rv 

w4 

^4 

rv 

rv 

*4 

rv 

rv 

nO 

rv 

c 

c 

c 

x 

c 

c 

<z 

o 

o 

c 

o 

o 

c. 

O 

o 

x 

o 

c 

c 

o 

o 

X 

c 

c 

o 

c 

o 

c 

o 

c 

♦ 

♦ 

*4- 

i 

1 

1 

i 

♦ 

• 

I 

• 

i 

• 

1 

• 

i 

1 

i 

i 

1 

t 

• 

» 

• 

i 

• 

i 

i 

• 

i 

U) 

La. 

u. 

x 

U.’ 

UJ 

x 

u. 

a) 

X 

X! 

X 

X 

IX 

UJ 

X 

X 

UJ 

X 

X 

UJ 

UJ 

UJ 

UJ 

X* 

UJ 

U' 

UJ 

UJ 

UJ 

ifl 

3 

3 

3 

X 

X 

no 

3 

3 

C 

*K» 

X 

X 

X 

•■4 

X 

X 

a 

fx. 

X 

rv 

O' 

fx- 

in 

a- 

a* 

X 

x 

o* 

f*. 

a- 

CT 

in 

3 

X 

r*~ 

(T 

rO 

*-* 

nO 

w** 

K> 

rv 

rv 

c 

X 

X- 

O' 

3 

rv 

X 

X 

c 

c 

X 

in 

a 

* 

cr 

a. 

X 

fN. 

rv 

<V 

—* 

X 

K> 

r*» 

X 

in 

3 

O* 

#-x 

•-» 

rv 

o 

<T 

O' 

X 

3 

X 

"O 

ST 

X 

rx- 

c 

in 

f*. 

rv 

e 

cr 

nO 

O 

o 

CT 

0 

rv 

rv 

o- 

a 

in 

fx. 

a 

3 

ro 

X 

X 

X 

rv 

O 

rv 

X 

re 

•— 

c 

X 

IT 

•— 

3 

«— 

V 

X 

X 

X 

rv 

y 

rv 

rv 

^x 

3 

•■X 

3 

hx 

ro 

•— 

f«C 

X 

3 

X 

3 

rv 

no 

X 

r*» 

3 

X 

X 

X 

no 

•- 

X 

a 

•-* 

X 

rv 

•3 

o 

rv 

in 

X 

3 

X 

X 

c 

rv 

o 

rx 

cr 

— 

a. 

K) 

— 

no 

3 

3 

ro 

rv 

x- 

X 

<n 

X 

V 

**r 

•V 

K, 

rv 

Oj 

— 

<V 

rv 

X 

— 

• 

» 

I 

ft 

• 

o 

rv 

v> 

^4 

<-U 

'V 

r\» 

rv 

3 

a 

rv 

-f 

rV 

_* 

*n 

rv 

3 

"O 

v> 

v; 

•n 

V> 

OJ 

V 

•r* 

rv 

•O 

rv 

•o 

o 

x 

X 

T- 

o 

c 

c 

C1 

X 

O 

“ 

c 

3 

o 

•3 

cr 

~ 

O 

O 

o 

3 

3 

O 

—> 

c 

o 

o 

♦ 

i 

• 

• 

• 

« 

< 

t 

• 

• 

• 

• 

1 

1 

i 

ft 

i 

ft 

ft 

ft 

• 

ft 

1 

ft 

1 

ft 

• 

1 

• 

ft 

uJ 

UJ 

JLl 

X 

OJ 

X 

jJ 

jj 

aJ 

UJ 

aJ 

Jj 

X 

JJ 

jj 

Xl 

JLj 

X 

X 

Jj 

J-» 

X 

Jj 

uJ 

uJ 

a.’ 

X 

X 

X 

_u 

X 

■X*. 

o 

n 

rv 

n 

T 

c 

fx* 

Ti 

'O 

»Xfc 

in 

c 

3 

3 

»— * 

X 

T 

3 

c 

r 

J* 

x> 

X 

V 

X 

*>r- 

X 

O 

w— 

LT 

X 

IT 

rv 

ro 

c 

c 

r^* 

•-» 

X 

C 

X 

in 

3 

IT 

rv 

X' 

o 

X 

r*x 

3 

X 

X 

c 

X 

rv 

> 

a- 

fx. 

c y 

no 

•— i 

O' 

px. 

in 

rv 

X 

o 

V*'. 

X 

o 

O 

rv 

a 

X 

X 

•— 

X 

w~ 

fx- 

X 

r\ 

rv 

o 

3 

•M 

o 

n*'- 

rv 

m 

rv 

rO 

a 

O' 

X 

X 

«o 

m 

rv 

a 

X 

C 

X 

i r 

X 

<T 

r*te 

nr, 

'O 

X 

X 

X 

X 

c; 

•— 

no 

X 

fO 

•-« 

rv 

c. 

O 

O' 

a 

i r 

O 

aC 

X 

in 

ro 

X 

o 

x 

n 

a 

in 

X 

rv 

«n 

a 

in 

K 

c 

•- 

3 

K*. 

rv 

X 

•— 

in 

tn 

3 

fx- 

rv 

c 

3 

3 

n 

r\ 

c 

3 

<r 

X 

•— 

X 

3 

c 

a 

s 

— 

rv 

*o 

X 

rfi 

X 

rv 

X 

n 

•O 

— 

— 

rv 

C* 

IT 

X 

3 

— 

— 

rv 

— 

r*» 

— 

3 

— 

3 

1 

1 

t 

i 

1 

» 

ft 

1 

1 

ft 

ft 

ft 

1 

rv 

rv 

— 

no 

rv 

rv 

rv 

rv 

ro 

K". 

V 

•— 

rv 

rv 

rv 

rv 

«— 

rv 

O* 

rv 

rv 

H 

V 

n 

«o 

o 

V 

c 

c 

o 

C 

W 

c 

C* 

X 

o 

C 

o 

O 

e 

o 

c 

X 

c 

X 

X 

W 

C 

X 

x 

X 

X 

C 

» 

1 

ft 

ft 

ft 

1 

i 

• 

• 

ft 

ft 

» 

ft 

ft 

l 

• 

ft 

ft 

ft 

t 

ft 

ft 

ft 

1 

ft 

• 

• 

ft 

ft 

ft 

UJ 

UJ 

u. 

u. 

UJ 

X 

X 

a. 

J- 

u. 

X 

X 

X 

X 

X 

X 

u. 

U' 

X 

X 

Ju 

X 

X 

U.‘ 

X1 

X' 

X 

•J- 

u 

X 

^r, 

•- 

c. 

C 

X 

ro 

X 

in 

tr> 

cr 

X 

rv 

X« 

<J 

r- 

rv 

3 

cr 

r\ 

3 

V 

3 

r 

or 

c 

s 

•— 

'X 

C 

3 

n 

O 

*o 

Xl 

V-4 

r*- 

O' 

•n 

X 

cr 

*** 

3 

X 

r«* 

r^ 

cr 

O 

rv 

N. 

X 

•o 

C 

X 

_rs 

~ 

rv 

IT 

r- 

OJ 

c 

in 

»0 

*r» 

T 

CT 

o 

rv 

OJ 

3 

X 

3 

(T 

3 

o 

c 

ro 

o. 

•o 

m 

n 

X 

X 

o 

3 

— 

3 

3 

n 

3 

iT 

fO 

-o 

X 

3 

X 

X 

m 

3 

3 

— • 

— 

3 

c 

•— 

a 

3 

■n 

v, 

no 

X 

3* 

3 

•o 

r- 

3 

no 

nr. 

rv 

n 

3 

?r 

•n 

c 

— 

3 

3 

at 

rv 

X 

X 

rv 

X 

•— 

rv 

OJ 

a 

O' 

re 

•O 

•- 

•n 

NO 

e 

PN. 

no 

r 

rv 

c* 

3 

«— 

O 

X- 

o 

X 

a 

X 

7 

IT 

rv 

»■« 

X 

rv 

rv 

"O 

X 

X’ 

cr 

no 

D 

3 

•O 

OJ 

CC 

n 

— 

in 

<n 

f«r* 

in 

• 

I 

l 

1 

1 

t 

9 

» 

1 

1 

ft 

1 

• 

« 

i 

i 

— «0  J no  3 


K J".  n 3 — — •* 


5 


II  C 


/ / /'/  a O' 


132 


u. 


Figure  46.  Example  Problem  Output  at  1.75  Milliseconds 
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Figure  47.  Summary  Output  for  Example  Problem 
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INTO  NOVA-2LTS  (NOVA-2L) . 
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I 


n n o n n o 


• I n d U T I'lPLTS  XBIS  PACE  IS  BEST  QUALITY  PRACTICABLE 

• RtlROECK  BLOC*  ntOE  OOPY  FURBISHED  TO  DDC 

• piirofck  I 0 0 U M 

• P.JROFCK  MN.RITC 
•RURDFCk  I j T P 

• O’lkDECK  BLAST. TPINT 

• PiJRDFCK  TNTa.PnST.vi7 
•PUPOECK  PRErt.vPRES 

• PIJBGF  BLOCK 
•PURGE  10DUB 
•PURGE  NIM.RITC 
•PURGE  1NTP 
•PURGE  BLAST. TPINT 

• PURGE  I'lTa.POSTA7 
•PURGE  PWFW. NP9ES 
•AUOEILE  INPUT, CWOVA 
•COWOECK  CLOAD 

COMMON  / C L 0 A l)  / PPl,PPO,TTO,TPPlME,AA,A'lN,OTTt,OTTO,A?, 

1 JL.MTIME.MLUAO,PT(ao)  .TTfPO)  , K T I ME  ( 1 0 , 1 0 ) , L T I ME  f u 1)  , I SP  ( <1  0 ) , 
a jLB(an  ,nps,  r tp (b,  io,  i o)  ,p*t if>,  to,  t o) , npx , npv,  *p(pp)  , rptagj , 

S I X I (?3)  . J Y J ( £3)  ,JLT  (10, 10)  , PR  T T ( 10, 1 0)  « L)  X 1 (?3)  ,0V)  (?3)  * 

« NGSUM,  del.  MAXOa,  'AAXO.PS 

•COnOEfK  ru«i 

common  /com/  Piaa.n 
level  ?,  p 

• CUMOECK  CO^P 

ci.xmom  /com?/  g(a?,n 
LEVEL  ?.  ') 

•PUPTECK  NOVA 
•PURGE  r TV  A 
•AHUF1LE  TNPUT.CBLANa 

• 0 F C K nova 

program  NOVA  (InPUT.OUTPLiT,  TAPFSririPuT,  TAPEfesOliTPtiT,  TAPElsbl3, 

i tapfio) 

This  IS  The  NCVA-PlTS VERSION  OE  nova.  The  AfcROOTNAMC  AND  blast 
routines  are  REPLACED  by  USER-DESIGNATE?)  PRESSURE  E'.JNCT  I C'.-jS. 

PB  I VISION  BUS  A L S'1  BEEN  NAPE  POP  BEADING  PRESSURE  LATa  Fhijm  TAPE. 

JUNE,  1978. 

• Call  Clfun 
•call  CNOVA 

c 

i format (Miai 
3 forma  r ( fe i a . 1 1 

» Format  ( a -1 1 a ) 

HCASe  s f* 

I .OuT  s 1 

k F » 3 1 . f. 

"PA')(S,n  '.CASES 
1)<»  -E'.ofS.n  I T1  TLC(  1 I,  T»1  ,?r  1 
•r  a sf  s • r a se  ♦ t 

*k  *»W  3 *) 

vtrial  a n 

--»o  <*.:  j jrn*pt*  Type,  a » •,  kos,  • pg'G 
if  ( < >5.E  . n «V  • s » 
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* 

( 


ir  . Ibis  pice  is  best  quality  praciicabui 

if  («qIm.lT.8)  „CHPT  s 1 FBCkB  OOPY  FURBISHED  1*0  DDC  

IF  (OAM.GT.2)  K U A K = K n A M - 10  0 

IF  fK0AM.LT.?)  HEAD  (5.PJ  FOAM 
IF  ( iriO'JT.EQ.O)  GO  TO  1400 
«9I  TE  (b,  AOOOj  (UlL£(I),l!l,?0) 

IF  (KTYPF.LT .b.OR.<TYPE.GT.7)  ICO*?  = b 
IF  (ICnMP.FQ.?)  **ITE  ( b » 500 0 ) 

IF  fNCHPT.EO.lJ  «KITE  (5,3100) 

IF  (MChPT.EQ.O)  I TF  ( b , 3 2 0 0 ) 

IFfKOAM.fO.O)  *RI  TE  (b.  3.300)  P9AM 
TF(KOAM.eQ.l)  WKlfF(b.3a('O1P0AM 

GO  TO  (300, 400.500, 500,700,800, 900, 1010, 1000, ] 020J ,K T YPE 
300  Aft  I TE  f b , 350  0 ) 

GO  TO  1050 
aOO  ARlTECb. 3bOO) 

GO  TO  1050 
500  *«ITE(b,37O0) 

GO  TO  1050 
bOO  *SITE(b.3800) 

GO  TO  1050 
700  .<.RITE(b,  3900  ) 

GO  TO  1050 
800  oRITE (b. 4000) 

GO  TO  1050 
900  NRI TE (b, 4 1 00  ) 

GO  TO  1050 
1000  *»IT6(b.a?O0) 

GO  TO  1050 
1010  A»ITE(b.a700) 

GO  TO  1050 
1 0?0  a«ITE (b,4800) 

1050  GO  TO  Cl  100, 1200.  1300) , KOS 
l 1 00  rtFITE(b.«300) 

GO  TO  1400 
1200  At?  I TE  (b,  4400  ) 

GO  10  1400 
l 400  Ad  1 Tt ( b, 4500) 

1400  NC ALL  = 2 

IF (KTYPE.GT.5)  CALL  0EP90B 
IF  (KTYPE.LT.bJ  CALL  OEPWOP 
IF  fKE99.GT.0J  GO  TO  lr>O0 
NC  ALL  = 1 
CALL  PINITfO) 

IF (Kl YPF.GT.5)  CALL  oeprob 
IF  (KTYPE.LT.bJ  CALL  OEPwrp 
IF  (KOS.ES.  1 ) G'J  Tu  IstiO 
IF (KEH9.GT.0)  Gu  TO  1 bOO 
NC ALL  = 0 
K OK  = 0 
CALL  PIN  II  ( 1 1 
IF  (KEHP.GT.O)  GU  U 170'. 

RT9I AL (11  = 1 .0 
1500  uTRIAl  = NTRIAL  ♦ l 

IF  (KOAM.LT. 2)  ,KlTE(h.4A0Pi  NC  ASt , I *L  , 9 T * I Al  ( 1) 

TF  (KIYPE.GT.5J  CALL  OEPFl’rt 
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1 80  0 
1700 

C 

?o,ir> 
3100 
3P00 
3 300 
3000 

3500 
3800 
3700 
3800 
3900 
<100  0 
01  00 
«poo 

0 300 

oooo 

0*300 

U800 

0700 

0800 

soon 


IBIS  PACE  IS  BEST  QUALITY  PRACTICABL* 

IF  (KrvPF.LT. 6)  C4LL  OFPm'P  ERO*  COPY  FURHISRED  TO  DDC 

IF (kFWR.NF.O)  GO  TO  lsno 
IF  (NCHPT.EO.O ) GO  TO  m«0 
CALL  RITFR  (COU.kNUI.,.'  TRIAL. «,*0><) 

IF (KOK .EG. 0)  Go  lu  1500 
JFfMCASF.LT.NCASfcS)  GO  10  100 
STOP 

FORMAT  11H1.30X,  17HN  0 V A - P L T S//lX,ao*0) 

FORMAT  (10H  ITtKAriJi'i  RllFi ) 

FORMAT  (??h  RESPONSE  RUM  ONLY*  NO  ITERATION) 

FCiRmaT  (oPh  NO  J a m a 3 F LEVEL.  PROBABILITY  OF  EXCEEDING,  F8.3) 

FORMAT  (5?m  CATASTROPHIC  OAmaGF  LEvEL.  “RDBArilLlTY  OF  EXCEEDING, 

1 F 8 . 3 ) 

FORMAT  (P8H0STNGLE-L4YFR  n£TAL  PANEL  ) 

FORMAT  ( ?OHOSI.MGLt-L  A YEP  PLASTIC  PA  Mtl.  ) 

FORMAT  ( ?SH0m0 xE YC0M8  mfTal  PAMEL  ) 

FORMAT  (P7HOhOn£YCCH8  PLASTIC  PANEL  ) 

FORMAT  (PRhomuL I I-LA yER  PLASTIC  PANEL  ) 

FORMAT  (30M0M£TAL  stringer  OR  longeron  1 
FORMAT  (15H0MPTAL  FRAME  ) 

FORMAT  (16H0PLAST1C  RING  ) 

FORMAT  (PlrtO STATIC  SOLUTION  ONLY) 

format  (PPPOOYYaMIC  response  ONLY) 

FORMAT  ( 37mOST  A T I C SOLUTION  A NO  DYNAMIC  RESPONSE) 

FORMAT  (IPHtCASE  NUMBER  IP/ 

Mum  TRIAL  NUMBER  13,  103,  18H  RANGE  FACTOR  = Ela.8) 

F0R«AT(1 IhOmETAl  RING) 

FORMA  T ( 1 3M0OI8  e»UC*Ll  vG) 

Format  (8RM0STRuC  T'jRAL  ELEMENT  DOES  DERIVE  ADDITIONAL  SUPPORT  Fr  O 
1 F iJSEL  AGE  SKTN) 

END 
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♦ piiwofc*  PI'i IT 

• Pitt-Rc  PI 'I  IT 
•AODFILE  I'jp  jr.CSfcTUP 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLI 
FRO*  QQrY  JCURMSfiED  TO  DOC 


• OFC*  PI ''*11 

Sl"»M1uriME  PIWITC'O 

• Call  C n ' j v a 

•CALL  CLuAO 

• CALL  C^U*<  1 
•CALL  Cnvi 
•CALL  CO^P 

0 I ME  VS  I O’J  10(100).  00(100)  , A I ( 1 004)  , MOPDt* (40) 

EQUIVALENCE  (10(1). 00(D) 

dimension  TTB  ( 1«,  4D  .PPTO  ( l «.ai ) , SP  («D  . 5PS  (an 
EQUIVALENCE  (Prtl (1, 1 , 1) ,PRT»(1 . 1) ) , (TTP(1 , 1 . 1 ) ,TTB(1 , 1 ) ) 
OATA  TPO/lOhFPFPFFFFFF/ 

C 

IF ( *.FQ.l 1 GO  ro  800 
PS  S 0.0 

IF (KQS.FQ.8)  GO  TU  150 
C 

C STATIC 

r 

READ (5.8000)  PS 
*'P  I TE  ( 6 » 8800 ) PS 
MU=  t 
PPPsPS 

IF  UTYPE.LT.b)  GJ  10  150 
IF  (KTYPF.LT.10)  GO  TO  50 
1)0  50  Isl.NMASS 
50  PB(  I)  s 0. 

GO  TO  150 

5 0 00  100  1st, NM ASS 
100  P«(I)  s PS 
150  PETiJft'J 
C 

C DYNAMIC 

C 

POO  IF  (KOS.EQ.l  ) GJ  T IJ  400 
READ  (5.8050)  mLJAJ 
rtPITE  (5.8400)  i J L U A 0 
GO  TO  (850, 500,  *>90, 6000)  , *LO*C 
P 5 0 PFAO(5,POnO)  PP1  ,PP).  TTO,  TPPI -IE.  AA  , AWN 
n •<  1 TF  (b,  8'jOO)  PP1  ,PP),  TTU,  T»P  I*-’*.  A A,  A,!N 
NL'SI 

IF(TPPIMF.cG.p.O)  GO  TO  *on 

PFT  I vfFsPPO*  ( 1 .1)  - T T **F  / T I P ) • * A it: 

PPttl^F  s PPPIME*c*P(-AA*iFSI'-|F/TT<il 
TT1sTPPT'At*PP1  / (PPl-PPPT  af  ) 

OTTIsl.o/TTl 
500  OT  TOsi .0/ TTO 
A7sAA*()TTr, 

a 0 0 wtH'Sto 

C 

SO"  wFAl)  (5.80501  JflMe 

Pf  AO  (5,8100)  ( T! ( I ) ,PT(I) . Isl .- TTVF) 

*PDE  ( h , P5°  0 ) 1 1 Ia».  (TT(l).PT(l),  I s 1 ,M  T I 4 E ) 


t 

] 

( 

t 

i 


y 

\ ' 
k 

. 


1 
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IHIS  PAGE  IS  BEST  QUALITY  PRACTICABL* 

JL  s ? raoi  (Xffl  PUftUISHEC  ID  DOC  - 

‘ PF  TUP'1 

C 

800  IF  (KTYPF.GT.5)  GO  TO  1000 
C PAMFLS. 

PF  a n (5,2050)  ijPx.mpy 

'-P1TF.  (8.2700)  yPX.nIPY 

feao  fs.pnor.)  ( XP  ( 1) . 1 = 1 .-'•py) 

API  T F (b.MOO)  (XP(I),Tsl,MPX) 

&fao  (5,?noo)  c yp( J) . Jsi  ,><py) 

•y  fc  1 1 E (8.3200)  (VP(J1 . Jsl  ,NPV) 
a-wTTF  (8,3300) 

DO  *2u  isi.r.px 

PfAD  (5,2050)  (KTIAE(J.I). Jsl.NPY) 

8?n  T T F (h.?«on)  (iUI«E(J.  I)  . Jsl  .WP*) 

OR  8«0  Isl.NPX 
DO  840  Jsl.NPY 
NTIve  s KTIMt(J, I) 

PFAJ  (5. 2000)  (TTP(K,J.I),Ksl.MTHE) 

*W  1 TE  (8.  3600)  I,  J,  (TTP(K.  J.  I)  ,KS1  ,MTIA£) 

*»ITE  (8.2900) 

«FAD  (5,2000)  (P8TU.  J.D.Ksl.NTIPE) 

H u 0 aPJTE  ( h , 3000)  (PrtUK,  J.  n,K  = 1 ,NTIMF) 
c SPATIAL  TMEPPOLATI  JN-F  X T,?  A Pf)L  A THIN.  INDICES  A9fc  10*EW  BOIinO. 

oo  poo  isi.mgt 

no  PfeO  III  = i.nPx 
IF  (XPCII  I)  .GT.xSCI  ) ) r,n  TO  8«0 
880  CONTINUE 
III  = N°  X 

880  if  Tin.GT.i ) nr  = nr  - t 

DXim  = ( XG  ( I ) - XP(ITT))/(XP(III+1)  - XP(IID) 

POO  I x H 1 1 r HI 

DP  980  J = l.NriT 
on  0?0  JJJ  s 1,MPY 
TP  (YP(JJJ).GT.XB(J) ) GO  TO  940 
9P0  CUMTIMiJH 

jjj  = npv 

P40  T F CJJJ.GT.l)  JJJ  S JJJ  - 1 

dyitj)  = fxp(j)  - ypujj)  )/rypf  jjjm)  - yp(JJ.F)) 

QhO  JYJ(J)  = JJJ 

•III  S (I 

on  080  I s 1 , M P X 
no  o 80  J=1  ,1-iPY 
9«C  JLI  ( J.  I ) S 2 
WE.TI.IPM 
C dEA^S. 

1 000  WF40  (5#?05r.)  IMPS 
'•«  I T £ r 8.34110)  JPS 
w r a 0 (5.2000)  ( 0P(  I ) . Isl . * PR) 

•iwTTc  (8.3500)  (SP(  [ ) , Tsi  ,f>.PS) 

aPTTF  (8,3300)  > 

w F AO  (5,2080)  (Lll  iS(I),  Is)  ,'JPS) 

'•«W  T T £ (8,2800)  (L  T MF  ( I ) , Is)  , MPS)  L 

no  1200  Isl, MPS 
NTT^E  s L TIMF.  (I) 

WfcAO  (5,20on)  (!!)«,  I)  ,«sl  ,*JTlvF) 
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■nPITE  (6.3700)  I , ( r M (K  . T)  . K = 1 .►  TT'-'P  ) * £AeE  Is  BSsT  QUALITY  PRACnCAfiLl 

PFAO  (S.POOO)  ( P9  T H ( < , T ) » K s 1 , MT  I Mg  ) *t0Jm  COPY  fUKMlSRED  10  DDC  

« W T T F (6.3*300) 

IPO <1  /iP[TE  (6.3000)  (PPTrt(K.T) ,*=1 .NTIMfc) 

SPATIAL  pJTfctfPUL  A 1 I J"i  - FKTOAPPLATIOM.  IvOICES  APfe  LJ>«E*  lUUNO. 

1P5<>  00  1500  1=1, WPS 
PSP  = SPfl) 

II  = FSP  «•  .ooooi 
01  I = FSP  - floa  T ( I l ) 

SPSX  = 0.0 

IF  (Il.EQ.O)  GO  Tu  laOH 

on  1 300  j= i • 1 1 

l 300  SPSX  = SPSX  ♦ 0S00CJ1 

laoo  IF  UI.LT.  NMASa*l  ) 3PS(I)  = SPSX  0 I T *0  S>00  ( 1 1 ♦ 1) 

1500  CONTINUE. 

00  lt>oo  .1  = P , NMA  S3 
isoo  osno(J)  = osoo(J-i)  * osoo(j) 

00  1 H 5 0 I S 1 , MM ASS 

on  1700  in  = i.nps 

if  (spsnin .GT.osoomi  go  to  isoo 

1 700  CUNT  I 'IIJE 
II  T = NP.S 

MOO  IF  (II  t.GT.t)  111  = 111-1 

osood)  = (osouti)  - sPsciim/rspsdiiM)  - spsimn 

1«50  ISP(I)  = III 

IF  (NLOAO.EO.-a)  GO  TO  7100 
00  1900  T=l, NPa 
1900  JLB(T)  = ? 

«C  TUPS 

L 0 A 0 OPTION  a - TAPE  INPUT  FWO*  TAPEtO. 

PPOGPAMMEO  FOP  7600  ONLY. 

60  0 0 M AXO  = 1191U 

VAXO?  = 5957 
PE.AO  (5.P0001  TIM1.SKIP 
T I '>*?  = TIM!  + TSTOP  ♦ oflti* 

NSKlP  = SKIP  ♦ .0001 
PE  AO  (5 » ?050 ) NG AGE 

PE.  AO  ( 5 , ?050  ) (NOPOtPCIl  , 1 = 1 .NGAGF) 

**  I TE  (6,  3900)  TMl,  SSKIP.NGAGE,  (WOPOEP(I)  , 1 = 1 ,MGARt) 

MGSU'-t  = 0 
00  6050  1=1, MG AGE 

IF  (NOPl)FP(I)  . G T . 0 ) N6SUM  = NGSiJM  ♦ i 
ftOSO  CONTINUE 

TF  ( KT.rPF.GT. 5)  GO  TO  7 n f . o 
PA  nElS. 

PE  A 0 (5,?050)  NPX,  JPY 
apite  (6.P700)  ,px , >jpy 

if  (NPX.GT.  1 . a 4 U . u = Y . G T . 1 ) r.o  Tn  »*on 
P C A 0 (5,?O00)  IXPll)  , 1 = 1 .IPX ) 

«E4  o ( 5 . P 0 «)  C ) ( YP(  I I , 1 = 1 ,6PY  ) 
iPITE  (6,3100)  IXP(I1 , 1=1 ,NPX) 

«9TTE  (6. 3? 0 0)  ( Y P ( I ) , T s 1 , vP  Y ) 


IHIS  PACE  IS  BEST  QUALITY  PRACTlCAfiLI 
PBO*  OOPY  PURMISHKD  TO  DDC 


IP  ( N3X*NPY.El).  1 ) S'J  TO  7100 
MU  = 0 

IP  (NPX*NFY.NE.  JGMJM)  CO  TO  8600 

ip  (npx.eo.i)  gu  ro  6600 

00  6500  1=1. NGT 
on  6300  III»1,nPK 
IF  (XPmiJ.GT.XGlI))  CO  TO  6400 
3 MO  CONTINUE 
III  = IMP  X 

6400  IP  (UI.GT.n  III  = III  - 1 

pxim  = rxGm  - xp(iii)) /(xpini  + n 

6500  IXim  S III 
6600  IP  (NPY.EQ.l)  GQ  TO  7100 
00  6900  Jsi.Nfll 

00  6700  JJJ  S 1,'jPY 
IF  ( YP ( J J J ) . GT . X8  ( J ) ) GO  TO  6800 
6700  CONTINUE 
JJJ  = IMP  V 

6800  IF  (JJJ.GT.l)  JJJ  = JJJ  - l 

Clt(J)  = f X Q ( J ) - Y3(JJJ))/(YP(JJJ+1) 
oOOO  JYJfJ)  s JJJ 
GO  TO  7100 
5E AMS. 

7000  PE  AO  (5,3050)  MP3 
.'.RITE  (6.3400)  NPS 
PFAO  (5,3000)  (SP( I) . 1st ,MPS) 
aPTTE  (6,35001  OJP(I).Tsl.NPS) 

Jl  = 3 

nij  = l 

IF  (NPS.EG.l)  GO  ro  7100 
IF  ( NP5 . ME . NGSUM ) GO  TO  8500 

NU  = 0 

GO  TQ  1350 

7100  OFL  = 0. 

TCV  = t.E-h 
T T v*  1 s TI-l/TCV 
kg  S 0 

*8  s o 

K*K  r 0 
NT  IMF  s 0 

PUFFER  IN  (10,1)  (10(11.10(100)) 

IF  (UNIT ( 10)1  7300. 7300. 8309 

7300  8UFFE3  IN  (10,1)  (11(11,10(100)) 

IF  CJNITMO))  7350,8100,3300 
7360  IP  mm.EfJ.T30J  3J  TO  -700 
MAORIS  s 10(17) 
o-P'i  I NT  s 10(13) 

L R s m.vOROS *NPi)  I N r 
»:a5  = N«OP0S*MS8  IP 

< < S < K ♦ 1 

«G  S VO»l)fR  (<K ) 

IF  (KS.GT.O)  <88  s <<<  ♦ 1 
IL  s M#OROS 

7300  PUFFER  IM  (10,1)  (Al (11,*I(LR)) 

IF  (UMlT(lO))  7400,8100,6300 


XP(III) ) 


Y?( JJJ) ) 


7 4 i)0 

C 


7SoO 
7 *>00 


7700 
7 7 50 
7900 

7000 


7950 

*000 


C 


IP  (DEL.EO.O.)  OtL  = ( A T (1  ♦■M.vOKRS)  - A I ( 1 1)  •$«  I **  TC  W 
LOCATE  FI ‘•'ST  TI>ic. 

00  7500  JsIL.L*. 

IP  ( A I { I-n  ,GE.  T I ^ 1 ) GO  TO  7 SO  o 
CONTI  NUF 
GO  TO  7^00 

TP  (KG.GT.O)  P ( < G , 1 ) = AI(I) 

T 1 = AI(I-1)*TCV 

IF  (NTIMC.EQ.O)  >» T I = dlMja-TIMl  *TCV) /DEL  ♦ 2 


IP  <NTIwE.GT.MAXO)  30  TO  9000 
J = t 

IF  fKG.EQ.O)  GO  TO  7900 
TL  r NaS  - L*  ♦ I 
IP  CIL.GT.O)  GO  TO  7900 
IL  : I ♦ N*3 
00  7700  IsTL.L*. n*S 
I A s I 
J s J ♦ t 

IF  (J.GT.NTIME)  GO  TO  7700 
P(KG.J)  = Aid) 

CONTINUE 

IL  = NAS  - L9  ♦ IX 

BUFFER  IN  (10,1)  Uim.AHLS)) 

IF  (UNITdO))  7900,9100.9300 

IF  (KG.EQ.O)  GO  TO  7900 

IF  ( J . GE . N T I Mg ) 30  TO  7 000 

00  9000  I r I L , L H , N a S 

T * = I 

J = I ♦ 1 

IF  (J.GT.NTIME)  50  TO  oooo 
IP  ( J.GT.maxO?)  GO  TO  7050 
P ( » G . J ) = AKI) 

GO  TO  9000 

Q(KG,.I-mAx02)  = AI(I) 

CONTINUE 

GO  TO  7750 

END  OF  GAGE  DATA. 
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ODD 


8100  IF  ( N03IJG.GT.0)  At?  I T E (6.UP00)  KK  , K 5 , |j£L  , T 1 , N T I 8F. , 

1 mao?D3. t'jpoi  ?r,  ( nci ) . tst  ,<?i 

IF  (KKK  .1  I .NGSU’-I)  GJ  TO  7 ? 0 0 

c 0 4 t ft  =?  e a 0 . 

820  0 K e A'  I ?0  10 

a W I T 6 (6,9800)  JEL.  MpF 
IF  ( mOB'JG.LT.?)  GO  TO  8?50 

no  8220  isi,imgsom 

Ml  = MptO  (NTplE,  M4X0?) 

MP  s NTTME  - 84X02 

A i?  I T £ (h,il900)  1 , (FfI..n,Jsl  . M T 1 ) 

IF  (.NT2.GT.0J  V-?ITE  (6.U900)  I . (Q  ( I , J ) , J si , wT? ) 
8PPO  CONTINUE 
32S0  BETl.tPN 


EPW02S. 


3 400 

API  TE 
GO  TO 

(6.9300) 

8900 

3aoo 

APT  T E 
GO  TO 

(b.aaoo) 

8900 

84X0 

6S00 

•API  TE 
GO  TO 

(6.0500) 

8900 

NPS,  '45SU* 

3600 

•1 « T T E 
GO  TO 

(b.ttSOO) 
890  0 

iiP  X , ^Pf  .NGSUM 

3700 

a F T T £ 
GO  TO 

(6,0700) 

8900 

*<, 4i 4GF 

3800 

aPITE 

(6.9000) 

8 90  0 

K E °R  s 
PET'JP'I 

? 

c 

2000  FORMAT (6F12. I ) 

20S0  Fot?«ftT  (61121 
21  00  FOQvaT  (2F 12. 1) 

2P00  FOB^AT  (2UH0ST4  I IC  P2E58HBE.  PST  s E15.6) 
2(0  0 Fij^viftT  (2(H007-gft  MC  L 0 4 0 C^NSTa.'ITS/ 


1 1 H 

PP  1 

r 

El  5.6/ 

1 1 H 

PPG 

s 

£15.6/ 

1 1H 

TTO 

s 

El  5.6/ 

1 1H 

TPM 

s 

El  5.6/ 

1 IM 

44 

s 

€15.6/ 

1 1 H 

4 f ■ “1 

s 

El S.6) 

2000  FOUvtftT  (21  M00»  J4'a1C  LOftH  option  la) 

2S00  F 0 B n a T ( 1 OHONtriHEw  JF  TP'ES  = Iu/28h  THfc,  SEC 
1 (PEIS. Si) 

270  r FOSuat  ( 2«H0.'iu  -Mt1?  Jr  I_nat;  STATIONS/ 


1 1 ph 

P8  00  F 0 9 » T 
P 90  0 FOBmA I 
4 0 0 0 F 0 * v ft  T 
(100  F 0 ri  ‘A A T 
(200  FlOivtal 
4(00  F OB  w ft  T 
( a 0 0 F 0 w v A T 
(SO  0 FO»vaT 


mp*  s i ( / 1 p h iPr  = MI 

( S X , 1 0 l 3 ) 

(i»h  p-?fso!<?£s  f»sn  = ) 

(SX.SEls.o) 

(19h0a-PuS1  1 1 1 <S  (pi)  s/(sx,sf  | 5 . S ) ) 
(2SH0Y-PUSIT1  j JS  IPI  JB  OF  G J s/(5x,5El5.6) ) 
(2frH0MoA8E«<  JF  T ft  a L E f \ T «?  t E S =1 
(panO'JUSrtt'B  J"  Lf'ftn  S Tft  TTONS/  1 ?.-i  nPS  = 
( 2 S h 0 E 4 6 u ■?  * 4 F \ T POSITIONS  : ✓ ( j*  . 5«i  1 5.  6 I ) 


PhFSSOhE# 


1 yj 
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r S I / 


1 OX  , hHf.  P X 


! 5 , 5 x , 6HnPy 


I 

/ 

( 

A 


? 


) 


format  ( t <ihoti<*'EcJ  ( j fE c » 


I 3/ 


i f5x,fc.ei5.*n 

3700  FORMAT  ClUHOTivtS  ( 5E  C 1 , 1 0 X . PH  I = 1 3 / ( 5 < . 6t 1 5 . ft ) ) 

3600  FI|B«»T  (19h  P9£:>ou9E3  (SFC1  s ) 

39  0 0 FORMAT  (9H0T4Pt  USE/  3<iH  START  1 I *fr‘ , SEC  (Tl»1)  = E15.6/ 


1 

34* 

SMP  FMEOJE'JCY 

(Ns* IP) 

= lb/ 

? 

3 4 * 

MO.  OF  GAGES  ON 

TAPE  (NO,  AGE) 

= 16/ 

3 

JSH 

LOCATION.  I)  OF 

GAGES  = / ( 5 X , 

1014)) 

aonn  fokvut  (?3HnraPt'  tvP'ir  rs  only  i dimensional/ 

1 ?»H  ETT HER  wPX  09  MPY  'H'ST  RF  1) 

a?nn  format  (IQHOOaTa  FJ9  Q 4 GE  nC.  Ta,  15m,  LOCATION  ID  la/ 

1 14H  T I MF  INTt9VAL  ElS.ft.  IXM,  START  T I >IF  tl5.ft/ 

? 16h  NIJMRER  OK  TIm£S  15/  ?feM  MiMREP  OF  /.QROS  PER  POINT  la/ 

? PRw  OF  POINTS  PFW  RECORD  15/  1*,9A19) 

4 3uO  F OR  v A T (?t,H0PA91TY  ERROR  ON  DATA  TAPE) 
aaoo  FORMAT  (?5m00aTa  EXCEEDS  TaMLE  SPACE  ?I5', 
a 50  0 FORMAT  (32honu4iE9  JF  ACTIVE  GAGES  IS  *9  ‘.NG  3IU1 
9 7 0 0 F 09  vi  a T (1?H0E'fJ  JF  TAPE  214) 

4R00  FORMAT  (RPhOTAPE  10  MAS  PEEN  R F A D / 

1 19H  TIME  INTERVAL  = E15.S/1SH  NO.  OF  TI'/ES  = Ift) 

9900  FORMAT  (5H0  I = I3.4X.9HP  (PSI)  a / ( 1 X . 1 OE 1 3 . U ) ) 

F NO 
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*P"NOFC<  P*F$S 
•PURGE  P9F4S 


THIS  PACE  IS  BEST  QUALITY  PRACTICABLI 
raOJI  OOfrl  EURUISiiED  10  DOC - 


•aoOFILE  INPUT, INTI 

• ofck  press 

St'SROJT  I N»e  PRESS 

• CALL  C N 0 V A 
*C«l.L  C L 0 A 0 

• CALL  r>?LK1 
•CALL  C0M1 

•call  comp 

DIMENSION  PRTTd(<41  ) 

DIMENSION  TT-S(U,ai),PPTP(ia.an 

F 00 1 VALENCE  (PR  r (l . 1 . I ) .PFTft u . n ) . ( TTP( 1 . 1 , I) . TTft Cl . 1 ) J 

equivalence  (pRTi(i»i).pprrsci)j 

c 

IFfNCALL.GI.O)  Gu  TJ  «0/)0 
7.7z  1 .o/PTRIAL(l  ) 

GO  TO  (SO.PPft.rtiJO,  . !M L 0 4 0 

SO  IF  ( T Ime.GE.  TPRI  -ifc)  30  TO  1 no 
PPP  = ZZ*PP1 *(1 .0  - r I ME  *0  TT 1 ) 

IF  fPPP.LT. 0.0)  PPPsO.O 
GO  TO  400 

100  IFfTIME.CE. IT 0)  GO  TO 

PPP  = PPf)*(1.0  - T IME*0T  TO)  *•  A YN 
PpP=ZZ*PPP*EXP(-AZ*l IMF) 

GO  TO  400 
^On  PPPsO.O 

00  TO  400 

C 

PPO  00  ?a  0 JrJL.NTI-fe 

IF  (TIMF.LE.TTIJ))  3n  TO  ?f>0 
P4*  CONTINUE 
JL  = NTIME 
PPP  s Z7*PT(JL) 
liO  TO  400 
P*0  IL  S J 

PPP  s PT(J-l)  ♦ (THE  - TTCJ-l  ) )*  (PT(  J)  - pT(J-t))/ 

1 ITT(J)  - TT( J-l)  ) 

PPP  s 17* PPP 
C 

ti  0 0 PPP  = PPP  ♦ PM 

if  (<type.lt.»>)  go  io  poor 

P*  s PPP 

)F  (<TYPE.E0.10)  PX  s 0. 

Oil  S4I)  ISl.MMftaS 
Sun  Ph ( J > s P X 
GO  TO  4 0 0 0 
C 

4(>(>  IF  (<TYPF.GT.SJ  GJ  TO  QOO 

C pa*FlS. 

00  *•><>  1=1, UP* 

00  - m 0 J s 1 , ' P Y 
C T 'lTF  RPflL  A TE  ON  T1<E. 

PPP  s 0.0 

I F ( TTMp.L  I . T TP  (1  , J,  I ) ) GO  TO  ano 
•It  = J L T ( J , I ) 
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fHIS  PAGE  IS  BEST  QUALITY  PRACTICABI^ 

NTT  VC  : KTIM6(J»I)  FROM  COPY  FURNISHED  TO  DOC 

00  8?0  *=JL,nT1  •' t 

K K = < 

IF  (TlME.LE.TTPd,  J,  ID  GO  TO  890 
«i(i  Continue 

JLT(J.I)  = M I v t 
PPP  = PRTfNTIVt, J,  I) 

GO  Til  MhO 
*«0  JL  = 

PI  = PHT( JL-1. J.  1) 

T 1 = TTPf U 

PPP  : PI  ♦ ( T I vt  - Tl )*(P9T (JL, J,I)  - PI  )/(  T TPfJL.J*!)  - 1 1 ) 

JLT(J»I)  = JL 
880  P P T T ( J » II  = PPP 
C INTERPOLATE  SPATIALLY. 

K = o 

00  mao  Isl.NGT 
11  = I XI ( n 
OX  = 0 X 1 (I) 

DO  880  J=1,N8T 

IF  fNUSE(J,n.£<J.O)  GO  TO  8*0 
K = < ♦•  1 
JJ  = JYJ(J) 

0 Y = 0 Y u J ) 

PI  = PRTT(JJ.II)  ♦ OWPRTTfJJd.il)  - P9  TT ( J J,  I J ) ) 

P?  = PPTTf  JJ,  I i t- 1 J ♦ OV*(PRTT(JJd,II-d)  - PPTT'JJ,  II*t)  ) 

PPP  = PI  ♦ DX*lP£  - PI) 

P4U)  = POP*ZZ  t PS 
880  CONTINUE 
GO  TO  9000 
C 9EAMS. 

900  00  950  1 = 1 , NPS 
PPP  = 0.0 

IF  (TTME.LT.TTed,  [)  ) GO  TO  950 
JL  = J L 8 ( I ) 

NT  I ME  s LTIME(I) 

00  910  K = J L , TT T I M (£ 

= < 

IF  ( TIME.LE.TTetN, I) ) GO  TO  9P0 
910  continue 

JL  3 ( I ) = ‘ITI^E 
ppp  = p -v t ( » r i «t,  n 
GO  TO  oxn 
9P0  JL  = 8K 

PI  = P^cTPfJL-l  , I) 

T 1 = TT9( JL-1 . Ii 

PPP  = Pi  ♦ (Tl-it-Il  )*(PSTP(JL.T)  - PI  )/(TTi(JL»  IT  - M) 

JLBd)  = JL 
950  P<tTT*A(I)  = PPP 
* = 0 

00  940  IsI.NmaSS  ^ 

II  S ISPIII  i-‘ 

0*  = OSOO(I)  K 

PPP  = PPTTB(II)  * JX*(P9TT*5(TI*1  ) - PRTT4UI)) 

9<i  0 PP  ( I)  = PPP*ZZ  ♦PS 
GO  TO  9000 
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I 

I 


c 

c 

c 

l 'tort 
c 


110  0 


I ->0  0 


IPPO 
IP  *0 
IPSO 
1 100 


1 0(»0 


C 

1SM« 


1M>0 


1700 


1 SOO 

C 


l -ISO 
l qoo 


1 *)po 


T a P K OHM  ll . 


THIS  PACE  IS  BEST  QUALITY  PRACTICABLE 
OOFY  yURHISHED  TO  DOC  


r f fcij.o)  i;  i ro  isoo 

miPn^M  lCAO. 

t'.n  lioo  k:  JL  > ‘V  r l «•£ 

Tt  s OEL*FLOAT(k-U 
if  r t i a f . l t . 71 j gu  ro  i ?r-o 
CllM  I'U.'E 
JL  S 'IT  I -IE 

IF  ( -IT  IMF  .LF  ."4  Hue  ) ppp  r P ( 1 , MTT'aF  ) 

IF  ( '»T  I-xF.GT  .'■'A  i Jd  l ppo  s <J(  1 , 'jT  IME-VA  XO?) 

CO  TO  1X00 

JL  = x 

T 1 s T 1 • OF  L 

T1  s (TTMF-Tn/OtL 

if  ( jl-i .gt.maao^j  jo  to  iapo 

hi  s p ( i . JL-n 

rp  ( iL.or.  -ooP)  «•)  r.i  i?3o 

HP  2 p ( l . JL) 

GO  T)  IPSO 

PI  2 1(1.  JL-1-aax  JP) 
pp  2 )(1,  JL-M*xi)P) 

PHP  s Pi  ♦ T 1 * ( PP»P  1 ) 

PPP  z PPPoZZ  ♦ Ho 
IF  (aType.LT. s)  (,u  T!)  onoo 
P»  s PPP 

IF  (KTYPE.Fij.lij  Ha  = o. 

no  laOo  (sl.NMAbo 
Pn (t)  2 PX 
GO  TO  HO 00 
wOM«nMIF'.)ffF'  L 0 A ij  . 
no  i (>oo  xsJl.mi  Ag 
T I s 0EL*FL04T(K-l  ) 

IF  firAE.LE.TlJ  GO  T;j  1700 
COMT I jOF 
JL  2 'jTI"e 
GO  TO  1 SO  o 
JL  2 < 

T 1=  T 1 - DEL 
T 1 2 (TIMF-Tl)/OtL 
IF  (KTVPE.GT.S)  Go  TJ  P«oo 
PANELS. 

IF  IJL.LT.MP't)  SJ  TJ  ISOO 
DO  1 K 5 0 k G 2 1 . ' i G a U A 

IF  ( '.  T T •’£  . LE  . a x y ? ) PVT  TP  (KG)  = P(KC,mT1  -'£) 

IF  ( ! T I xE  ,G  T.  "4  J P«TTH(<S)  2 ‘.  ( <S,  VT  I v£ --HXGP  ) 

ClV  T I -nlE 
hi)  TO  P0('O 

on  isso  k c = l « .» o x 

IF  f JL-1  .GT.'-Axy*)  G'.j  TO  lHPo 

PI  2 P(KG.JL-I) 

TF  ( JL.GT.VAXOP)  G)  Tl  l«?0 
PP  s P(KG.JL) 
r-n  to  iss o 

HI  ; C(KG»JL*1  *;;M  ) ? > 


- 
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IPXP  P2  = T(KG.  JL-X’A*»2) 

1»*50  PPTT-t(KG)  : P!  ♦ Tt*(PP  - PI) 


Jins  PAGE  IS  BEST  QUAIT1^  rIUCIlCAB_4 
OUPK  JUFttUSHSl;  iO  VOQ  - 


JOiin  k s 0 

00  ? '-5 «■» 0 Isl.NGr 

IF  (NPX  .F  ).  1 ) GO  r ) 2050 

ii  = ixim 
ox  = oxt  (I) 

205"  00  ?mo  jn.vPT 

IF  ('iiJSE  ( J#  I ) .feO.OJ  50  TP  ??no 

K S K ♦ I 

IF  (VPY.FG.l)  Gu  TO  2110 
jj  S JY.KJ) 

0 Y s 0 Y 1 ( J ) 

PI  s PPTTS(JJ) 

PPP  : PI  ♦ py*(Pf#TT4(JJ*n  - PH 
P 4 t K ) a PPP*ZZ  ♦ PS 
GO  TO  2*01 

2100  IF  (J.GT.I)  GO  TO  2500 
PlsPPTTPCII) 

PPP  a (PI  ♦ OX«(Pi<rmiIM)  - P I ) ) * 7Z 
2200  P 4 ( K ) S PPP  ♦ Pa 
2 300  CONTINUE 

GO  TO  oonn 
: BEAMS. 


2400 

IF 

CJU.LT.NTIHE)  go 

TO  2500 

no 

2450  *G  = 1 , .yGSU* 

IF 

CyTIME.LE.mAX02) 

P«TTM*G) 

IF 

(MI'KE.GT.vi4Xu2) 

PVTTP ( KG ) 

2 450 

COMTIVIUE 

GO 

to  ?hon 

2501 

00 

2550  Kf,sl  ,.'jG5o*> 

IF 

( JL-1 .GT.^4X02) 

30  TO  2520 

PI 

= P(KG.JL-l) 

IF 

( Jt.GT .V4XD2)  GO 

TO  2530 

P2 

s P(KG.JL) 

GO 

TO  2550 

252  o 

PI 

s 0(kG. JL-1- *4X02) 

2530 

P? 

s 3(KG. JL-MAX02) 

3550  PWTTB(KG)  S Pt  ♦ T 1 * ( P2-P 1 ) 

2bon  no  p 7 oo  isum^abs 

II  = ISP(I) 

ox  = osoom 

PI  s PPTTH  (II) 

2700  PH(I)  = 7 Z * ( P l ♦ i)X*(PSTT«(TI»l  ) 


PI)  ) 


Pa 


on  oo 


rft'Tli* 
fe  ••  0 
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